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We present a systematic study of how vortices in superfluid films interact with the spatially vary- 
ing Gaussian curvature of the underlying substrate. The Gaussian curvature acts as a source for 
a geometric potential that attracts (repels) vortices towards regions of negative (positive) Gaus- 
sian curvature independently of the sign of their topological charge. Various experimental tests 
involving rotating superfluid films and vortex pinning are first discussed for films coating gently 
curved substrates that can be treated in perturbation theory from fiatness. An estimate of the 
experimental regimes of interest is obtained by comparing the strength of the geometrical forces to 
the vortex pinning induced by the varying thickness of the film which is in turn caused by capillary 
effects and gravity. We then present a non-perturbative technique based on conformal mappings 
that leads an exact solution for the geometric potential as well as the geometric correction to the 
interaction between vortices. The conformal mapping approach is illustrated by means of explicit 
calculations of the geometric effects encountered in the study of some strongly curved surfaces 
and by deriving universal bounds on their strength. 
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I. INTRODUCTION 

In superfluid helium, vortices form when the he- 
hum is rotated rapidly or when t here is turbulence 
(|Tillev and Tillevl . ll99d : IVinenl . ll969f) . Though such vor- 
tices are similar to the vortices that make up a vortex 
street behind the wings of an airplane or to the funnel 
clouds of tornadoes, they are only an Angstrom or two 
across ( Guvon et all . l200l[ ) . A more essential difference 
is that the vortices in a superfluid do not need a con- 
stant source of energy to survive. In fact, a vortex is 
long-lived because the strength of its flow is fixed by the 
quantization of angular momentum. Thus, the dissipa- 
tive mechanisms of a conventional fluid are absent. 

In this article, we focus on forces that the vortices ex- 
perience as a result of geometric constraints, with an em- 
phasis on those encountered in thin layers of liquid he- 
lium wetting a curved substrate with spatially varying 
Gaussian curvature. As a result of the broken trans- 
lational invariance of the underlying curved space, the 
energy of a single vortex with circulation quantum num- 
ber rii at position u^ includes both a divergent term 
and a position dependent self-energy, Es{ui), given by 
(IViteUi and Turn er. 2004) 



Es{u,) = -TrKn^Uaiu^), 



(1) 



C. Van der Waals Attraction on a Curved Surface 



where K 



Psti^ 



is the superfluid stiffness expressed in 



terms of the "^He atomic mass, m, and the superfluid 
mass density, Ps- The potential UQ{\ii) is obtained from 
solving a covariant Poisson equation with the Gaussian 
curvature, G(ui), acting as a source 



V2[/G(ui) = G(ui). 



(2) 



Vortices (and anti- vortices) are attracted (repelled) to 
regions of negative (positive) Gaussian curvature. These 
geometric interactions, while more exotic, are similar in 
origin to boundary-vortex interactions and can be suit- 
ably treated by the method of conformal mapping. 

Similar ideas naturally arise in a variety of soft- 
matter systems which have been confined in a thin 
layer wetting a curved substrate. The specific form 
of the resulting geometric interactions depends on the 
symmetry of the order parameter as well as on the 
shape of the substrate. Examples that have been 
studied both theoretically and experime ntally include 
colloidal crystals on curved interfa ces ("B ausch et al\ . 
l2003t iBowick et al\ . l2000l : IVitelU et al 200d), columnar 
phase s of block co-polymers (JHexemeiJ : ISantangelo et al\ . 
I2OO7I) as well as thin laye r s of nematic liquid crys- 
tals (iFernandez-Nieves et a?.l.l2007l: JPark and Lubenskvl . 



ll996HViteni and Nelsonl . 12004120061 ). Fueled by the drive 
towards technological applications b ased on the no t ion of 
self-assembly di r ected by geometry ( DeVries et al\ . 120071 : 
iDinsmore et all l2002t iNelsonl . 12002^ the study of these 
frustrated materials aims at predicting how the non- 
uniform distribution of curvature of the underlying sub- 
strates induces an inhomogeneous phase in the curved 
monolayer. An understanding of the resulting macro- 
scopic properties can be built from a mesoscopic descrip- 
tion cast in terms of the energetics of the topological de- 
fects which often exist even in the ground state and play 
a crucial role in determining how the material melts or 
ruptures. The advantage of this approach stems from the 
huge reduction in degrees of freedom achieved upon re- 
expressing the energy stored in the elastic field in terms of 
a few topological defects, rather than keeping track of the 
state of all the microscopic components, e.g. individual 
particle positions or molecular orientations. This step 
allows one to carry out efficient computa tional studies 
(JBowick et al\ . I2OOOI : iHexemer et al\ . l2007l ) and provides 
a suitable starting point for analytical work in the form 
of effective free energies derived from continuum elastic 
theory. 

Many of the mathematical techniques employed in this 
article to study vortices in curved superfluid films find 
application in the soft matter domain, in particular in 
those contexts where bond-ori entational order is impor- 
tant ( Vitelli and Turneii . l2004f) . In flat space both super- 
fluid and liquid crystal films can be describe d, as a first 
appr o ximat ion , by an XY modcl(Nclson an d Kosterlita . 
I1977I ). Both liquid crystal disclinations and vortices are 
modeled as a Coulomb gas of charged particles interact- 
ing logarithmically. However, the quantum nature of the 
problem considered in the present work introduces a fun- 
damental difference between these two classes of systems 



that is best illustrated by contrasting the angle of the 
liquid crystal director with the phase of the superfluid's 
collective wave function. The former represents the ori- 
entation of a vector (with both ends identified in the case 
of a nematic) that lives in the tangent space of the sur- 
face while the latter is a quantum mechanical object that 
transforms like a scalar since it is defined in an internal 
space. This subtle difference resurfaces upon consider- 
ing the distinct curved-space generalizations of the XY 
model that apply to each of these two systems. 

The free energy functional T^ to be minimized for 
the case of orientational order on a surface with points 
label ed by the coordinates u = (ui,M2) reads ( Davidl . 
1 19891) : 



K 



(3) 
where gap and g indicate the metric tensor and its de- 
terminant while r^Q (u) is a connection that compensates 
for the rotation of the 2D basis vectors Ea(u) (with re- 
spect to which 9{u) is measured) in the direction of Ua 



( Kamienl . l2002l ). Since the curl of the field f^gf u) is equal 
to the Gaussian curvature G(u) (|Davidl . ll989l ). the inte- 
grand in Eq. ([3]) never vanishes because ^^^(u) ^ dad 
on a surface with G(u) y^ 0. As the substrate becomes 
more curved, the resulting energy cost of geometric frus- 
tration can be lowered by generating disclination-dipoles 
in the ground state even in the absence of topological 
constraints. 

The connection rio,(u) is a geometric gauge field akin 
to the electromagnetic vector potential, with the Gaus- 
sian curvature playing the role of a magnetic field. If 
topological defects are present, they appear as monopoles 
in the singular part of daO{u). In analogy with elec- 
tromagnetic theory, their interaction with the Gaussian 
curvature arises mathematically from the cross-products 
between daO{u) and the geometry induced vector poten- 
tial ria(u), see Eq. ([3]). As a result of this interaction, 
disclinations in a liquid crystal are attracted to regions of 
the substrate whose curvatu re has the same sign as th e 
defect's topological charge (jPark and Lubenskvl . Il996f ). 
whereas vortices in a superfluid favor negatively curved 
regions independently of their sense of circulation. The 
anomalous coupling between vortices and Gaussian cur- 
vature introduced in Equations (P) and ^ originates 
from the distortion of the flow pattern by the protrusions 
and wrinkles of the surface. 

For a disclination with topological index rn (de- 
fined by the amount 9 increases along a path enclosing 
the defect's core) the geo metric potential Ey{ui) reads 
(IVitem and Turneii 120041) 



S,(u,) = 27ri^n, (1-^) Ug{u,), 



(4) 



where K is the elastic stiffness and Uq is the same po- 
tential defined in Eq. ^. Note that the anomalous 
coupling also contributes to determine the energetics of 



liquid crystal disclinations but, in this case, the gauge 
coupling, which is linear in n, is stronger for small n. 

To understand the physical and mathematical origin 
of these distinct coupling mechanisms, note that in the 
ground state of a ^He film, the phase 0{u.) can be con- 
stant throughout the surface so that the corresponding 
energy vanishes. In a system with geometric frustration, 
the gauge coupling between defects and the underlying 
curvature is mediated by the deformed ground state tex- 
ture that exists in the liquid crystal layer prior to the 
introduction of the defects simply as a result of geomet- 
rical constraints. Once a defect is introduced it interacts 
with these preexisting elastic deformations. Unlike the 
case of orientational order considered previously, no ge- 
ometric frustration exists in the superfluid film. The su- 
perfluid free energy J^g to be minimized is a simple scalar 
generalization of the familiar flat space counterpart 



K 



.f. = y / d'u^g''''do,0{u)di3eiu) 



(5) 



The crucial point is that no connection fla (u) is necessary 
to write down the covariant derivative for this simpler 
case of a scalar order parameter. Therefore the ground 
state is given by 9{x.) equal to a constant. There is no 
preexisting texture for a vortex to interact with, and so 
another mechanism is required to explain the coupling of 
vortices to geometry. 

In the following sections, we will employ the method of 
conformal mapping to demonstrate that when an isolated 
vortex is placed on a curved surface it feels a force as 
if there were a smeared out topological "image charge," 
jointly proportional to the vortex's own circulation and 
the Gaussian curvature across the substrate. Such an 
imaginary topological charge distribution produces a real 
force analogous to the force on an electrostatic charge due 
to its mirror image in a conducting surface. 

The method of conformal mapping may seem, prima 
facie, a surprising route to derive a coupling between 
vortices and geometry, because the free energy of Eq. 
([5|) is invariant under conformal transformations that in- 
troduce a non-uniform compression of the surface while 
keeping local angles unchanged. This invariance property 
at first seems to rule out the possibility of a geometrical 
interaction! This apparent contradiction can be seen by 
choosing a special set of (isothermal) coordinates that 
can always bring the two dimensional m etric tensor i n 
the diagonal form ga/3(u) = e2'^("U„/3(u) (|Davidl . fl989l ). 
The result of this step is to eliminate the geometry depen- 
dence from the free energy of Eq. (O since the product 
g°'^{u)y/g = Sap{u) and J^s reduces to its counterpart 
for a planar surface, where there is of course no geome- 
try dependence. 

An interaction between vortices and geometry violates 
this conformal symmetry of the free energy from which 
it emerges, but in fact the conformal symmetry is not 
an exact symmetry when vortices are present. Analo- 
gous subtleties frequently arise in the study of fields that 
fluctuate thermally or quantum mechanically, due to the 



occurrence of a cut-off length scale below which fluctu- 
ations cannot occur. A conformal mapping is a strange 
type of symmetry that stretches lengths and thus does 
not preserve the microscopic structure of a system. At 
finite temperatures, the discreteness o f a system, such 
as a thermally fluctuating membrane ( Polvakovl Il98lh 
which is actually made up of a network of molecules, can 
have an important effect because the fluctuations excite 
modes with microscopic wavelengths. This produces vi- 
olations of the conformal symmetry at every point of the 
surface. In a superfluid at zero temperature, however, 
short wavelengths not describable by the continuum free 
energy JF^ are excited only in the cores of vortices. Ob- 
taining a finite value for the energy necessitates the re- 
moval of vortex cores of a certain fixed radius in the local 
tangent plane, so a conformal mapping is not a symmetry 
in the neighborhood of a vortex. However the amount by 
which this symmetry fails can be calculated in a simple 
form (intriguingly independent of the microscopic model 
of the cores) in terms of the rescaling function w(u) evalu- 
ated at the locations of the vortices, where the symmetry 
fails. Rather than ruling out the possiblity of a geometric 
interaction, a realistic treatment of conformal mapping 
becomes a powerful mathematical tool for deriving these 
interactions, a technique which is relevant especially to 
other branches of theoretical physics su ch as the study o f 
scattering amplitudes in string theory (jPolvakovlll987l ). 

While the free energy, J^s, of the curved superfluid layer 
in Eq. ^ does not exhibit a geometric gauge field, rotat- 
ing the sample at a constant angular velocity leads to an 
energy of the same form as Eq. Q . The resulting forces 
exerted on the vortices compete with the geometric in- 
teractions to determine the equilibrium configurations of 
an arrangement of topological defects. This simple idea 
is behind some of the experimental suggestions put for- 
ward in this article to map out the geometric potential by 
progressively increasing the rotational speed while mon- 
itoring the equilibrium position of say a single vortex on 
a heliu m coated surface shaped like the bottom of a wine 
bottle ( Voll et all . 120061 ). Since the position dependence 
of the force induced by the rotation is easily calculated, 
one can read off the geometric interaction by simply as- 
suming force balance. 

The theory of curved helium films also helps build in- 
tuition for the more general case of vortex lines confined 
in a bounded three dimensional region such as the cavity 
shown in cross-section in Fig. [TJ'V. The vortex, drawn as 
a bold black line, can be pinned by the constriction of 
the container. The classic problem of understanding the 
interaction of the vortex wit h itself and with the bump 
as the superfluid flows past ( Schwara . Il98ll ). is of cru- 
cial importance in elucidating how vortices can be pro- 
duced when a superfluid starts rotating despite the ab- 
sence of any friction. A possible mechanism, known as 
the "vortex mill" , assumes that vortex rings break off 
a pinned vortex line, while the pinned vortex remain s 
in place (JGlaberson and Donnellvl . Il966t ISchward . Il990t ) . 



The common route to studying vortex dynamics in three- 




FIG. 1 Cross sections of two regions in which one can study 
vortex energetics. A) A region with nonparaUel boundaries. 
The vortex is pushed to the right. This tendency can be 
interpreted either in terms of a drive toward a shorter length 
or as the local induction force due to the curvature in the 
vortex enforced by the boundaries. B) A cross-section of a 
constant thickness layer of helium bounded above by air and 
below by the substrate. The vortices keep a constant length 
and remain straight while moving around. Hence there is no 
local induction force/thickness-variation force to overwhelm 
the geometrical forces that we focus on. 



dimensional geometries is the "local induction approxi- 
mation" which assumes that each element of a vortex 
experiences a force determined only by its local radius of 
curvature. This simplifying assumption omits any long 
range forces experienced by the vortices as they interact 
with the boundaries (or among themselves). In the op- 
posite limit of films with uniform thickness (which can be 
thought of as special types of bulk supcrfluid regions with 
two parallel boundaries, as in Fig. [Tj3), all the forces ex- 
erted on the vortices are long-range. This is the regime 
of interest to our investigation. 

This article is organized in two tracks. The first, com- 
prising sections HIl- IIVl is phenomenological in nature and 
emphasizes intuitive analogies between the (non-linear) 
geometric forces and conventional electrostatics, simple 
illustrations of the main results and experimental ideas. 
The second track, sections IVl- IVIl is more technical and 
presents a unified derivation of the geometric potential 
by the method of conformal mapping and its application 
to the study of complex surface morphologies. 

The first track starts with a review of superfluid dy- 
namics that can be used to relate the anomalous coupling 
to hydrodynamic lift. In Sec. III.Al the geometrical force 
is evaluated, using a mapping between the geometric po- 
tential studied here and the familiar Newton's theorem 
that allows an efficient calculation of the gravitational 
field for a spherically symmetric mass distribution. An 
intriguing consequence of this analogy is that vortices on 
saddle surfaces can be trapped in regions of negative cur- 
vature leading to geometrically confined persistent cur- 
rents as discussed in section III.BI Section III.CI relates 
this observation to Earnshaw's theorem from electrostat- 
ics. Upon heating and subsequently cooling a curved 
superfluid fllm, some of the thermally generated defects 
can remain trapped in metastable states located at the 
saddles of the substrate. The existence of such geometry 
induced vortex hysteresis is conjectured in section III.DI 
In section IIII.Al we derive the forces experienced by vor- 



tices when the vessel containing the superfluid layer is 
rotated around the axis of symmetry of a curved sur- 
face shaped as a Gaussian bump. The dependence of 
single and multiple defect-configurations on different an- 
gular speeds and aspect ratios of the bump is studied in 
sections IIII.BI and IIII.CI The Abrikosov lattice of vor- 
tices on a curved surface is discussed in section IIII.DI In 
realistic experimental situations the thickness of the su- 
perfluid layer will not be uniform and additional forces 
will drive vortices towards thinner regions of the sample. 
The strength of these forces is assessed in section ITVl and 
related to spatial variations of the film thickness due to 
gravity and surface tension. A short discussion of choice 
of parameters for the proposed rotation experiments fol- 
lows in section IIV.BI The relevance of our discussion 
to experiments performed in bounded three dimensional 
samples is addressed in section IIV.DI 

The second track starts with a general derivation of the 
geometric potential by the method of conformal mapping 
in section IV. Al The computational efficiency of this ap- 
proach is illustrated in section FV.BI where the geometric 
potential of a vortex is evaluated on an Enneper disk, a 
strongly deformed minimal surface. We show that chang- 
ing the geometry of the substrate has interesting effects 
not only on the one-body geometric potential but also 
on the two-body interaction between vortices. In section 
IV. CI we use conformal methods to show how a periodic 
lattice of bumps can cause the vortex interaction to be- 
come anisotropic. In section IV. Dl we demonstrate that 
the quantization of circulation leads to an extremely long- 
range force on an elongated surface with the topology of 
a sphere. The interaction energy is no longer logarithmic, 
but now grows linearly with the distance between the two 
vortices. As we demonstrate, the whole notion of split- 
ting the energy in a one body geometric potential and a 
vortex-vortex interaction is subject to ambiguities on de- 
formed spheres. Section [V.EI provides some guidance on 
how to perform calculations in this context by choosing 
a convenient Green's function among the several avail- 
able. Finally, in section IVII we present a discussion of 
some general upper bounds which constrain the strength 
of geometric forces. The conclusion serves as a concise 
summary and contains a table designed to locate at a 
glance our main results throughout the article including 
the more technical points relegated to appendices but 
useful to perform calculations. 



II. FLUID DYNAMICS AND VORTEX-CURVATURE 
INTERACTIONS 

We start by writing down the collective wave function 
of the superfluid as 



vi/(u) = J^iM e'«(") 



7714 



(6) 



where u = {ui,U2} is a set of curvilinear coordinates 
for the surface, m is the mass of a *He atom and ps 



is the superfluid mass density that we shall assume to 
be constant in what follows. To obtain the superfluid 
current we can use the standard expression Jq(u) = 
2^ {^da^* - **9q*) showing that the superfluid ve- 
locity is given by 



z;a(u) = da9{u). 

1714 



(7) 



The circulation along a path C enclosing a vortex is given 
by 



du°'i 



(8) 



where the quantmn of circulation, k = — , is equal to 

9.98 10~* m^ s~'^ and the integer n is the topological 
index of the vortex. The free energy can be cast in the 
form 



1 h^ 



(fu^g'^Pdo.edpe , 



(9) 



where g"^ is the (inverse) metric tensor describing the 
surface on which the superfluid layer lies and g is its 
determinant. We will often use the superfluid stiffness 



K 



Psh^ 



m^ 



(10) 



This expression for the free-energy can be parameter- 
ized in terms of the vortex positions once the seemingly 
divergent kinetic energy near a vortex core is correctly ac- 
counted for. As is well known, the radius-independence 
of the circulation about a vortex implies that the veloc- 
ity in its proximity is given by — —, which leads to a 
logarithmic divergence in ([5]). The energy stored in an 
annulus of internal radius r^ and outer radius Tout reads 



J^out 



(11) 



which diverges as rin ^ 0. A physical trait of superfluid 
helium prevents this from happening: it cannot sustain 
speeds which are greater than Uc, the critical velocity. 
Thus the superfluidity is destroyed below a core radius 
of a ~ _2_ This breakdown may be modeled by excis- 
ing a disk of radius a around each vortex and by adding 
a constant core energy tc to account for the energy as- 
sociated with the disruption of the superfluidity in the 
core. 

Starting on the flat plane, the interaction of two vor- 
tices can now be determined. Superimposing the fields 
of the two vortices and integrating the cross-term in the 
kinetic energy of Eq. © leads to a Coulomb-like inter- 
action, Vij ~ ~2TTKninj In — — ^ in addition to vortex 
self-energies. In deducing the force between the vortices 
from this expression, it is useful to assume that a does not 
vary significantly with position. The justification for this 
simplification is that the background flow due to other 



vortices only gives a fractionally small correction to the - 
flow near each vortex, and therefore barely affects where 
the critical velocity is attained. 

For the more complicated cas e of a curved surface, with 
a very distant boundary (see ( Vitelli and Nelsonl . 12004 ) 
for the discussion of effects du e to a boundary at a finite 
distance), we found in Ref. ("Vitelli and Turnen . 12004 ) 
that the energy including both single-particle and two- 
particle interactions is. 



i<j 



(12) 

apart from a position-independent term (given for a dis- 
tant circular boundary of radius R by T^i^ifiiYhi— + 
N^, with N the total number of vortices and ec the core 
energy of one of them). The pair potential Vij = r(ui, Uj) 
is expressed in terms of F, the Green's function of the co- 
variant Laplacian defined by: 



V2f(u,v) = -<5,(u,v) 



(13) 



Note that the covariant delta function Sc includes a fac- 
tor of -y= so that its integral with respect to the "proper 

area" ^duidu2 is normalized. Equation (J13p determines 
the Green's function up to a constant provided that we 
assume additionally that the Green's function is symmet- 
ric between its two arguments. The constant is fixed by 
assuming that at large separations, the Green's function 
approaches the Green's function of an undeformed plane. 
This expression shows that vortices behave like electro- 
static particles, with charges given by 27rni and coupling 
constant K. 

The single-particle potential [/g(u) is the "geometric 
potential" defined in Eq. ([2|). This potential entails a 
repulsion VttKUg of vortices of either sign from positive 
curvature and an attraction to negative curvature. The 
following analogy with boundary interactions and image 
charges is useful. The flow-field is modified by having to 
conform to the curvature, leading to an image charge of 
the vortex which is spread continuously over the surface, 
with density — 2^G(u). Just like the image of a vortex in a 
circular boundary has an equal and opposite circulation, 
the continuous image of the vortex in the curvature has 
a charge density proportional to the number of quanta n 
in the vortex. 

This point of view may be connected to fluid mechan- 
ics by analyzing the streamlines on the bump and in the 
presence of a circular boundary as illustrated in Fig. [2l 
Streamlines are tangent to the direction of flow and their 
density is proportional to the local speed. Only an in- 
compressible velocity field (div v = 0) may be described 
by streamlines, since incompressibility ensures that any 
closed curve has an equal number of streamlines enter- 
ing and exiting. This condition is satisfied for superfluids 
(far below the critical speed) since minimizing ^ leads 
to V • V6' = 0, or div V = according to Eq. 0. Thus 
the flow is both irrotational (the circulation around any 



curve, not enclosing a vortex, is 0) and incompressible: 

(14) 



div V = 
curl V = 0. 



The former relation implies that we may write 
V = curl xn; 



(15) 



so that the streamlines are equally spaced level curves 
of X- (For example, around a vortex, the radii of the 
successive streamlines forms a geometric sequence, r(l — 
e) where e sets the ratio between flowline density and 
speed.) 

Now consider, as an illustration, the flow field on the 
slope of the bump represented in Fig. [3 The curves 
must spread out to go over the bump, leading to a lower 
velocity above the vortex. By Bernoulli's principle (true 
for irrotational flows), this creates a high pressure that 
pushes the vortex away from the bump. Note, however, 
that the actual motion of a vortex is more subtle: Al- 
though the gradient of the energy points away from the 
bump, a vortex (disregarding friction) always moves at 
right angles to the gradient of the energy. Thus a vortex 
in the absence of drag forces actually circles around the 
bump, in the same direction as the fluid flows around the 
vorte x. Dissipation converts th e motion into an outward 
spiral! Ambegaokar et a/.l . Il978l ). We will not study the 
dynamics. 

A convenient mathematical formulation of the problem 
of determining the flow pattern of a collection of vortices 
is obtained by introducing the scalar function x(u) which 
satisfies 



V^x 



(u) = - ^ 27rni(5c(u, Ui) ee -cr(u). (16) 



The sum can be described as a singular distribution of 
surface charge. This relation follows from the circulation 
condition, 2'Kni — § V0-dl, which can be rewritten as the 
integral of the flux of Vx through the boundary, § Vx • 
hdl, by using Eq. P^ . In analogy with Gauss's law, 
there must therefore be delta-function sources for x ^t 
the locations of the vortices, as described by Eq. P^ . 
Solving Eq. (fTB|) in terms of the Green's function gives: 



hrii 



,(u) = ^^r(u,u.) 



(17) 



The flow due to a given vortex is proportional to its 
"charge" 27™^. The energy as a function of the posi- 
tions of the vortices, Eq. (H^), can now be derived by 
integrating the kinetic energy in the flow determined by 
Eq. (fT6|) for each placement of the vortices. 

We will begin by discussing applications of Eq. P^ . 
saving its derivation until later. Interestingly, the energy 
of the vortices can be described by a differential equation 
analogous to Eq. (fTB|) for the flow. We choose one vor- 
tex and flx the positions of all the others. We take the 



Laplacian of Eq. (|T2|) with respect to the position of the 
chosen vortex and use Eq. ^ and Eq. p^ . The energy 
as a function of the chosen vortex £'(ui) satisfies: 



-^1^ --'"■) -?=<-) 



(18) 



The notation ai stands for the delta function charge dis- 
tributions of all the vortices with the exception of the i**^, 
so that ai{u) = ^ ,j27mj(5(u — Xj). The self-charge 
term which we have had to remove (so that Ui can be 
substituted in place of u as in Eq. p5|) ) is replaced here 
by a spread-out charge proportional to the removed term. 



A. Anomalous force on rotationally symmetric surfaces 

On an azimuthally symmetric surface the force on a 
defect can be found by exploiting Gauss's law. This is 
analogous to the familiar "Newton's Shell" theorem that 
predicts the gravitational field at the surface of a sphere 
surrounding a spherically symmetric mass distribution by 
concentrating all the enclosed mass at the center. 

Points on an azimuthally symmetric two dimensional 
surface embedded in three dimensional Euclidean space 
are specified by a three dimensional vector R(r, 0) given 
by 



R(r,, 




(19) 



where r and are plane polar coordinates in the xy plane 
of Fig. [21 and h{r) is the height as a function of radius; 

2 

e.g. h{r) — Hqc ^''3 for the Gaussian bump with height 
ho and spatial extent '^ tq. It is useful to characterize 
the deviation of the bump from a plane in terms of a 
dimensionless aspect ratio 



ro 



(20) 



The metric tensor, gap, is diagonal for this choice of co- 
ordinates. In general, g^^ = r^, g^r = 1 + h'(r)^ , and for 
the Gaussian bump we have 



ffa/3 



1 
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a^ -T exp 

'n 









(21) 



Note that the g^^ entry is equal to the flat space result r^ 
in polar coordinates while g^r is modified in a way that 
depends on a but tends to the plane result grr = 1 for 
both small and large r. 

The Gaussian curvature for the bump is readily found 
from the eigenyalues of the second fundamental form 
( Dubrovin et aLl , Il992l ): e.g., for the Gaussian bump. 
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(22) 





FIG. 2 Left, the flow around a vortex situated on the side of a Gaussian bump, calculated with the methods described in this 
paper. The low density of flow lines above the vortex indicates a lower velocity and thus a higher pressure, leading to the 
repulsion represented in Eq. p2p . — V(— ttKii'^Ug). Right, the analogous flow around a vortex in a disk with a solid boundary. 
The attraction to the boundary is also seen to result from high speeds, since the flow lines are compressed in between the vortex 
and the boundary. 




FIG. 3 (a) A bumpy surface shaped as a Gaussian, (b) Top 
view of (a) showing a schematic representation of the positive 
and negative intrinsic curvature as a non-uniform background 
"charge" distribution that switches sign at r = ro. The vary- 
ing density of -I- and - signs tries to mimic the changing cur- 
vature of the bump. 



Note that a controls the overall magnitude of G{r) and 
that G{r) changes sign at r = ro (see Fig. [SJj). The 
integrated Gaussian curvature inside a cup of radius r 
centered on the bump vanishes as r — > oo. The pos- 
itive Gaussian curvature enclosed within the radius tq 
(see Fig. [3]) approaches 2tt for a ^ 1, half the integrated 
Gaussian curvature of a sphere. We shall show below 
that there is always more positive than negative curva- 
ture within any given radius, for an azimuthally sym- 
metric surface. It will follow that the force on a vortex 
is repulsive at any distance. 

In general an individual vortex of index rii confined on 
a curved surface at position u^ feels a geometric interac- 
tion described by the energy 



E{u,) = -TTKnfUG{vii). 



(23) 



For an azimuthally symmetric surface such as the bump 
represented in Fig. [31 we can derive Newton's theo- 



rem as follows. Define E — —VUg so that the co- 
variant radial component of 'Ei is Ej. — —drUa- Then 
—V^Ug — div E — -j=d,^/gg^^Er, and if we integrate 
both sides of Eq. ^ out to r. 



2n^g^'-Er = 



^drdc^Gir) 



(24) 



so that Er has a simple expression in terms of the net 
Gaussian curvature at a radius less than r. Now E^ is 
the "covariant component" of the geometrical "electric 
field" , not the actual field, which would be obtained by 
differentiating with respect to arclength rather than the 
projected coordinate r. Therefore the magnitude of E is 
^E=, which, rephrasing Eq. (|M)) . obeys this generalized 

version of Newton's theorem: 

The magnitude of E is —-^ — times the integrated Gaus- 
sian curvature. 

(Recall that g^'^ ~ -^ and g = det gap = rgrr-) Note 
that the force on the vortex is F = —VE = —nKnfEi ac- 
cording to Eq. (pS)) , so that the geometrical force is pro- 
portional to the integrated Gaussian curvature divided 
by the distance of the vortex from the axis of symmetry 
of the surface; the expression for the force on the vortex 
obtained in the next section by integrating the Gaussian 
curvature is 



_Kti 1 , 



(25) 



if rii = ±1. Note that this force is always repulsive since 
the integrated curvature is positive. 

The geometric potential can now be expressed explic- 
itly by integrating E,. = — 9,[/g with the aid of Eq. (I24p : 



UG{r) = - 



dr 



1 



■ exp 



-1 



(26) 




FIG. 4 Plot of the interaction energy E{r) = —TrKUcir) 
between a singly quantized vortex and a Gaussian bump with 
a = 1. The energy is measured in units of K and the radius 
is measured in units of vq. Note that the force points away 
from the bump and has its maximum strength near ro. 




FIG. 5 An azimuthally symmetric substrate and its down- 
ward projection on a flat plane. The shaded strip surround- 
ing P is more stretched than the one surrounding Q despite 
their projections onto the plane having the same area. As 
discussed in the text it follows that the energy stored in the 
field will be lower if the center of the vortex is located at Q 
rather than P. 



The resulting potential UGir) vanishes at infinity. Its 
range and strength are given respectively by the linear 
size of the bump and its aspect ratio squared (see Fig. 

a 

We now summarize an intuitive argument that ex- 
plains why the energy of a vortex on top of the bump 
i s greater than the energy of a vor tex that is far away 
(JHalpcrin': Vitelli and Turner , 12004 ). Fig. [3 reproduced 
from (| Vitelli and Turnerl . 12004) . focuses on a rotationally 
symmetric bump coated by a helium film (of a constant 
thickness) . We can estimate the energy of a vortex on top 
of the bump by comparing the situation to a vortex on a 
plane, illustrated vertically below the bump. Rotational 
symmetry implies that the superfluid phase is given by 
9 = (f), the azimuthal angle. The velocity depends on the 
rate of change on the phase according to Eq. ([7]) , so since 
an infinitesimal arc of the circle concentric with the top 
of the bump has size rd(f>, the velocity is — . Here r is 
the radius of the circle measured horizontally to the axis 
of the bump. This calculation shows that the velocity, 
and thus the energy density, are the same at any point 
on the bump and its projection into the plane. However, 
the energy contained in the tilted annulus on the bump 
stretching from r to r + dr is greater than the energy in 
the annulus directly below it because, though the energy 
density is the same, the annulus's area is greater. Hence 
a vortex on a bump has a greater energy than a vortex 
in a plane, whether it is the vortex at P' in the projec- 
tion plane, or the vortex at Q which is very far from the 
bump. This reasoning indicates a repulsive force, since 
the vortex lowers its energy by moving away from the 
bump. 

The intuitive argument suggests that the Gaussian cur- 
vature should appear in the force law, as in Eq. ([2]); 
in fact, it is a widely known fact that the sign of the 
Gaussian curvature of a surface determines how fast the 
circumference of a circle on the surface increases, rela- 



tive to the circumference of a circle on the plane, as a 
function of the radius. Of course, for less symmetric sur- 
faces, comparing the energy on a curved surface to that 
on a flat reference plane directly below it will be more 
complicated, since symmetry and the constant circula- 
tion constraint do not force the energy densities to be 
equal. Thus, a simple vertical projection will not set up 
a monotonic relation between energies. The conformal 
mapping technique which we use in Sec. |V]is a variation 
on the idea of comparing a "target" substrate to a sim- 
ple "reference" surface which is in principle applicable to 
arbitrary surfaces, and furthermore not only allows one 
to compare energies, but also to calculate them quantita- 
tively. The technique can also be used to give a concise 
derivation of Eq. ([2|) . 

Such an intuitive argument applies only for az- 
imuthally symmetric surfaces. For less symmetric sur- 
faces, comparing the energy on a curved surface to that 
on a flat reference plane directly below it will be more 
complicated, since symmetry and the constant circula- 
tion constraint do not force energy densities to be equal 
at corresponding postions. Thus, a simple vertical pro- 
jection will not set up a monotonic relation between en- 
ergies. The conformal mapping technique which we use 
in Section |V] is a variation on the idea of comparing a 
"target" substrate to a simple "reference" surface which 
is in principle applicable to arbitrary surfaces, and fur- 
thermore not only allows one to compare energies, but 
also to calculate them quantitatively. The technique can 
also be used to give a concise derivation of Eq. ([2]) . 



B. Vortex-trapping surfaces 

In order to illustrate the consequences of the curvature- 
induced interaction for different surface morphologies, we 




FIG. 6 Plot of vortex trapping surface. 



study a "Gaussian saddle" surface (suggested to us by 
Stuart Trugman) for which the geometric potential has 
its absolute minimum at the origin. 

First, we show that an alternative design for a vor- 
tex trap geometry fails because of the long-range na- 
ture of the curvature-induced interaction. Fig. [3] shows 
that bumps have negative curvature on their flanks; it 
might seem possible that a well-chosen bump would have 
enough negative curvature to hold a vortex. However, 
a vortex cannot be held by nearby negative curvature 
alone; it also feels the positive curvature at the center of 
symmetry because, according to Gauss's law applied to 
the azimuthally symmetric region, the force is due to the 
net curvature contained in any circle concentric with the 
top of the bump. This curvature is given generally by 



G 



1 



--dr 



1 



(27) 



and the net curvature within radius r^ is thus 

Vff dr dcj, G{r) = 27r(l - \ ) 

= 27r(l- cos 6l[r.„]) (28) 

where 6 is the angle between the surface at the location 
of the vortex and the horizontal plane. This formula also 
describes the cone angle of a cone tangent to the sur- 
face at radius r^; it can also be derived from the Gauss- 
Bonnet theorem which implies that the net curvature of 
a curved region depends only on the boundary of the re- 
gion and how it is embedded in a small strip containing 
it; thus the net curvature of the cone (concentrated at 
the sharp point of the cone) is the same as the net cur- 
vature of the bump which it is tangent to. Since this 
curvature is always positive, the defect is always repelled 



from the top of an azimuthally symmetric bump^. To 
find a way to trap a vortex, one must therefore investi- 
gate some non-symmetric surfaces. The curvature-defect 
interaction energy on a generic surface is mediated by the 
Green's function of the surface, Eq. ((13]), as can be seen 
by solving Eq. ^■. 



E[vl) = Kit fd^u T{u, v) G(v), 



(29) 



for a singly quantized vortex. One such surface, which 
we will treat perturbatively in the amount of deformation 
from flatness (Section IV. Bl treats a different confining 
surface exactly) is the Gaussian saddle represented in 
Fig. [5] and described by the height function 



h\{x,y) = — (x^ - Aj/2) e 



''o 



(30) 



where the exponential factor was i ncluded to make the 
surface flat away from the saddle ( Trugmanh . Here, A 
is a parameter which we later vary to illustrate the non- 
locality of the interaction. The leading order contribution 
to the curvature-defect interaction (for a << 1) is of the 
same order a^ as the curvature corrections to the defect- 
defect interaction (calculated in Appendix[A| . In fact F is 
multiplied in Eq. pS)) by the Gaussian curvature G'(u'), 
of order a^ . Thus, for a single defect, it is sufficient to 
use the flat space Green's function T fiat 



^flat{x 



,y-y') = -7^ log y/{x-x')^ + {y-y')^) 
Zn 

(31) 
for calculating the defect-curvature interaction. Further- 
more, the Gaussian curvature that acts as the source of 
the geometric potential can be calculated from the usual 
second-order approximation: 



G\ix,y) 



dx^ dy"^ 



\dxdy) 



(32) 



This function is plotted in Fig. [71 and its sign is repre- 
sented in the middle frame of FigurefTUl The graph of the 
vortex-curvature interaction energy for this surface. Fig. 
m shows that a vortex is indeed confined at the center of 
the saddle; the energy graphed in this figure is given by 



E\{x, y)Ki K-K I dx' dy' Tfiatix - x',y- y') Gx{x' , y'), 



(33) 



for A = 1. For realistic film thicknesses and a of order 
unity (see Sec. IIV|) . the depth of the well is about 50 



^ More specifically, all bumps with azimuthally symmetric embed- 
dings repel vortices from their tops. The negative curvature 
cones discussed in Sec. I VII have an internal azimuthal symmetry 
but their three dimensional embeddings are not symmetric. 
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FIG. 7 Plot of the curvature of the A = 1 vortex trapping 
surface. 




FIG. 8 Plot of the geometric potential for the A 
trapping surface. 



1 vortex 



Kelvin! We have found that the energy associated with 
a vortex at the origin is less than for any other position. 
Of course, the configuration with a vortex at the origin 
cannot beat the configuration with no vortices at all! The 
latter has zero kinetic energy; when the vortex is at the 
origin, the energy is positive provided that Eq. p^ is 
supplemented by the position-independent contribution 
ttK In — . This term is always necessary for comparing 
configurations with different numbers of defects, as when 
one studies the format ion of a vortex lattice at i ncreasing 
rotational frequencies! Campbell and Zifj . [19791 1. 



C. Negative curvature which does not trap 

In this section, we shall discuss what happens when 
the parameter A of the saddle surface is increased; Fig. 
ini illustrates such a surface corresponding to A = 17. To 
give a hint of what causes the equilibrium to change its 




FIG. 9 A saddle surface with A — 17; this parameter value is 
just large enough to destabilize a vortex at the center. 



character. Fig. [TOl shows the sign of the curvature for the 
Gaussian bump, the saddle with A = 1, and the saddle 
with A = 17. 

In the graph of the defect-curvature interaction energy 
with A = 17, one notices that the origin is an unstable 
equilibrium position for the vortex. We will derive the 
exact value of A where this instability first occurs below. 
However, symmetry considerations alone show that the 
origin is a stable equilibrium point when A = 1, as Fig. 
[5] shows. One might be tempted to argue from Newton's 
theorem that a vortex at a small enough radius r is always 
attracted to the origin by the negative curvature at radii 
smaller than r. However, the asymmetry of the saddle 
surfaces invalidates Newton's theorem and positive cur- 
vature more distant from the origin than the vortex might 
be able to push the vortex toward infinity. This does not 
occur for the saddle surface with A = 1; although the 
rotational symmetry needed for Newton's theorem is ab- 
sent, the surface does have order four symmetry, under a 
90 degree rotation combined with the isometry z — > — z. 

Upon expanding the defect-curvature interaction en- 
ergy about the origin, we obtain 



Eq + ax + by + ex + 2dxy + ey + 



(34) 



This energy must be invariant under the symmetries of 
the surface (without a sign change). Order two symme- 
try implies that the linear terms vanish, so the center 
point is an equilibrium. The order four symmetry (ap- 
parent in Fig. HOB ) implies that it is either a maximum 
or a minimum (a quadratic function with a saddle point 
has only 180 degree symmetry). In more detail, 90° rota- 
tional symmetry, given hy x —f y, y— >— x, implies that 
c — e,d = 0. Since the Laplacian of E at the origin is 
proportional to minus the local curvature, 2c = 2e — c+e 
is positive, so the origin is a local minimum. Without the 
order four symmetry the negative 
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(A) 



(B) 



(C) 



FIG. 10 Plots of the sign of the curvature, with white for positive curvature. (A) is for the Gaussian bump, and (B) and (C) 
are for the saddle surfaces with A = 1 and 17 respecitvely. Because of the lack of symmetry in the third figure, the center point 
becomes a saddle point of the energy; the vortex is pushed away by the strong positive curvature in the ellipsoidal regions at 
positive and negative y. 



curvature only ens ures that c+e > 0. Earnshaw ' s theo- 
rem o f electrostatics (JEarnshawl [l842l : iJeansl . Il927l : IScottl . 
[195 9) states that an electric charge cannot have a stable 
equilibrium at a point where the charge density is zero 
or has the same sign as the charge. The charge cannot 
be confined by electric fields produced by electrostatic 
charge distributions in a surrounding apparatus. (This 
theorem motivated the design of magnetic and electro- 
dynamic traps for trapping charged particles in plasma 
physics and atomic physics.) The argument provided 
here can be generalized to give the following converse rule 
based on discrete symmetries (whereas Newton's theorem 
applies only for continuous azimuthal symmetry): 

If P is a symmetry point of a charge distribution with 
rotation angle — , and rn > 3, and the charge density at 
P is positive, then P is a point of stable equilibrium for 
particles of negative charge. 

This is the formulation for electrical charges in two di- 
mensions; for vortices, the sense of the rotation of the 
vortex does not matter of course, since the vortex inter- 
acts with its own image charge distribution. Hence if the 
curvature at P is negative, then a vortex will be trapped 
there. 

Similar reasoning can be used to show that a gener- 
alization of Eq. ([30]) . the "Gaussian Monkey Saddle" 

given by h{x, y) = -%5R(a; — iy)^e ^''o , traps vortices in 
an energy well of the form E = Ea + ^^^ht + {cnst. + 

cnst. cos &0)r^ -I- . . . . The reasoning needs to be modified 
because the curvature at the origin of the monkey-saddle 
is zero and the trapping is due to the negative curvature 
near the origin. 

At a point of low symmetry (such as the origin in Eq. 
(I30p when A 7^ 1), the character of an extremum depends 
on the charge distribution elsewhere, since the previous 




FIG. 11 The geometry-defect interaction energy of a vortex 
on the saddle surface with A = 17. One notices a slight insta- 
bility in the x direction. 



argument only implies that c -I- e > 0. Fig. [TUD suggests 
that a vortex at the origin is destabilized by its repulsion 
from the positive curvature above and below the origin, 
which is not balanced by enough positive curvature to the 
left and right. In fact, more detailed calculations show 
that the range of A for which the origin is an energy 
minimum is y/&b — 8 < A < -s/BS -I- 8; the origin is a 
saddle point outside this range, as is just barely visible 
for the case of A = 17 in Fig. [TT] (Likewise, for negative 
values of A, the origin is a maximum when VGS — 8 < 
—A < \/65 + 8, but a saddle point outside this range.) 

These results follow by changing the integration vari- 
ables to £^ = x — x' ,rj = y — y' in Eq. (j33p and then expand- 
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ing to second order about the origin {x,y) = (0,0). The 
integral expressions for second derivatives of the energy 
can be evaluated explicitly, 



Ex{x,y)^K'K[ 



2 1 + A2 - 6A 



16 



T^ A^ - 1 

4 Arl 4 Ar} 



r,„2i-^' 



Go)] (35) 



'0 







where Gq = — 4A-^ is the curvature at the origin. In 

Appendix IBJ we determine the geometric potential for 
arbitrary x and y in (unwieldy) closed form. 

D. Hysteresis of vortices and trapping strength 

The geometrical interaction has its maximum strength 
when the Gaussian curvature is the strongest. However, 
the geometric charge (i.e., integrated Gaussian curva- 
ture) of any particular feature on a surface has a strength 
roughly equivalent at most to the charge of one or two 
vortices. Eq. [21] therefore suggests that the force on a 
vortex due to a feature of the surface is less than the 
force due to a couple vortices at the same distance. Pre- 
cise limits on the strength of the geometric interaction 
will be stated and proven in Section IVIl for arbitrary 
geometries. 

As a consequence the geometric interaction has its 
most significant effects when the number of vortices is 
comparable to the number of bumps and saddles on the 
surface, so that the geometrical force is not obscured by 
interactions with the other vortices. This is a recurring 
(melancholy) theme of our calculations, to be illustrated 
in Section IIII.DI for arrangements of vortices in a rotat- 
ing film. The current section illustrates the point by 
discussing hysteresis on a surface with multiple saddle 
points (i.e., traps). If a vortex- free supcrfluid film is 
heated, many vortices form in pairs of opposite signs. 
When it is cooled again, positive and negative vortices 
can remain trapped in metastable states in the saddles, 
but even with the strongest curvature possible, the argu- 
ment above suggests that not more than one vortex can 
be trapped per saddle. 

The effectiveness of the defect trapping by geometry 
is determined mainly by the ratio of the saddle density 
to vortex density. As shown in the previous section, the 
geometric energy near the center of a vortex trap with 
90° symmetry is given by 



E{r) « -K\G,\r\ 



(36) 



The force on the vortex found by differentiating the en- 
ergy reads 



F{r) « --K\G^\r. 



(37) 



Eq. (I37p shows that the trap pulls the vortex more and 
more strongly as the vortex is pulled away from the cen- 
ter, like a spring, until the vortex reaches the end of the 



trap at a distance of the order of ro where the force starts 

2 

decreasing. Since Gq ~ % (which is valid for a small as- 

pect ratio a), "the spring breaks down" when the vortex 
is pulled with a force greater than 



F„ 



Firo) 



Kg' 
ro 



(38) 



Let us consider a pair of saddles separated by distance 
d. It is possible that one vortex can be trapped in each 
saddle even for a small a provided that d is large enough. 
Remote vortices do not interact strongly enough to push 
one another out of their traps. The Coulomb attrac- 
tion or repulsion of the vortices must be weaker than the 
breakdown force of the trap Fmax, i.e., -^^ > ^. The 
minimum distance between the two saddles is therefore 






ro_ 

2 ■ 



(39) 



Let us find the maximum density of trapped vortices 
that can remain when the helium film is cooled through 
the Kosterlitz-Thouless temperature. Let us suppose 
there is a lattice of saddles forming a bumpy texture like a 
chicken skin. Suppose bumps cover the whole surface, so 
that the spacing between the saddles is of order tq. Then 
not every saddle can trap a vortex; the largest density of 
saddles which trap vortices is of the order of l/(i^„j„, so 
the fraction of saddles which ultimately contain vortices 

2 

is at most 1° ex a^ . Note that not as many vortices can 
be trapped if they all have the same sign, since the in- 
teractions from distant vortices add up producing a very 
large net force. On the other hand, producing vortices of 
both signs by heating and then cooling the helium film re- 
sults in screened vortex interactions which are weaker and 
hence less likely to push the defects out of the metastable 
states in which they are trapped. 



III. ROTATING SUPERFLUID FILMS ON A 
CORRUGATED SUBSTRATE 

A. The effect of rotation 

Suppose that the vessel containing the supcrfluid layer 
is rotated around the axis of symmetry of the Gaussian 
bump with angular velocity 17 = f2 z, as might occur 
at the bottom of a spinning wine bottle. The container 
can rotate independently of the supcrfluid in it because 
there is no friction between the two. However, a state 
with vanishing supcrfluid angular momentum is not the 
ground state. To see this, note that the energy, Erot^ in 
a frame rotating at angular velocity Jl is given by: 



Er 



E L n 



(40) 



where E is the energy in the laboratory frame and L 
is the angular momentum. Hence Erot is lowered when 
L • J7 > 0, that is, when the circulation in the super- 
fluid is non- vanishing. This is achieved by introducing 
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quantized vortices in the system (see Eq.®), whose mi- 
croscopic core radius (of the order of a few A) is made 
of normal rather than superfiuid component. The energy 
of rotation, L • il, corresponding to a vortex at position 
X, y on the bump can be evaluated from 



Lz = Ps / dxdy^/g{x~y){xvy - yv^) . (41) 

Js 

Upon casting the integral in Eq. (|4ip in polar coordinates 
r, (f) and using the identity 



{xvy - yvx) = rcj) 



we obtain 



L. 



R 



dry/g{r) (p du°'Va 
Jc 



(42) 



(43) 



where R is the size of the system. The line integral in 
Ea. (|^5)) of radius is evaluated over circular contours of 
radius r centered at the origin of the bump. The circu- 
lation vanishes if the vortex of strength n at distance r„ 
is not enclosed by the contour of radius r: 



du°'va = nK6{r — r„) 



(44) 



Upon substituting in Eq. (|43p . we obtain 



R 



Lz = npsK I dr^/g{r) 
npsK 



27r 



{A{R)-A{r,)) 



(45) 



where A{R) is the total area spanned by the bump and 
A{rv) is the area of the cup of the bump bounded by the 
position of the vortex. Thus, after suppressing a con- 
stant, the rotation generates an approximately parabolic 
potential energy Eii{r) (see Fig. p^ ) that confines a vor- 
tex of positive index n close to the axis of rotation as in 
flat space: 



En{rvj = n A(r„) , 



TO4 



(46) 



where a consta nt has been n eglected. One recovers the 
flat space result (jVinen!, Il969!) by setting a equal to zero. 
Ea. (|45p has an appealing intuitive interpretation as the 
total number of superfiuid atoms beyond the vortex, 
— {A{R) — A{r)), times a quantum of angular momen- 
tum h carried by each of them. The closer the vortex is 
to the axis, the more atoms there are rotating with the 
container. 

Above a critical frequency ili, the restoring force due 
to the rotation (the gradient of Eq. (|1S|) ) is greater than 
the attraction to the boundary. The energy of attraction 
to the boundary is approximately 7riirin(l — j^), where 
we assume the aspect ratio of the bump is small so that 
the flat space result is recovered. Upon expanding this 



boundary potential harmonically about the origin and 
comparing to Eq. (|^5)) . one sees that 



ni 



li?2- 



(47) 



Above r^i, the origin is a local minimum in the energy 
function for a single vortex, though higher frequencies are 
necessary to produce the vortex in the first place. What 
determines the critical frequency for producing a vortex 
is unclear. There is a higher frequency fl'i ~ ';S^^''^^^ 
at which the single vortex actually has a lower energy 
(according to Eq. (|15)l ) than no vortex at all, but critical 
speeds are ra rely in agreement with the measured values 
( Vinenl . Il963l ). In the context of thin layers, it is likely 



that a third, much larger critical speed ficrit ^ mRD ' ^^ 
necessary before vortices form spontaneously, where Dq 
is the thickness of the film (see Sec. llV.Bp . 



B. Single defect ground state 

The equilibrium position of an isolated vortex far from 
the boundary is determined from the competition be- 
tween the confining potential caused by the rotation and 
the geometric interaction that pushes the vortex away 
from the top of the bump. The energy of the vortex, 
E{r), as a function of its radial distance from the center 
of the bump is given up to a constant by the sum of the 
geometric potential and the potential due to rotation. 



E{r) 
K 



-nUcir) 



Mr) 
A2 



(48) 



where we have ignored the effects of the distant boundary, 
boundary effects are discussed in the next section. The 
"rotational length" A is defined as 
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A helium atom at radius A from the origin rotating with 
the frequency of the substrate has a single quantum of 
angular momentum. The geometric contribution to E(r) 
( see Fig. [T2l) varies strongly as the shape of the sub- 
strate is changed. The rotation contribution to E(r) 
confinement (see Fig. [T51) varies predominantly as the 
frequency is changed; near the center of rotation, where 
the substrate is parallel to the horizontal plane, the ro- 
tational contribution barely changes as a is increased. 

As one varies a (fixing vq and Q) there is a transition to 
an asymmetric minimum. In fact. Fig. 1 141 reveals that for 
a greater than a critical value ac the total energy E{r) 
assumes a Mexican hat shape whose minimum is offset 
from the top of the bump. The position of this minimum 
is found by taking a derivative of Eq. P5|) with respect 
to r: 
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FIG. 12 Plot of minus the geometric potential —Ucir) for 
a — 0.5,1,1.5,2. The arrow indicates increasing a. The 
radial coordinate r is measured in units of A and ro — X. 




FIG. 13 Plot of the area of a cup of radius r for a = 
0.5,1,1.5,2. The arrow indicates increasing a. The radial 
coordinate r is measured in units of A and ro — A. 
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FIG. 14 Plot of E{r) measured in units of _ft' = 
varied. In these units, the thermal energy ksT is less than 0.1 
below the Kosterlitz-Thouless temperature, for 200A films. 
The radial coordinate r is measured in units of A and ro — ^■ 
Note that this plot is a 2D slice of a 3D potential. For a < Oc, 
E{r) is approximately a paraboloid while, for a > Uc, we have 
a Mexican hat potential. 




Now ^ can be shown to equal 2TTr\/l + h^ by differenti- 
ating Eq. 135] and ^^, which is the same as FgVI + h'^ 
can be evaluated by substituting for Fq from Eq. [25l 
This leads to an implicit equation for the position of the 
minimum, r„j, namely 



— = sm( ) 

A ^ 2 ' 



(51) 



Here d[r\ defined in Sec. III.Bi is the angle that the tan- 
gent at r to the bump forms with a horizontal plane. A 
simple construction allows one to solve Eq. ([5T|) graphi- 
cally by finding the intercept (s) of the curve on the right- 
hand side with the straight line of slope j on the left- 
hand side (see Fig. [T6|) . A brief calculation based on 
this construction shows that for a > ac = -r'^, there 
are two intercepts: one at ?' — (the maximum) and 
one at r = r™, the minimum; whereas for a < ac only 
a minimum at r = exists exactly like in flat space. 
It is possible to go through this second order transition 
by changing other parameters such as the rotational fre- 
quency. See Figs. [14] and [15] for illustrations of how 
the transition occurs when the shape of the substrate 
is varied. More details on the choice of substrate pa- 



t2 

FIG. 15 Plot of Eir) in units of — ^ versus r as ro is varied. 
The aspect ratio is kept fixed a.t a — 2 while the range of 
the geometric potential (corresponding to the width of the 
bump) is varied so that ro = 0.2, 0.4, 0.6, 0.8, 1 in units of 
A. As ro decreases, the geometric force becomes stronger, so 
the system goes through a transition analogous to the one 
displayed in Fig. 1141 



rameters are given in Sec. IIVI Once these parameters 
have been chosen, changing the frequency would likely 
be more convenient; Fig. \W\ shows how the equilib- 
rium position of the vortex varies. If the vortex position 
Tm can be measured precisely as a function of Q and if 
there is not too much pinning, then the geometrical po- 
tential can even be reconstructed by integrating Ug — 

-It 2^r,^in'Wi + h'{r^{n'W^ 

which follows from Eq. \UU\ 



■{n')dn' + cnst. 



C. Multiple defect configurations 

As the angular speed is raised, a cascade of transi- 
tions characterized by an increasing number of vortices 
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FIG. 16 Graphical method for determining equihbrium po- 
sitions of one vortex. The equilibrium position is at the in- 
tersection of sin -^ and j. If we fix ro and set a — 1, the 
rotational frequency will control the position of the vortex. 
The four lines correspond to Q - 
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occurs just as in flat space. In order to facilitate the 
mathematical analysis we introduce a conformal se t of 
coordinates {7^(r),(/)} (see (IVitelli and Nelsonl . [20041) for 
details). The function Tl{r) corresponds to a nonlinear 
stretch of the radial coordinate that "flattens" the bump, 
leaving the points at the origin and inflnity unchanged: 
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Note the unwonted appearance of the geometric potential 
Uc{r) playing the role of the conformal scale factor; this 
surprise is the starting point for our derivation of the 
geometric interaction in Section IV. Al The free energy 
of Ny vortices on a bump bounded by a circular wall at 
distance R from its center is given by 
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The Green's function expressed in scaled coordinates 
reads 
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where (j)i is the usul polar angle and the dimcnsionless 
vortex position ti is defined by 



' " TZiR) ■ 



(55) 



Eq. ([55)1 is now cast in a form equivalent to the flat space 
expression apart from the third term which results from 



the curvature of the underlying substrate and vanishes 
when a = 0. However, we emphasize that the Green's 
function T^ also is modified by the curvature of the sur- 
face and thus depends on a. 

The contributions from the second term and the nu- 
merator of the Green's functions in the first term account 
for the interaction of each vortex with its own image and 
with the images of the other vo rtices present on the bump 
(see (IVitelh and Nelsonl I200I ). If i? > ro, and aU the 
vortices are near the top of the bump (i.e., Vi ^^ vq) then 
these boundary effects may all be omitted when deter- 
mining equilibrium positions, as the forces which they 
imply are on the order of K^, small compared to the 
intervortex forces and geometric forces, which have a typ- 
ical value of — . 

'"0. . 

Let us imagine rotating the supcrfiuid, so that each 
vortex is confined by a potential of the form Eq. (|46|l . 
In fiat space, the locally stable configurations usually 
invol ve concentric rings of vortices (Campbell and Zifil , 
Il979[ ). In particular, there are two stable configurations 
of six vortices. The lower energy configuration has one 
vortex in the center and five in a pentagon surrounding it. 
The other configuration, six vor tices in a hexagon, has a 
slight ly higher energy, and Ref. (JYarmchuk and Packardl 
Il982l ) saw the configuration fiuctuating randomly be- 
tween the two, probably due to mechanical vibrations 
since thermal oscillations would not be strong enough to 
move the vortices. (The experiment used a Dq — 2 cm 
high column of superfiuid; if one regards the problem 
as two dimensional by considering fiows that are homo- 
geneous in the z direction, ps — -D0P3 is so large that 
K = -^ps is on the order of millions of degrees Kelvin.) 
There are no other stable configurations. However, on the 
curved surface of a bump, there are several more config- 
urations which can be found by numerically minimizing 
Eq. ([55)1 ; the progression of patterns as a increases de- 
pends on how tightly confined the vortices are compared 
to the size of the bump, as illustrated in Fig. [171 If the 
vortices are tightly confined, the interactions of the vor- 
tices (which are different in curved space) stabilize the 
new vortex arrangements. If the vortices are spaced far 
apart, the geometric interaction between the bump and 
the central vortex causes a transition akin to the decen- 
tering transition in the previous section. 

For example, if fi = 9-^, then at a = 0, the five 
off-center vortices start out in a ring of radius .6ro. This 
pentagonal arrangement (see Fig. [TTK) is locally stable 
for a < ai — 2.7. However, for a > a2 = 2.1, another 
arrangement with less symmetry is also stable (see frame 
B of Fig. [T7| . and above ai it takes over from the pen- 
tagon. For a2 < a < Qfi, both arrangements are locally 
stable, with the asymmetric shape becoming energeti- 
cally favored at some intermediate aspect ratio. (There 
is also a third arrangement which coexists with the less 
symmetric arrangement for the larger aspect ratios, seen 
in the frame C of Fig. [T71) In the plane, the configu- 
ration labelled B, for example, is unstable, because the 
outer rectangle of vortices can rotate through angle e. 
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decreasing its interaction energy with the two interior 
vortices while keeping the rotational confinement energy 
constant. (That the interaction energy decreases can be 
demonstrated by expanding it in powers of e.) Because 
the Green's functions are different on the curved surface 
(they do not depend solely on the distance between the 
vortices in the projected view shown), figures B and C 
are stabilized. 

At lower rotational frequencies, the equilibria which 
occur are even less symmetric. For Vi = -^, the vor- 

tices form a pentagon of radius tq when the surface is 
flat. This pentagon is far enough away that it has a mi- 
nor influence on the central vortex, which undergoes a 
transition similar to the one discussed in Sec. IIII.BI At 
a'l = 1.4, the central vortex moves off-axis (the transi- 
tion is continuous), causing only a slight deformation of 
the pentagon (see frame D of Fig. \T7\ . As for the sin- 
gle vortex on a rotating bump, the geometric potential 
has pushed the central vortex away from the maximum, 
and the other vortices are far enough away that they are 
not influenced much. At higher aspect ratios, the fig- 
ure distorts further, taking a shape similar to the one 
which occurs for il = 9-^, but offset due to the geo- 
metric interaction. For these two rotational frequencies 
the hexagonal configuration is less stable than the pen- 
tagon; it will not take the place of the pentagon once 
the pentagon is destabilized. The hexagon is of course 
metastable for nearly flat surfaces. 
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FIG. 17 Arrangements of 6 vortices that can occur on a 
curved surface. A circle of radius ro is drawn to give a sense 
that the confinenient is tighter in the top row (fi — 9-^) 

mrg 

than in the bottom row {fl = -^). The upper row shows the 

patterns which occur at large angular frequencies (17 = 9-2-^). 

The transition from the pentagon to the rectangle with two in- 
terior points is discontinuous, and there is a range of aspect ra- 
tios 2.1 < a < 2.7 where both configurations are metastable. 
The third configuration is nearly degenerate with the second 
configuration. The lower row shows the configurations which 
occur for fi — — r as a increases. The first transition is con- 

mrg 

tinuous and caused by the central vortex's being repelled from 
the top by the geometric interaction. The third configuration 
is similar to the second large Q configuration but the effect of 
the geometric repulsion is seen in its asymmetry. 



D. Abrikosov lattice on a curved surface 



As in Section IILDI the geometric potential will have 
significant consequences only when the number of vor- 
tices near each geometrical feature such as a bump is 
of order unity. As an example, consider the triangular 
vortex lattice that forms at higher rotational frequencies 
(ft ^ — ^ is the criterion for a large number of vortices 

mrg 

to reside on top of the Gaussian bump). In flat space, 
the vorte x number density is approximately constant and 
equal to (JTillev and Tillevl . Il990f ) 
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At equilibrium, the force exerted on an arbitrary vor- 
tex as a result of the rotation exactly balances the force 
resulting from the interaction with the other vortices in 
the lattice and from the anomalous coupling to the Gaus- 
sian curvature. We can determine the distribution on a 
curved substrate by making the continuum approxima- 
tion to Eq. P^ . The sum of delta- functions a gets re- 
placed by 2711/ (r) and the self-charge subtraction can be 
neglected in the continuum approximation. The Gaus- 
sian curvature can be neglected because it is small com- 
pared to the large density of vortex charge. Upon apply- 
ing Gauss's theorem to the vortex charge distribution in 
an analogous way to Section [ILAI we find that the force 



on a vortex at radius r is given by 

1 f^ 

F^ = - 4:7:'^Ki'{r'yy/lTl^dr' (57) 

^ Jo 

while the rotational confinement force, obtained by dif- 
ferentiating Eq. (|45|) . is 



2nhn 
Fn = -ps r. 



(58) 



Balancing the two forces leads to an areal density of vor- 
tices. 
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Eq. ([59)) has a succinct geometric interpretation: the 
vortex density i'{r) arises from distributing the vortices 
on the bump so that the projection of this density on the 
xy plane is uniform and equal to the fiat space result. The 
superfiuid tries to mimic a rigidly rotating curved body as 
much as possible given that the fiow must be irrotational 
outside of vort ex cores as for the case o f a rotating cylin- 
der of helium (JTillev and Tillevl . Il990l ). To check this, 
first notice that the approximate rigid rotation entails a 
flow speed of f2r at points whose projected distance from 
the rotation axis is r. Hence, the circulation increases 
according to the quadratic law f v • dl — 2TTilr'^ . Since 
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this quantity is proportional to the projected area of the 
surface out to radius r, the discretized version of such a 
distribution would consist of vortices, each with circula- 
tion K = ^^, with a constant projected density — as 
in flat space. This result can be generalized with some 
effort to any surface rotated at a constant rate, whether 
the surface is symmetric or not. 

The geometric force has to compete with the interac- 
tions among the many vortices expected at high angu- 
lar frequencies. More precisely, the maximum force at 
radius tq according to Eq. (|25p is of order -^ while 
the force due to all the vortices Eq. (|57p is of order 
^ (2^)(7rrJj^(2^,.(o)) ^ 2tt^ KroiyiO) . The last expression 

greatly exceeds -^ in the limit of high angular velocity. 
The geometrical repulsion leads to a small depletion of 
the vortex density of the order of one vortex in an area 
of order irr^ . 

The vortex arrangements produced by rotation are 
reminiscent o f Abrikosov lattices in a superconductor 
( Vinenl . Il969f ) . In fact an analogy exists between a thin 
film of superconductor in a magnetic field and a rotat- 
ing film of superfluid. A major difference between bulk 
superfluids and bulk superconductors is that the vortices 
in a bulk superconductor have an exponentially decay- 
ing interaction rather than a logarithmic one because 
of the magnetic field (produced by the vortex current) 
which screens the supercurrent. The analogy is more ap- 
propriate in a thin superconducting film, where the su- 
percurrents (being confined to the film) produce a much 
weaker magnetic field. In fact, Abrikosov vortices in a 
superconducting film exhibit helium-like unscreened log- 



arithmic interactions out to length scales of order ^ = jj 
where A is the bulk Lon don pe netration depth and D is 
the film thickne ss (see ( Pearl , 1 1964) and, for a review, 
section 6.2.5 of ( Nelsonl . 12002 )1. Our results on helium 
superfluids without rotation therefore apply also to vor- 
tices in a curved superconducting layer in the absence of 
an external magnetic field. Curved superconducting lay- 
ers in external magnetic fields can be understood as well 
by replacing the magnetic field by rotation of the super- 
fiuid. Let us review the analogy between a container of 
superfluid helium rotating at angula r speed fl an d a su- 
perconductor in a magnetic field H (jVinenl . Il969l ) . Note 
that in Eq. gOl S is given by ips // d'^u^ {VOf and 
h\J9 is the momentum in the rest frame, p, although we 
are working in the rotating frame (the frame in which a 
vortex lattice would be at rest). For helium, the momen- 
tum in the rest frame p is related to the momentum in 
the rotating frame p' by the "gauge" transformation 



mr X $7 



(60) 



Similarly, in the case of a superconductor the momentum 
p in the absence of a magnetic field is related to the 
moment um p' in the pres ence of the field by the familiar 
relation (|Tinkhaml . 119961 ) 



where A is the vector potential. Comparison of Eq. ([60]) 
and Eq. (pT|) suggests a formal analogy between the two 
problems. 
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Eq. ([61]) establishes a correspondence between the an- 
gular velocity Q, and the magnetic field H that allows to 
convert most of the relations we derived for helium to the 
problem of a superconducting layer, with the identifica- 
tion 
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Of course, one should keep the external magnetic field 
small so that a dense Abrikosov lattice does not form, 
since (as for superfiuids) when there are too many vor- 
tices, the curvature interaction is overcome by the vortex 
interactions. 



IV. EXPERIMENTAL CONSIDERATIONS 

Vortices in bulk fiuids are extended objects such as 
curves connecting opposite boundaries, rings or knots. A 
vortex interacts with itself and with its image generated 
by the boundary of the fluid. However, if the vortex is 
curved, such forces (the three-dimensional generalization 
of the geometric force) are usually dominated by a force 
which depends on the curvature of the vortex called the 
"local i nduction force." This force has a strength per unit 
length (|Saffmanl . fl99l of 



J LI A = TT^paKln — , 
rw^ Ka 



(64) 



where p^ is the bulk superfluid density and k is the cur- 
vature of the vortex at the point where this force acts. 
This force is in danger of dominating the long range forces 
because of the core size appearing in the logarithm. 

"Two-dimensional" regions are a special case of three- 
dimensional regions in which two of the boundaries are 
parallel and at a distance Dq much less than the radius of 
curvature of the boundaries. The two dimensional super- 
fluid density is given by ps — P3D0, and the interactions 
of the vortices should be captured by the two dimensional 
theory described in this paper once this substitution is 
made. A discrepancy will occur, however, if the bound- 
aries of the film are not exactly parallel because the vor- 
tices are forced to curve in order to meet both boundaries 
at right angles. In this case, there is a force which is a 
relic of the local induction force (see Sec. lIV.Cp . 
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where rg is the relevant curvature scale. According to this 
formula, vortices are attracted to the thinnest portions of 
the film. We will need to ensure that the thickness of the 
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film is uniform enough so that this force does not dom- 
inate over the geometric interactions we are interested 
in. 

There is a maximum film thickness for which the ge- 
ometric force is relevant. The most stringent require- 
ment arises from demanding that the van der Waals 
force causes wetting of the surface with a sufficiently uni- 
form film. Van der Waals forces compete against grav- 
ity, which thickens the superfluid at lower portions of 
the substrate, and surface tension, which thickens the 
superfluid where the mean curvature of the substrate is 
negative. Both gravity and surface tension thin the film 
on hills and thicken it in valleys, but if the film is thin 
enough, the van der Waals force can keep the nonuni- 
formity very small. Section IIV.BI discusses the critical 
speeds for the nucleation of vortices in thin films, which 
are typically hig her than those required i n long thin ro- 
tating cylinders ( Yarmchuk and Packardlll982^ . We as- 
sume that vorticity is not created from scratch, but from 
pinn ed vortices presen t even before the rotation has be- 
gun (iTillev and Tillevl . \imh . Finally in Sees. [IVXTl and 
lIV.DI a comparison is made between forces on vortex lines 
in three-dimensional geometries and on point vortices in 
two dimensions. 
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FIG. 18 Definition plot for a laminating film. h{x) is the 
height of the substrate above the horizontal surface at a point 
X = (a;, y), and hL(^) is the height of the upper surface of the 
film. -D(x) is the thickness of the film which (if the film has a 
slowly varying thickness) is given by (/i_l(x) — /i(x)) cos^(x) 
where 6'(x) is the local inclination angle of the substrate. 



be determined by minimizing the free energy F, 
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A. The Van der Waals force and thickness variation 

We start by providing an estimate of the variation in 
the relative thickness 



Dt-Do 
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for a liquid layer which wets a bump and apply it to 
thin helium films. Dt denotes the thickness on top of 
the bump and Do is the thickness far away. The wetting 
properties of very thin films {^ 100 A) of dodecane on 
polymeric fibers of approximatel y cylindrical shape have 
been thoroughly investigated in (JQuere et a/.l . 119891 ). We 
start by reviewing a theoretical treatment of the statics of 
wetting on rough surfaces by (Andelman et aLl . [l9883 . A 
film on a solid substrate that is curved has a mean curva- 
ture determined by the shape of the substrate, unlike in 
the case of a large drop of water on a non- wetting surface. 
By choosing an appropriate shape, the drop can adjust 
its mean curvature (and thereby balance surface tension 
against gravity). The shape is therefore described by a 
differential equation. A thin film on a solid substrate, in 
contrast, has approximately the same curvature as the 
substrate that it outlines. 

Consider a film that completely wets a solid surface. 
The surface itself is described by its height function /i(x), 
where x denotes a pair of Cartesian coordinates in the 
horizontal plane below the surface (see Fig. [T5|) . The 
height function for the liquid- vapor interface /il(x) can 



where 7, pa, and ^ are respectively the liquid- vapor sur- 
face tension, the total mass density, and the chemical 
potential (per unit volume). (Note that V here is not 
the covariant gradient for the surface; it is the gradient 
in the xy plane.) The second term describes the gravita- 
tional potential energy integrated through the thickness 
of the film. The second and fourth terms model the force 
between the helium atoms and the substrate assuming for 
simplicity a non-retarded van der Waals interaction. The 
last term involves an integral over interactions between 
pairs of points, one above the helium film and one in the 
substrate, but with no points in the liquid helium itself. 
This is equivalent to including interactions between all 
pairs of atoms contained in all combinations of the va- 
por, liquid and solid regions, as long as w(r) ~ — ar" ^ 
where a is th e appropriate com bination of parameters for 
these phases ( Andelman et 0/.I , [l988). 

Minimization of Eq. ([S7)) leads to a differential equa- 
tion for h]^{x.) that is a suitable starting point for evalu- 
ating the profile of the l iquid- vapor interface numerically 
( Andelman et a/.l . ll9 88). In what follows, we will instead 
work within an approximation valid when Dq <C rQ,ho] 
in this case, the curvature of the film is fixed. The 
local film thickness is described by -D(x) = (/il(x) — 
/i(x))/^l + |V/i(x)2)|, see Fig.dUl). We need to deter- 
mine how each contribution to the free energy per unit 
area at a point u is changed by an increase in thickness 
SD{x). 

First let us consider the variation of the van der Waals 
energy in order to understand how this attraction sets the 
thickness of the film. When the film thickens by 6D over 
a small area A of the film (centered at a;, 2:), the change in 
the van der Waals energy is given by —ASDIl{x) where 
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the disjoining pressure is 



n(x) = 
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dz'w(V(x-x')2 + (z-z')2). 



(68) 
For a film on a horizontal surface at /i = 0, the surface 
area and gravitational potential energy do not increase 
when the film is thickened. The equilibrium thickness is 
determined by balancing the variation of the chemical 
potential contribution, —fiA6D, against the disjoining 
pressure, giving fi — —Il{Do). The disjoining pressure 
obtained by integrating Eq. ([55]) for a flat surface is 
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Ah — TT^a is the Hamaker constant for the solid and 
the vapor intera c ting a cross a liquid layer of thickness 



Do (llsraelachvilil 119851 ). One sees that a negative value 
of Ah ~ TT^Q: is necessary for wetting. The equilibrium 
thickness is 



Dn = 




(70) 



(For example, liquid "^He on a CaF2 surface has Ah ~ 
— 10~^^ J, and has a liquid-vapor surface tension of 3 x 
10~^ J/m^.) When there is a bump on the surface, Eq. 
(|70p gives the equilibrium thickness far from the bump. 
Note that both Ah and fi are negative in this expression. 
Increasing /i therefore increases the thickness of the film 
as expected. 

Now let us continue by considering the effects of grav- 
ity and surface tension for a curved substrate. The in- 
crease in gravitational potential energy is psghSD, just 
because there is an additional mass per unit area of the 
fluid psSD at height h. (The additional elevation from 
adding the fluid at the top of the fluid that was already 
present can be ignored if the layer is very thin.) The 
variation of the surfa ce tensio n energy can be related to 
the mean curvature ( Kamienl . 120021 using the fact that 
the area of a small patch of the liquid vapor interface 
A{'x) (at a distance D from the substrate) is related to 
the corre sponding area of t he solid surface Ao{x) by the 
relation (JHvde et all Il997l) 

A(x) = Ao(x) [1 + 2iJ(x)D(x) + G(x)D2(x)] . (71) 

The second term is proportional to the mean curvature 
H — ^{ki + K2) oi the surface, where we use the con- 
vention that the principal curvatures ki, K2 are positive 
when the surface curves away from the outward-pointing 
normal. The last term, proportional to the Gaussian cur- 
vature, can be ignored relative to the previous term since 
it is smaller by a factor of — . The mean curvature of the 
upper surface of the fluid is nearly the same as for the 
substrate, so the energy required to increase the thick- 
ness of the film is 2^H{x.)SD. For example, at the top of 
the bump, an increased thickness leads to an increased 



area, so surface tension prefers a smaller thickness there. 
Gravity also thins the film at the top of the bump so that 
vortices are attracted to the top. 

Now we must balance these forces against the disjoin- 
ing pressure. The flat space form of the disjoining pres- 
sure is not significantly altered by the curvature of the 
substrate for very thin films. According to Eq. (p5|) . the 
disjoining pressure is the sum of all the van der Waals 
interaction energies between the points of the substrate 
and a fixed point at the surface of the film. The integral 
(evaluated in Appendix [C| for a point at a distance D 
away from the substrate shows 
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The curvature correction in the second term of Eq. (|72p 
arises (when 7J > as at the top of the bump) because 
the surface bends away from the vapor molecules which 
interact only with the very nearest atoms of the solid 
substrate. This effect is small if D^ « r^ and will be 
neglected in our estimates. 

We can now collect the various contributions to set up 
a pressure balance equation that allows us to estimate the 
relative change in layer thickness e defined in Eq. (j66p . 
This equation reads: 
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Apart from the lengths tq, /iq and Dq inherited from 
the geometry of the system, it is convenient to define 
three characteristic length scales 5, g and Ic, obtained by 
pairwise balancing of the first three terms of Eq. ((73|) : 
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(74) 



The last relation in E g. (17^ defines the familiar capillary 
120011) below which surface tension 



length (jGuvon et al\ 
dominates over gravity while the first and the second give 
the length scales involving the disjoining pressure. For 
"'He on CaF2, (5 ~ 10 A(for most liquids it is one order 
of magnitude less), g ~ 0.7 fj,m and Ic — OAmm. 

Upon substituting Eq. ([70)1 in Eq. ([75]) . we obtain an 
approximate relation between I?(x) and Dq, 
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This relation leads to an estimate of the relative change in 
layer thickness, e, valid for thin films (if we take a ~ 1), 
namely 



roS^ \ ^Ic' 



(76) 
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Now the thickness- variation force (Eq. ([65]) ') is small 
compared to the coupling to the geometry only if 



e < 



21n^' 



(77) 



where we estimated the maximum of the geometrical 
force to be ^'^" , see Eq. (^5]) . This hmit on e leads 
in turn to an upper bound for the film thickness Dq for 
each choice of tq. Assuming a ^ 1 and splitting into 
cases according to the size of the bump gives the limits: 
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(78) 


^ , ,^' ^l for ^0 > Ic 


(79) 
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For smaller bumps, surface tension plays the main role in 
creating thickness variation and for larger bumps, grav- 
ity has the largest effects. In order for the vortices to be 
easily observable, Dq should be as large as possible. Eqs. 
(175)1 and ([75]) show that the film can be made thickest 
(while retaining its approximate uniformity) when grav- 
ity and surface tension have comparable effects, rg ~ Ic- 
The optimal thickness (calculated with a = 1 and for 
the Gaussian bump; similar numbers are optimal for the 
saddle surface) is about 150 A at rg — .5 mm.^ A smaller 
value of To might be required if there is a lower limit Qmin 
on the rotational frequency as discussed in the next sec- 
tion, requiring a slightly thinner film, about 100 A. 

The restriction on the film thickness is the most serious 
obtacle to studying two-dimensional superflows experi- 
mentally. The method of observing vortices described in 
(JYarmchuk and Packard! . Il982^ requires the vortices to 
be long enough to be able to tr ap an observable number 
of ele ctrons. If the method of ( Yarmchuk and Packardl 
Il982l ) turns out to be unsuitable for thin films and an al- 
ternative method cannot be found, one might also study a 
saturated superfluid layer confined between two solid sur- 
faces. In this case, none of the considerations on wetting 
are relevant, and the "film" could have a large thickness. 
The two solid surfaces would have to be parallel to one 
another, and hence not congruent. (Congruent surfaces 
displaced by a fixed distance in the vertical direction lead 
to a thickness varying as cos 0(x) where 9 is again the in- 
clination angle.) The two surfaces would have to be very 
accurately shaped in order to make the film uniformly 
thick. 

Another concern is that a vortex may be pinned to 
an irregularity on the substrate strongly enough that it 



^ Interestingly, even without rotating the bump, there are two radii 
where the vortex could rest for a film thickness of 200 A. Then 
the confinement due to the varying film thickness and the ge- 
ometric repulsion are comparable in magnitude, producing an 
equilibrium off-center position for the vortex. 



will not move to the location favored by the geometrical 
force. On the other hand, the geometrical force is much 
stronger than random forces due to thermal energy. Even 
with a film as thin as hundreds of Angstroms, the geo- 
metric force is very strong. With ps — PsDq — .2g/ccDo, 
which assumes that the superfiuid density is not depleted 
too much by thermal effects or by the thinness of the 

fc2 

film, the value oi K = ^^^ is about 40 Kelvin. Since 
the potential wells which trap the vortices on the saddle 
surface or on the rotating Gaussian bump have depths 
on the order of Ka^ the geometric force will be strong 
enough to prevent the vortex from wandering out of the 
trap due to thermal Brownian motion except very close to 
the Kosterlitz-Thouless transition where ps is depleted. 



B. Parameters for the rotation experiment 

A practical consideration can limit the thick- 
ness of the film further. In the experiments of 
( Yarmchuk and Packardl . 1 1982f) . uniform rotations slower 
than flmin = 1 rad/s or so were hard to attain. The 
maximum value of fi for which the geometrical force can 
displace the vortex as in Sec. IIII.BI can be determined 
by taking the limit r,„ -^ in Eq. ([51]) : this shows that 
for a given tq, the vortex is displaced from the top of the 
bump only if 
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Hence, if details of the rotational apparatus require that 
fi > 1 rad/s, then the size of the bump should be less 
than .1 mm (for a = 1). Consequently Eq. (j78p requires 
even thinner films than for the case discussed in the pre- 
vious section where the bump size is determined purely 
by the natural forces of gravity and surface tension and 
is equal to the capillary length. 

An additional issue is whether it is possible to cre- 
ate just one vortex reliably. In practice, the num- 
ber of vortices in a rotating film is a property of 
a history-dependent metastable state whose dynamics 
are not completely und erstood (see Chapter 6 of Ref. 
()Tillev and Tillevl . Il990f )). We can discuss this in the 
context of a flat rotating disk. The critical frequency for 
local stability of a single vortex at the center of the disk 
is ilstabi = ^Ri ■ There is a barrier to creating such a 
vortex, so the frequency initially would have to be raised 
up to a higher frequency in order for a vortex to form 
at the boundary and then move in to the center of the 
rotating helium. The critical frequency for creating vor- 
ticity at the boundary is rJcrit- It is not clear what this 
critical speed is. Perhaps fJcrit for a thin disk of helium is 
determined by the height D of the disk; the critical linear 



speed"^ would then be R^cnt ^ ^j^ ( Vinenl . [l963( ). Now 



At the critical speed, vortices are believed to form by break- 
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because Jlcrit is so much greater than ^stabi, a lattice of 
many vortices will form. According to Eq. ([55]) . there 

will be about 



mnR-' 



R 



In — of them. By slowing down 



h D 

to just above ^stabi, where the rotational confinement 

is not strong enough to confine more than one vortex 

(even metastably), one would hope to retain just a single 

vortex. 



C. Films of varying thickness from the three-dimensional 
point of view 

In a chamber of arbitrary shape the length of a vortex 
line crossing through the chamber changes as the vortex 
moves, and therefore knowledge of the core energy per 
unit length is crucial for determining the forces experi- 
enced by the vortex. The energy of a length D isolated 
vortex in a cylinder of radius R including the core energy 
is given by 
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where a = ae fs""^ is the core radius rescaled to account 
for the effects of the core energy per unit length, ecs- 

A vortex line connecting two approximately parallel 
parts of the boundary feels a force as a consequence of 
the variation of Eune ■ Let the two nearly parallel bound- 
aries be S and S' as illustrated in Fig. ^^ Let us define 
D{x) to be the length of the line which is perpendicular 
to iS at X and which extends to x' on S' (see Fig. UM . If 
a vortex is attached to 5 at x then since the fluid flow is 
fastest at the vortex, the energy depends primarily on the 
thickness of the helium film where the vortex is located. 
This thickness is approximately given by D{x) no matter 
how the vortex connects the two surfaces (unless it is ex- 
tremely wiggily). A useful fact is that, if the two surfaces 
are at a constant separation, then the line between x and 



ing away from pinned vortices. This implies that the crit- 
ical velocity is determined by the shortest dimension H of 
the chamber which is perpendicular to the flow of the normal 
fluid. Th e critical velocity is estimated i n the remanent vortic- 
ity theorvl lGlaberson and DonnellvUfQGBT) by assuming that ring 
vortices break off of pinned vortices stretching across this short- 
est dimension and then expand until they hit the surfaces of the 
superfluid and break into two line vortices with opposite circu- 
lations. The critical speed is the speed at which a ring vortex 
of size H would expand due to the motion of the fluid rather 
than contract. The critical velocity is determined by balancing 
the Magnus force of the fluid moving past the vortex (which ex- 
pands the vortex) against the attraction of the vortex for itself, 
and it is Vn = —2— In ^r- ■ Thus the critical speed for a narrow 
cylinder is determined by setting H equal to the radius R of the 
cylinder so uJc = Vc/R = r~fv7 In — , as seen experimentally in 
l|Yarmchuk and Packardl . Il982 ). For a thin film on a disk, the 
shortest dimension is the height, so H = D. 




FIG. 19 Illustration of the definition of -D(x), for nonparallel 
surfaces. The distance between S' and 5 is measured along 
the segment which is perpendicular to S at x. The segment 
meets <S' (obliquely) at a point which we call x'. A segment 
displaced to the right is also drawn, illustrating the derivation 
ofEq. ([82)1. 



x' is perpendicular to both surfaces. This can be derived 
by differentiating, -D(x)^ = |x — x'p. Upon displacing x 
shghtly, 2D5D = 2(x - x') • (<5x - (5x') = 2(x' - x) • Sx' , 
since the line connecting x to x' was chosen to be per- 
pendicular to S. Expressing the inner product in terms 
of the angle P indicated in Fig. [191 this reads, 



SD= |(5x'|cos/3« \Sii\cosp. 



(82) 



In order to calculate the kinetic energy of a vortex in a 
curved film of varying thickness, we will derive a more de- 
tailed form of Eq. ([5T|) for a curved film. First, the core 
energy due to the local disruption of superfiuidity has 
the same form, ec3-D(x) where -D(x) is the thickness of 
the film at the location x of the vortex. The major com- 
ponent of the energy is given by J^ps JJJ dCrdrd(l)^ 
where it is assumed that the flow is parallel to the bound- 
aries and roughly independent of ^, the coordinate nor- 
mal to S and approximately normal to S' . We integrate 
over ^ and divide the remaining two-dimensional integra- 
tion into two parts: 
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r>Lth 



where we are using polar coordinates centered on the lo- 
cation of the vortex. Here, Lth is the distance over which 
D varies appreciably so that Lth ~ ^^jj- We regard the 
thickness as a constant in the second term because the 
integral starts far enough away from the vortex that it 
does not depend on the speciflc thickness of the film at 
the location of the vortex. The first term may be ap- 
proximated by replacing D{r, cj)) by -D(x) (the thickness 
of the film at the location of the vortex, where the energy 
is very big). Therefore, 



E « e,3i?(x) + ^p3D{x) In 
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energy of the distant flow. (84) 
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FIG. 20 A figure for reference in determining tiie curvature of 
a vortex. If the upper and lower surface S and S' are nearly 
but not quite parallel, the curvature of a vortex connecting 
them can be estimated simply. Extend the tangents at the 
end-points X, Y of the vortex until they intersect at the cen- 
ter O of curvature of the vortex (which happens because the 
vortex is perpendicular to S,S'). The radius of curvature, 
OX, can then be found by using Eq. (|82p and noticing that 
OX = XX' tan /3 and sin/? « 1. 



The force obtained by taking the gradient of the kinetic 
energy reads 
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where Fg represents the variation in energy of the long 
range portion of the flow, which includes the geometric 
forces. 

We now show that the thickness variation force can 
be interpreted as the Biot-Savart self-interraction of the 
curved vortex. The net Biot-Savart force points towards 
the center of curvature of the vortex, and so favors shrink- 
age of its length. We will use the "local induction approx- 
imation" to the Biot-Savart self-interaction, given in Eq. 
(IM)l . This expression results from integrating the interac- 
tion of a particular element of the vortex line with nearby 
elements; the peculiar dependence on the core radius a 
arises from the need to place a cut-off in the diverging 
interaction for nearby elements. 

To see this, note that if the opposing boundaries 5, S' 
are not at a constant distance, the vortex must curve in 
order to connect them, because it meets both boundaries 
at right angles. The resulting curvature of the vortex, k, 
is equal to ^^ (see Fig. [20l), so that the force per unit 



length given by Eq. ([M]) . multiplied by D, agrees with 
the force derived by the energy method, Eq. ([55]) . On the 
other hand, for films of constant thickness, the surfaces 
S, S' can be connected by straight vortices. According to 
the Biot-Savart law a straight vortex has zero interaction 
with itself. In addition, there is still an interaction energy 
with its image, which is where the geometric force, the 
third term in Eq. I85i comes from. When the helium has 
boundaries, a distribution of "image vorticity" beyond 
the boundaries can be su pplied to simul ate the effects of 
the boundary conditions (|Saffmanl . ll992( ). The integral of 
the Coulomb interaction with the distribution of surface 
curvature (Eq. (j29p is reminiscent of the long-range inte- 
gral of the Biot-Savart interaction with the distribution 
of vorticity beyond the surface. 
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FIG. 21 Sketch of the experimental geometry from Ref. 
(jVoll et alV '2006). This experiment may give some indica- 
tion about the difference in energy between a vortex stretch- 
ing from A to B and one stretching from C to D. A wire 
i s placed a. l ong A B and a vortex is formed around it. Ref. 
(jVoll et al\ . I2OO6I ) found that this vortex depins more easily 
than a vortex in a tube with a flat bottom. This suggests that 
the energy of a vortex (when there is no wire) would decrease 
by moving from AB to CD, but estimates seem to rule out 
this explanation. 



D. Vortex Depinning 

Because of the divergence of the force described by 
Eq. ([65|l in the limit a ^ 0, it will be hard to see ef- 
fects of the geometric f orce without film s of very nearly 
uniform thickness. In ( Vol! et all . l2006l ) an experiment 
is described in which a vortex extending along a wire in 
a helium-filled tube leaves the wire more easily when the 
wire is connected to a bump on the bottom of the tube 
than when the bottom is flat. In this case, any geometric 
repulsion from the bump should be negligible in compar- 
ison to the extra energy associated with the stretching of 
the vortex. (See Fig. [^TJ) Geometric energies should be 

of order Rps-^ where R is the radius of the tube and ap- 
proximately the height of the bump, but the vortex has 
to stretch an additional length on the order of R in order 
to leave the bump. Therefore, the extra kinetic energy 
from Eq. ([5T|) . Rp:^-^ In— , probably overwhelms the 
geometric effects. Thus, the energy barrier for detach- 
ment of the vortex line should increase when the wire is 
attached to a bump. 

This suggests that depinning is not simply caused by 
thermally activated barrier crossings (which would oc- 
cur at the rate proportional to e~^^/'^). In fact, for 
the length scale of the bump in the experiment (a few 
millimeters), the additional energy barrier due to the 
stretching of the vortex line is millions of degrees kelvin 
at the temperature where the vortex depins! The de- 
pinning may instead d epend on "remanent vorticity" 
( Tillev and Tillevl . ll990f ) in the form of extra pinned vor- 
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tices stretching from the top of the bump to the cyhnder's 
walls. This would decrease the energy barrier because the 
vortex could leave the wire by attaching itself to some of 
the vortices that already exist. 

In a more nearly two-dimensional geometry, Eq. [85] 
shows that observing the equilibrium positions of vor- 
tices in a film of varying thickness could shed light on 
the value of the core size and core-energy and whether 
they are fixed functions of temperature as most mod- 
els imply. The structure of vortex cores is still not 
well understood, and the re are several alternative models 
( Tillev and Tillevl . 119901 ). Experiments on the geometric 
force in contexts where the two contributions to Eq. ((85|) 
are comparable could give information about the effective 
core radius a. By allowing the film thickness to vary by 
about e — ^j- ~ j^^ , see Eq. ([77|) , one ensures that the 
vortex core size has a decisive effect on the equilibrium 
positions of the vortices (comparable to the effects of the 
long range forces). 



V. COMPLEX SURFACE MORPHOLOGIES 

Up to this point, our discussion has been confined 
to rotationally symmetric surfaces and slightly deformed 
surfaces for which the electrostatic analogy and pertur- 
bation theory can be successfully employed to determine 
the geometric potential. To investigate geometric effects 
that arise for strong deformations and for surfaces with 
the topology of a sphere, we adopt a more versatile geo- 
metric approach based on the method of conformal map- 
ping, often employed in the study of complicated bound- 
ary problems in electromagnetism and fluid mechanics. 
This approach also sheds light on the physical origin of 
the geometric potential. A concrete goal is to solve for 
the energetics (and the associated flows) of topological 
defects on a complicated substrate T, the target surface, 
whose metric tensor we denote by grab- This is accom- 
plished by means of a conformal map C that transforms 
the target surface into a reference surface R, with metric 
tensor gRab- The computational advantages result from 
choosing the conformal map so that i? is a simple surface 
(e.g. an infinite flat plane, a flat disk, or a regular sphere) 
that preserves the topology of the target surface. Figure 
[22] represents a complicated planar domain denoted by T 
which can be mapped conformally onto a simple annulus 
labeled by R. We will introduce all the basic concepts in 
the context of this simple planar problem before turning 
to the conformal mapping between target and reference 
surfaces which is represented schematically in F ig. 1231 
Such mappings can always be found in principle (JDavidl . 
I1989D . 

The conformal transformation will map the original 
positions of the defects on T, denoted by u, onto a 
new set of coordinates on R denoted hy U = C(u). 
In what follows, capital calligraphic fonts always indi- 
cate coordinates on the reference surface. For sufficiently 
small objects near a point u the map will act as a sim- 



ilarity transformation; that is, an infinitesimal length, 
dsT = \/ gTabdu°-du^, will be rescaled by a scale factor 
e"(") which is independent of the orientation of the length 
onT: 



dsR = e"(")dsT, 



(86) 



where dsR — yj gRAsdhl'^dlJ^ . This result in turn im- 
plies a simple relation between the metric tensors of the 
two surfaces: 



gRAB 



e ^ 'gTAB, 



(87) 



where we have assumed for simplicity that the coordi- 
nates used on the target surface are chosen so that cor- 
responding points on the two surfaces have the same co- 
ordinates U"^ . 

We will demonstrate that, once the geometric quantity 
uj{u) is calculated, the geometric potential of an isolated 
vortex interacting with the curvature is automatically de- 
termined. For multiple vortices, the energy consists of 
single-vortex terms and vortex-vortex interactions. On a 
deformed sphere or plane, the geometric potential reads 



Ei{ui) = -TmfKuj{ui), 



where K is the stiffness parameter defined in Eq. ([11 
For a deformed disk, there are boundary interactions not 
included in Eq. ([55]) . (We will not consider multiply 
connected surfaces here, but the single particle energy on 
a multiply connected surface has additional contributions 
which cannot be described by a local Poisson equation.) 
The interaction energy is 

V 

E2{u„Uj) = -2TTn,njK\n^^, (89) 

where Vij is the distance between the two image points 
on the reference surface. When the reference surface is 
an undeformed sphere (the other possibilities are a plane 
or disk), Vij is the distance bet ween the points aloni 
a cho rd rather than a great circlet Lubenskv and Pros 
Il992f ). We will show below that on a deformed plane uj 
is equal to Uq, but from now on we will use lu instead; 
the two functions are conceptually different, and are not 
even equal on a deformed sphere. 

Equation ([55]) is derived by the method of conformal 
mapping in section fV. Al and its computational efficiency 
is illustrated in section IV. Bl where the geometric poten- 
tial of a vortex is evaluated on an Enneper disk, a min- 
imal surface that naturally arises in the context of soap 
films, but whose geometry is distorted enough compared 
to flat space that it cannot be analyzed with perturba- 
tion theory. Changing the geometry of the substrate has 
interesting effects not only on the one-body geometric 
potential but also on the two-body interaction between 
vortices. In section IV. CI we use conformal methods to 
show how a periodic lattice of bumps can cause the vor- 
tex interaction to become anisotropic. In section FV.DI we 
demonstrate that the quantization of circulation leads to 
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an extremely long-range force on an elongated surface 
with the topology of a sphere. The interactfon energy is 
no longer logarithmic as in Eq. (|89p . but now grows lin- 
early with the distance between the two vortices. Indeed 
the charge neutrality constraint imposed by the compact 
topology of a sphere, blurs the distinction between geo- 
metric potential and vortex-interaction drawn in Equa- 
tions (|55|) and ([59)1 . This is most easily seen by bypassing 
vortex energetics on the reference surface (which is after 
all an auxiliary concept) and opting for a more direct re- 
statement of the problem in terms of Green's functions 
on the actual target surface coated by the helium layer. 
The interaction energy now reads 



Reference 



Target 



E2{ui,Uj) ^ 4TT^Kn^njr{ui,Uj), 



(90) 



and the single-particle energy takes the form of a self- 
energy 

E[{u,) = -nn^KUciu^) = nn^K /"/" G(u')r(ui, u^d^u', 

(91) 
where F is a Green's function for the surface that gen- 
eralizes the logarithmic potential familiar from two di- 
mensional electrostatics. For a deformed plane the two 
descriptions of the interaction energy are equivalent since 
the Green's function on a deformed plane can be obtained 
by conformal mapping, 



F(u„u,) = -;lln^. 
zvr a 



(92) 



We will see that the expressions for the single-particle 
energies are also equivalent. In contrast, for a deformed 
sphere, we show in section fV. El and appendix IDJ that the 
two formulations do not agree term by term { E[ ^ Ei 
and E'2 ^ E2), although the combined effect of one- 
particle and interaction terms is the same (up to an ad- 
ditive constant). Both self-energies and interaction ener- 
gies include effects of the geometry and explicit formu- 
las are provided on an azimuthally symmetric deformed 
sphere. Finally, in section FvH we present a discussion of 
some general upper bounds to which the strength of ge- 
ometric forces is subjected (even in the regime of strong 
deformations) which are useful in experimental estimates 
and which illustrate a major difference between electro- 
static and geometric forces: The former can always be 
increased by piling-up physical charges but the latter are 
generated by the Gaussian curvature that can grow only 
at the price of "warping" the underlying geometry of 
space. Too much warping either leads to self-intersection 
of the surface or a dilution of the long-range force. 



A. Using conformal mapping 

We start by proving a simple relationship between the 
total energies (including self- and interaction parts) Ex 
and Efi of two corresponding vortex-configurations on 




FIG. 22 The flow in a wiggily annular region (the target sub- 
strate T), obtained by conformally mapping the flow from a 
circular annulus (the reference substrate R). Since every pair 
of radial spokes in the first picture comprises the same energy, 
this is true in the conformal image as well. The small region 
in the constriction on the right manages this by compensating 
for its small area by having a high flow speed. 



the target and reference surfaces respectively: 



N 



ET = Eii-TiKY,u{vi,). 



(93) 



The right-hand side of Eq. ([M|) can be calculated for 
the reference surface and then subsequently decomposed 
into single- vortex and vortex- vortex interactions; several 
examples are worked out in detail in Sections IV. Bl and 
IV. CI The general approach is best illustrated by con- 
sidering the planar flow in the complicated annular con- 
tainer shown in Fig. [22l which can be tackled by confor- 
mally mapping it to a simpler circular annulus. 

The flow in the reference annulus is clearly circular, 
and it has the same - dependence as for a vortex. A 
crucial property of conformal transformation allows us 
to transplant this understanding of the reference flow to 
the target annulus. This property concerns the following 
mapping of the stream function x (see Eq. (llSp ) from 
the reference surface to the target surface: 



Xt(u) = Xb{U) 



(94) 



Once the coordinate change U = C(u) is found, the pre- 
scription to determine xt provided in Eq. ([M)) guarantees 
that the corresponding flow will be irrotational and in- 
compressible, as required. Visually, all we are doing is 
taking the streamlines on the reference surface, which 
are level curves of XR^ ^^id mapping them by C~^ to the 
target surface. We can informally state this first property 
of conformal maps as follows: 

1) The conformal image of a physical flow pattern is 
still a physical pattern. 

Note that any multiple of the mapped stream function, 
axR{C{\i)), corresponds to an irrotational and incom- 
pressible flow but in this case the rates of flow in the tar- 
get and reference substrate are different. Only the choice 
a = 1 ensures that both flows have the same number 
n of circulation quanta around the hole (or around each 
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vortex if some are present.) This follows from another 
basic property of conformal maps: 

2) In flow patterns related by a conformal map, accord- 
ing to Eq. \94^ , circulation integrals around correspond- 
ing curves are the same. 

This property follows from the fact that the contri- 
bution to the circulation from each element of the con- 
tour, V • dl, is conformally invariant. The infinitesimal 
length dl and the velocity v scale inversely to each other 
under conformal transformations and the angle between 
them is preserved by the map. To understand why, note 
that flow lines are compressed together (stretched apart) 
when they are mapped to the target space if the confor- 
mal parameter w is greater (less) than zero. As a result, 
the velocity increases (decreases) by the same factor as 
distances are decreased (increased). This heuristic argu- 
ment is confirmed by noting that the velocity is given by 
the covariant curl of the stream function, see Equations 
(fTS]) and ([94]) . By definition, the cova riant curl carries 
a mu ltiplicative factor of g^^^^ = e'^ (JDubrovin et all 
Il992l ). hence 



vt = e vr. 



(95) 



On the other hand, Idl^l = e~"|dlij| rescales in the op- 
posite way, as indicated in Eq. ([5S)) . 

The problem of finding the energy in the deformed an- 
nulus can now be reduced to a simple-rotationally sym- 
metric problem by appealing to a third property of the 
conformal mapping: 

3) The kinetic energies in corresponding regions are 
equal, provided the regions do not contain vortices. 

The proof of this statement relies on the previous dis- 
cussion: the kinetic energy contained in an element of 
the area of the surface dA is ^psv'^dA. By Eq. ((95| . 
v^ scales as e^" whereas dA scales by e~'^'^ making the 
energy conformally invariant. Figure [22] illustrates picto- 
rially that the energy density in the original flow on the 
target surface is smoothed out and simplified, its varia- 
tions being replaced by variations in the conformal scale 
factor. 

We now return to the energetics of flows containing 
point vortices. The starting point of our analysis follows 
from the defining property of a conformal map, namely 
that a conformal image of a small figure has the same 
shape as the original figure, while a larger shape becomes 
distorted (consider Greenland, which has an elongated 
shape, but appears to round out at the top in a Mercator 
projection, which is itself a conformal map). To quantify 
the size limits, note that if a shape has size /, uj changes 
by about IWlu across the shape. Thus, as long as 
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the mapping rescales the shape uniformly. The right- 
hand side is ordinarily of order L, the curvature scale of 
the surface. As a result we can conclude that 



4) The circular shape of the streamlines near a micro- 
scopic vortex core on a substrate of slowly varying cur- 
vature is preserved. 

On a deformed substrate with a flow induced by vor- 
tices, the flow speeds will increase or decrease not just 
depending on the distance to the vortex, but also de- 
pending on the shape of the surface. For example, the 
vortex on top of the bump in the example of Sec. III.AI 
has a flow that decays more slowly with distance than 
in flat space. Also, for a vortex well off to the side of a 
bump, if the bump's height h is larger than its width 2ro, 
it turns out that the flow pattern penetrates only up to 
an elevation of about ro up the side of the bump. 

The method of conformal mapping elucidates these ge- 
ometrical rearrangements of the flow pattern. To find 
the fiow pattern around the vortices at positions Uj, we 
find the flow pattern around vortices at the correspond- 
ing positions Ui on the reference surface and then map 
these streamlines onto the target surface by Eq. (|94p . 
The energies are not equal in this case, in spite of prop- 
erty 3. Property 3 does not apply to a region contain- 
ing vortex cores, because we would have to suppose the 
area of the cores on the reference surface to be greater 
by e^" and the energy in the cores to be smaller by a 
factor of e^'^, in order for the conformal relation Eq. 
^PsV^dAj- = hPsVJ^dAfi to continue to hold. In contrast, 
the core radius is fixed by the short-distance correlations 
of the helium atoms and the core energy is related to the 
interaction energy of the atoms. 

The vortex cores are not significantly affected by the 
curvature of the substrate; moreover, the whole flow pat- 
tern in the vicinity of the core is nearly independent of 
the location of a vortex. We observe that each vortex has 
a "dominion," a region where the flows are forced by the 
presence of the vortex to be 
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The leading term has the same form as one expects for a 
vortex in a rotationally symmetric situation, and the ef- 
fects of geometry are accounted for by Sv ; by dimensional 
analysis, this error is of the order of ^^ where L is the 
radius of curvature of the substrate (or possibly the dis- 
tance to another vortex or to the boundary, whichever is 
shortest). Therefore we can introduce any length I « L 
and note that I is then a distance below which the effects 
of curvature do not have a signiflcant effect (compared 
to the diverging velocity field). The geometry correction 
gives a contribution to the energy within this radius that 
is also small, as seen by integrating the kinetic energy 
over the annulus between a and I (using Eq. ([97]) ): 
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(98) 



where a is the core radius and e^ the core energy. The 
error term is quadratic in -^ because the integral over the 
cross-term from squaring Eq. (|97p cancels. 
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FIG. 23 Comparing the energy of T and R by splitting T up into portions I and O and using different maps to map tfiem to 
R. Left, a target surface, wliicli lias the topology of an infinite plane and is distorted by a three dimensional lump. Right, the 
reference plane (shown in the plane of the page). / consists of the interior of radius I disks in T, and is mapped rigidly to the 
heavily demarcated disks in R. O consists of the gray portion of T, and is mapped conformally to the gray portion of R. I 
and O contain the same kinetic energy as their images in R, but the images do not fit together perfectly. Thus the difference 
in energy between the flows in T and in R is determined by calculating how much energy is contained in the annuli which are 
either mismatched or covered twice. 



When we wish to find the kinetic energy of the super- 
flow, these near-vortex regions are thus the simplest to 
account for, as their energy is nearly independent of their 
positions relative features such as bumps on the substrate 
as long as 
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Since the target and reference configurations have the 
same number N of vortices, the energies contained within 
radius I of the vortices are the same: 



Ei, = E^, = T:NK\n- 
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In order to find the forces on a set of vortices, we need 
to account for all the energy of the vortices in regions 
away from the vortices where the flow has been affected 
by the curvature. Let us imagine cutting the target sur- 
face up into an inner region / (the union of the radius I 
disks around each vortex) and an outer region O (con- 
sisting of everything else) as illustrated in Fig. [^ We 
can map / to the reference surface by simply translating 
each of the disks so that they surround the vortices on R. 
The modifications to the flow are all in O; for example, 
streamlines there are deformed from their circular shape. 
These irregularities can be removed (or at least simpli- 
fied) by using the conformal mapping to map O to R, 
just as in the case of the annulus illustrated above. The 
inhomogcncity of the conformal map compensates for the 
irregularity of the flow. Some portions of O are expanded 
and some are contracted, but its circular boundaries are 
small enough (see Eqs. (1^^ , (IM)) ) that they are simply 
rescaled into circles of different radii (property 4) : 
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Now we have mapped / and O from the target surface 
to regions on the reference surface which contain the same 



energy. But these images of the regions on R do not flt 
together. The conformal map on O stretches or contracts 
each hole in it, to circles of radii e"'"'^/. These stretched 
edges (corresponding to the solid black circle in Fig. [^5]) 
do not fit together with the images of O, which have been 
moved rigidly from the target surface. The energies are 
related by 
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We must correct for the gaps and overlaps between the 
two image regions on R in order to relate the last expres- 
sion to En. If Lu{ui) > 0, there is an annular gap near 
vortex i; using Eq. ([57)1 (since this gap is part of the flow 
controlled by this vortex): 
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Summing all these contributions gives our desired result, 
Eq. dMI. 

We emphasize the energy difference is not produced 
within the cores, or anywhere near the vortices. In fact, 
the fact that the energies on T and R differ is a result 
of assuming that there is no change in the flow within a 
macroscopic distance / of the vortex. The scale I ^ L 
only has to be small compared to the geometry of the 
system and has no relation to the atomic structure of the 
core. On the other hand, taking I as small as possible 
has an elegant consequence: Eq. ([M)) actually has an 

error which is O (k (■^) j , smaller than O (k (-^j) j as 

predicted at flrst. Taking smaller values of I gives a more 
accurate result, since the conformal mapping (which does 
not suffer from the error in Eq. (|98)) ) is used to calculate 
the energy of a larger portion of the flow pattern. 

Now Eq. (|55|) shows that the position-dependent scale 
factor, a'(u), plays the role of a single-particle energy. 
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Additionally, this single particle energy can be regarded 
as the "geometric potential," since it turns out to be re- 
lated to the curvature in a way analogous to how the 
electrostatic potential is related to the charge. The func- 
tion to depends on the shape of the boundaries and on 
the curvature of the surface. A varying scale factor is 
necessary to map between surfaces with different distri- 
butions of curvature (such as planes with and without 
bumps). (A constant scale factor only rescales the curva- 
ture.) The curvature therefore depends on the variation 
oi uj. In fact it can be shown that 



Gt{u, v) = e2"("'")Gfl(J7(w, v), V{u, v)) 

+ , , ^, zdaJdet{gT,cd)gfdbUJ. 

y^det{gT,cd) ^ 

(105) 

The second term is the Laplacian of w as a function on 
the target surface^. The correspondence with electro- 
static potentials, with Gt{u, v) as a source, is clear if the 
reference surface is a plane or disk and Gu = 0. Then 
Eq. p05)) reduces to Eq. ^. 

For a deformed sphere or plane, single particle poten- 
tials come entirely from the second term in Eq. (|93p . The 
reference energy, corresponding to vortices on a sphere or 
plane cannot favor one position over another because of 
the homogeneity of these reference surfaces. The first 
term, Ej^, leads to the vortex- vortex interactions which 
depend only on the separation of the vortices on the ref- 
erence surface, again by symmetry; this energy is 



En = ^Ti K 
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where F depends on the reference surface: 
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where 7 is the angle between the two points. (2i?sin ^ is 
the chordal distance between the points, not the geodesic 
distance along the surface, as one might have guessed for 
the natural generalization of a Green's function to curved 
space.) The detailed derivation of the second formulation 
of the vortex interaction energies in terms of the Green's 
functions on the deformed surface (see Eqs. ([90|l and 
JT])) is contained in Appendix iDl 



* As a check of this identity, imagine reversing the roles of the 
reference and target surfaces. Then u) should be regarded as a 
function on the reference surface. This changes the Laplacian by 
a factor of e^" (because gji is replaced by Qt)- Also, the sign 
of LD should be reversed. Rearranging the equation now brings it 
back into the original form with T and R switched. 



B. Vortices on a "Soap Film" Surface 

There are experiinental and theoretical motivations for 
studying substrates shaped as minimal surfaces. An ex- 
ample of a minimal surface is easy to make by dipping a 
loop of wire in soap; the spanning soap film tries to min- 
imize its area. Vortices can be studied on a helium film 
coating a solid substrate whose surface has the shape of 
such a film. Such surfaces are characterized by a vanish- 
ing mean curvature, i?(x), so the contribution of the sur- 
face tension 2jH(ti) to the thickness variation equation, 
Eq. ([75)1 is drastically reduced. From the mathemati- 
cal point of view, there is a widely known parametriza- 
tion due to Weierstrass (|Hvde et all Il997[ ) which readily 
leads to an exact expression for the geometric potential 
of a vortex on such a surface. 

Weierstrass's formulae, which give a minimal surface 
for each choice of an analytic function i?(C), read: 
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The correspondence between this parametric surface and 
the complex variable (^ = X + iy is a conformal map, and 
the conformal factor can be expressed in terms of i?(C)- 
Therefore, the analysis of vortices on such a surface is 
not difficult at all when the X, 3^-plane is used as the 
reference surface. 

As an example, let R{C) be equal to L^ where L con- 
trols the size of the target surface. Then the surface 
produced is given in parametric form by 



X — L 



X o 



-y' 



yx^ + y 



z = L{x^ ~ y^ 



(111) 



We consider a superfluid film coating only a circle of 
radius A about the origin of the X,y plane because the 
complete surface has self-intersections. This surface can 
be called the Enneper disk and is illustrated in Fig. [M] 
The figure illustrates that the left and right hand sides of 
the saddle fold over it and would pass through each other 
if allowed to extend further while the front and the back 
would eventually intersect each other underneath the 
saddle. The former pair of intersection curves correspond 
to the two branches of the hyperbola X^ — 3{y^ + !)■ 
When the reference surface is curved into the Enneper 
surface, the X axis bends upward so that the branches 
map to the same intersection curve in the yz plane. (The 
other intersection curves are obtained by exchanging X 
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FIG. 24 Vortex and its streamlines on an "Enneper Disk" 



and 3^.) Since the points where these hyperbolae are clos- 
est to the origin are (±-\/3, 0), (0, ±\/3), a non-self inter- 
secting portion of the Enneper surface results as long as 

Now we explicitly calculate how a single vortex inter- 
acts with the curvature of such a surface by using Eq. 
(IM|) . (We will use conformal mapping instead of Eq. 
(IM)) since the latter equation does not hold on a surface 
with the topology of a disk, as w does not satisfy the 
Dirichlet boundary conditions which are implied by such 
an expression.) The metric obtained from piip is given 
by 

dx^ + dy^ + dz^ = L^{dX^ + dy^){l + X^ +y^f. (112) 

(Surprisingly, this metric is rotationally symmetric. This 
implies that the surface may be slid along itself with- 
out stretching, but with changing amounts of bending.) 
Hence 



^Enneper — — In L{1 +TZ ), 



(113) 



where 7^^ = X"^ +y'^. According to Eq. dM]), this indi- 
cates that the vortex should be attracted to the middle 
of the surface, but of course this force competes with the 



boundary interaction Kt: In 



A^-TZ^ 



which tries to pull the 



vortex to the edge. This expression for the boundary in- 
teraction is obtained from the familiar formula for the 
energy of a vortex interacting w ith its image in a flat ref- 
erence diskl Vitelli and Nelsonl . |2004 ). The total energy 
is then 



E = K[Khi 
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7^1n(A2-7^2)(l+7^2)]. (114) 



As long as j4 > 1, the central point of the saddle is a 
local minimum and this condition is compatible with the 
requirement A < vS for non-self-intersecting disks. Fig. 



[24] shows the flow lines of a vortex forced by the geomet- 
ric interactions towards the center of an Enneper surface 
with A = 1.5. 

In general, conformal mapping allows us to express the 
energy of a single vortex on a deformed surface with a 
boundary in the form: 

E^TTKiln '±L±-u;{u)], (115) 

aA 

where Tl{u.) refers to the image of a defect at u under a 
conformal map to a flat, circular disk of radius A. The 
Green's function method cannot be used to determine 
the energy of defects on a surface with a boundary. Al- 
though the conformal factor w(u) satisfies the Poisson 
equation, Eq. ^, it cannot be expressed as the integral 
of the curvature times the Green's function (as in Eq. 
dnH), since cj does not satisfy simple boundary condi- 
tions. In any case, the first term in Eq. (|115p has no 
general expression in terms of lo either. Interestingly, the 
total single-particle energy satisfies a nonlinear differen- 
tial equation (the Liouville equation): 



V^^(u) = -7ri^G(u) 
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This result can be derived by using Eq. (|105p to cal- 
culate the Laplacian of the first term and using V'^ = 
g2ii;(u) 1 a -^ a |-Q calculate the Laplacian of the second 
term. E{u) also satisfies an asymptotic boundary condi- 
tion: 
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as d, the distance from u to the boundary, approaches 0. 
Together the differential equation and the boundary con- 
dition should determine the total geometrical and bound- 
ary energy of a single vortex, although the nonlinear Eq. 
(|116p is difficult to solve. 



C. Periodic surfaces 

In this section, we illustrate how a periodically curved 
substrate distorts the flat space interaction energies be- 
tween vortices, besides generating the single-particle ge- 
ometric potential. This effect is shown to be a conse- 
quence of the action of a conformal map which generally 
will map the target surface into a periodic reference sub- 
strate with different lattice vectors from the vectors of the 
target substrate. According to the general relation Eq. 
((92)) . the long-distance behavior of the Green's function 
is given by the logarithm of a distorted distance. 

Consider a surface with a periodic height function z = 
h{x,y), i.e., say h satisfies 

h{x + \„y + n,) = h{x,y) fori = {1,2}, (118) 

where i labels the two basis vectors, which are not as- 
sumed to be orthogonal. Figure (|^5|) shows the corre- 
sponding periods (Ai,/ii). A conformal mapping can be 
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FIG. 25 An illustration of the period lattice of the target 
surface, projected into the xy-plane. The coordinates Ai,/ii 
and A2, /i2 describe the basis vectors, which we do not assume 
to be orthogonal. 



chosen to preserve the fact that the substrate is periodic 
but not the actual values of the periods, which are there- 
fore given on the reference substrate by two new pairs 
denoted by (Ai,Mi). In other words, we suppose that a 
tesselation of the target substrate by congruent unit cells 
is mapped to a set of congruent unit cells on the reference 
substrate. Then the map transforming the original coor- 
dinates (x, y) into the target coordinates (<¥, 3^) satisfies 



X{x + Ai, y + ^i) = X{x, y) + A^ 
y{x + A,, y + iii) = y{x, y) + M,. 



(119) 



There is no simple formula for the new set of primi- 
tive lattice vectors {Ai,Mi) for the reference space. In 
some cases, though, precise information can be derived 
from the fact that the (Aj,Mj) share the symmetry of 
the topography of the original substrate. For example, if 
the lattice is composed of bumps which have a 90° rota- 
tional point symmetry, then the reference lattice will be 
square. On the other hand, the topology of the periodic 
surface with a square lattice shown in the contour plot 
of Fig. ([26|) does not posses a 90° rotational symmetry 
and hence its conformal image will have a rectangular 
lattice. 

To get an idea how the conformal mapping behaves 
macroscopically, we try to decompose it into a linear 
transformation C with matrix coefficients {A, B; C, D} 
and a periodic modulation captured by the functions 
£,{x, y) and ri{x, y) that distort the X and y axes respec- 



X{x,y) = Ax + By + i{x,y) 
y{x, y) ^ Cx + Dy + ■q{x, y). 



(120) 



(This decomposition is justified by the self-consistency 
of the following calculations.) The matrix coefficients of 
the linear transformation can be determined by requir- 
ing consistency with Eq. (|119p . Start by evaluating the 
left hand sides of the two Equations (|120p at the posi- 
tions {x + \i,y + ^i} which are shifted by the two pairs 
of periods {Ai,^i}, so that the right hand sides become 
X{x,y) + Aj and y{x,y) + M,, according to Eq. (|119p . 
Then subtract the resulting equations from the corre- 
sponding unshifted Equations (|120p . for each value of i. 
We then obtain two pairs of equations 



AX, 
CX, 









for z = {1,2}, 



(121) 



where we have used the fact that the periodic functions 
S,{x,y) and r]{x,y) are unchanged when shifted by the 
periods. We can now solve the four equations of Eq. 
(|12ip simultaneously for A, B, C, and D to see that the 
linear transformation matrix C reads 



C = 



A B 
C D 



Ai A2 
Ml M2 



Xi A2 

Ml M2 



(122) 



(Now we can justify the original decomposition, Eq. 
(|120p . by defining C by Eq. (jl22[) and defining £,{x,y) 
and ri{x,y) as the discrepancy between the conformal 
map X{x,y),y{x,y) and the linear map C{x,y) as in 
Eq. p20p . We can then check that £,{x,y) and ri{x,y) 
are periodic functions of the coordinates.) 

The linear transformation can be used to calculate ap- 
proximately the long distance behavior of the Green's 
function 



^{x,y;x',y') 



-^HiAxf + iAyr] 

— -InUAAx + BAy)^ 
An 



+ {CAx + DAy)^], (123) 

where we used the fact that the periodic functions £,{x,y) 
and 77(2:, y) are bounded and hence negligible in com- 
parison to Aa; and Ay at long distances. This expres- 
sion illustrates the fact that the matrix C captures the 
long distance lattice distortions induced by the conformal 
mapping, apart from the additional waviness described 
by ^(x,y) and r]{x,y). The linear transformation deter- 
mined by C is by itself typically not conformal, meaning 
that it generates an anisotropic deformation of the tar- 
get lattice which does not preserve the angle between the 
original lattice vectors. 

The deformation of the lattice is controlled by the 
curvature of the substrate. To spell out this connec- 
tion and allow an explicit evaluation of the long-distance 
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FIG. 26 A periodic surface with a square lattice but not 
square symmetry. This surface is iUustrated by its contour 
plot. It has the two reflectional symmetries but not the 90° 
symmetry of a square lattice; hence, generically, the confor- 
mal image will have a rectangular lattice. 



Green's function in Eq. (|123p . we explicitly evaluate 
the matrix elements Ly in terms of the height function 
h(x^ y) of a gently curved (or low-aspect-ratio) surface, 
one for which h{x,y) < {A.,,/iJ, {S,{x,y),'q{x,y)} < 1 
and {A - 1, D - 1, S, C} < 1. The new set of (isother- 
mal) coordinates X and y , used to implement the con- 
formal transformation, are found by solving the Cauchy- 
Riemann Equations (jA3|) 



dyX = ^^g^^'d^y - ^g^^'dyy, 



which, upon substituting from Eq. ()120|1 and making the 
small aspect ratio approximation discussed in Appendix 
\X[ reduce to 

A + d,(^D + dyrj+^{hl~hl) (124) 

B + dy£_^-C-d^r] + h^hy. (125) 

We now proceed to show that these equations do not 
have solutions unless the lattice is distorted, that is to say 
the matrix C cannot be the identity for a generic periodic 
function h{x, y). Note that the periodicity of ^(x, y) and 
•q{x,y) implies that the integral of either one over any 
unit cell, e.g.. 
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is independent of the quantities (a, b) by which the unit 
cell is shifted. Upon differentiating the integral with re- 
spect to a, one obtains 

// dxdyUx,y)=0. (127) 

J Jccll 



Similarly, the averages of rj^jriy, and ^j, are also equal to 
zero. Hence, upon averaging Eq. (|125p over a unit cell 
we obtain the key relations 



A-D = -<hl-hl> 

2 X y 

B + C = < h^hy >, 



(128) 



which prove our assertion that L cannot be a simple di- 
lation or rotation. Some shear is naturally introduced by 
the non trivial metric (or height function) of the under- 
lying surface. 

The C matrix is undetermined up to a dilation and a 
rotation but this is of no consequence to the determina- 
tion of the Green's function. In fact, to find the Green's 
function, note that Eq. (|128p allows us to write the ma- 
trix coefficients in terms of two undetermined constants 
ei and £2 that will drop out of the final answer: 
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so that consistency with Equations (|128p is guaranteed. 
(The variables ei and €2 parameterize an overall in- 
finitesimal scaling (by 1 -t- ei) and a rotation (by angle 
£2) respectively.) Substitution of these equations into 
Eq. (|123p gives the desired long-distance behavior of 
the Green's function purely in terms of derivatives of the 
height function, which we assume to be known: 
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(133) 



This is the central result of this section; it can also be ap- 
plied toJiitera£U£ns_betwreendisclinations in liquid crys- 
tals ( Vitelli and Nelsonl. |2004| ) and dislocations in crys- 
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tals () Vitelli et aLl . 120061 ). The anisotropic correction to 



the Green's function, captured by the second and third 
term, suggests that a distorted version of the triangular 
lattice of vortices expected on a flat substrate may form 
when the helium-coated surface is rotated slowly enough 
that there is only one vortex to several unit cells. How- 
ever, the actual ground state is likely to be difficult to 
observe, as the geometric potential will try to trap the 
vortices near saddles as discussed in Section Hi. Bl 



D. Band-flows on elongated surfaces 

In this section, we show that the quantization of cir- 
culation can induce an extremely long-range force on a 
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FIG. 27 A capped cylinder; a cylinder of length 2H is closed 
off by hemispheres at the north and south poles of radius R. 
The circulation around every lattitude is the same. 



stretched-out sphere (such as a surface with the shape of 
a zucchini or a very prolate spheroid). We first demon- 
strate the main result in the context of a simple example 
before presenting a general formula for the forces expe- 
rienced by vortices on azimuthally symmetric surfaces. 
Details are presented in Appendix [Ej Consider a cylin- 
der of length 2H and radius R « H with hemispherical 
caps of radius R at the ends, depicted in Fig. [271 and 
imagine a symmetric arrangement of a vortex (n=l) and 
an anti-vortex (n=-l) at the north and south poles re- 
spectively. Extrapolating our intuition from flat space 
suggests that the energy of the vortex and anti- vortex is 
2-kK\\i — , where D is the distance between the vortices. 
However, more careful reasoning shows that the energy 
grows linearly rather than logarithmically with D. The 
reason is that, unlike in flat space, the velocity field does 
not fall off like the inverse of the distance from each vor- 
tex. Note that the azimuthal symmetry of the arrange- 
ment of the vortices implies that the flow is parallel to 
the lines of latitude of the surface. Since the circulation 
around each latitude must be — , the flow speed on the 
cylindrical part reads 
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The kinetic energy of this part of the flow is 
[4:TrRH]^Psv'^ = 27r-^^. Since this cylindrical part of 
the flow forms the main contribution to the kinetic en- 
ergy when H >> R we find that the energy of a vortex- 
antivortex pair situated at opposite poles is linear, 
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(The exact expression also includes a near- vortex energy 
of approximately 2nK\n —.) 

In contrast, when the vortices forming the neutral pair 
are across from each other on the same latitude, the afore- 
mentioned long-range persistence of the velocity field is 
absent because the vorticity is screened within a distance 
of order R. The resulting kinetic energy follows the fa- 
miliar logarithmic growth 



9 p 

Eequator ~ 271 K In . 

a 



(136) 



More generally, consider an azimuthally symmetric 
surface described by the radial distance, r{z), as a func- 
tion of height, z, as indicated in Fig. [^5] A. If the north 
and south poles of the surface are at Zs and z„, then 
r{zs) = r{zn) = since the surface closes at the top 
and bottom. A point on such a surface can be identi- 
fied by the coordinates {(j), a) where </> is the azimuthal 
angle and a is the distance to the point from the north 
pole along one of the longitudes such as the one shown in 
Fig. 1281 'V. In Appendix |E1 we develop an approximation 
scheme which rests on the observation that the flow pat- 
tern becomes mostly azimuthally symmetric if ^ << 1, 
even if the vortices break the azimuthal symmetry of the 
surface because they are not at the poles. If a pair of 
ji ~ ±1 vortices are present at different heights zi^2, then 
the fluid in the band between them flows almost horizon- 
tally and at a nearly ^-independent speed (except for 
irregularities near the vortices) while the fluid beyond 
them is approximately stagnant (see Fig. [28]) . Along 
any latitude inside the band the circulation is exactly — , 
while it is zero above and below it. These properties ap- 
proximately determine the flow away from the vortices 
since the asymmetric irregularities near the vortices de- 
cay exponentially, giving a speed of 
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in between zi and Z2, the locations of the vortices, and 
zero elsewhere. We describe this approximation as the 
"band model" . This expression shows that constrictions 
in the surface cause the speed to increase. For the cylin- 
der with spherical caps and arbitrarily placed vortices, 
Eq. (|137p shows that the speed is approximately con- 
stant within the band. (It increases within a distance on 
the order of R from the vortices, which are on the edges 
of the band.) 

The kinetic energy can be determined approximately 
by noting that the energy in a thin ring on the surface 
between the vortices (extending from the longitudinal ar- 
clength a to a + da) is 
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where the flrst factor represents the area of the ring since 
a is the geodesic distance along the surface and the sec- 
ond factor represents the included kinetic energy. The 
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On the capped cylinder, this force is independent of the 
positions of the vortices. Even on an arbitrary elongated 
surface, a noteworthy feature is that the force on vortex 
1 does not depend on the position of vortex 2! This force 
can be explained with the familiar phenomenon of lift: 
the vortex is on the boundary between stationary and 
moving fluid, so there is a pressure difference due to the 
Bernoulli principle. 

Approximating the flow pattern generated by multiple 
vortices in a similar fashion requires only minor modi- 
fications of the previous argument. In a low-resolution 
snapshot of the flow, the point-vortices would appear as 
circles of discontinuity in the velocity field that go all the 
way around the axis (the analogue for a layer of super- 
fiuid of a two dimensional vortex sheet). If the vortices 
are labeled in order of decreasing z a loop just below the 
^th vortex contains 
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units of circulation above it. Approximate azimuthal 
symmetry of the flow then implies that, 
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FIG. 28 The flow on an azimuthally symmetric surface, de- 
scribed by the coordinates (z, r) and an azimuthal angle (j) 
(not shown). A) The surface is defined as the surface of rev- 
olution of the curve r{z) in the r — z plane. The other two 
images show the flow on an ellipsoid and compare the flow 
pattern predicted by the band model (B) and the exact solu- 
tion determined by conformal mapping (C). The flow lines in 
the band between the two vortices become close to horizontal 
and are approximately azimuthally symmetric. Beyond the 
vortices, they are spaced far apart indicating a vanishingly 
small speed for a greatly elongated surface. 



flow is zero past the two vortices, so the total energy is 
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For the capped cylinder with a constant r(z) the integral 
is TT-^y^, where D is the vertical distance between the vor- 
tices. This expression generalizes the result for vortices 
at the poles of the surface, Eq. (|135p . However, when the 
vortices are on opposite ends of an equator, they are too 
close for the asymmetries to be neglected. The energy in 
this case is calculated in Appendix[El The force Fi^band, 
experienced by vortex 1 can be determined by taking the 
gradient of Eq. (|138p . If vortex 1 is moved downward, 
the band of moving fluid shrinks so the energy drops. 
The force on the vortex therefore reads 
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4> (for zi< z < zi+i), (141) 



a natural generalization of Eq. p37p that is proved in 
Appendix [E] 

Conformal mapping can be employed to justify (with- 
out detailed calculations) the decay of the nonaz- 
imuthally symmetric parts of the flow that are not de- 
termined by the quantization condition. We sketch the 
basic reasoning here by focusing (for simplicity) on the 
flow pattern near the equator of the surface, at a distance 
CTeg from the north pole. The conformal transformation 
that maps the elongated sphere onto a regular reference 
sphere with coordinates 9,$ reads (see Appendix [E]) 
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The upper and lower halves of the elongated sphere can 
be mapped to the upper and lower hemispheres by choos- 
ing appropriately between the two values of 6 that corre- 
spond to a given value of sin 0. Near the equator the inte- 
gral can be approximated by ^ since r varies slowly. 

Teg 

Suppose the vortices are far from the equator, at a dis- 
tance greater than kreq for a large k. Then the vortices 
above the equator are mapped exponentially close (at a 
distance less than e^'^) to the sphere's north pole. Like- 
wise vortices on the southern half of the surface map 
exponentially close to the sphere's south pole. We have 
thus reduced the task of finding the flow due to a com- 
plicated arrangement of vortices to a symmetric case. In 
fact, after mapping the flow on the long, thin surface to 
the reference sphere, nothing can be resolved beyond a 
pair of multiply quantized vortices at the north and south 
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poles containing N and —TV units of circulation respec- 
tively, where N is the total circulation number of all the 
vortices above the equator. Since the image vortices are 
very close to the poles, their flow pattern on the reference 
surface is approximately azimuthally symmetric near the 
equator. When mapped back to the elongated surface, 
the flow retains its approximate azimuthal symmetry in 
the region around the equator, completing our argument. 
A similar argument proves the approximate azimuthal 
symmetry of the flow near lines of latitude other than 
the equator; one simply adjusts the conformal map in 
Eq. (|142p so that another latitude of the target surface 
is mapped to the equator of the reference sphere. 

Now the geometrical force derived in Eq. ()139p can 
compete with physical forces such as those induced by 
rotating an ellipsoid about its long axis with angular ve- 
locity ft = fl z. Let us extend the treatment of rotational 
forces on curved substrates, introduced in section [TTTl to 
the case of an ellipsoid described by the radial function 
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Let us use the aspect ratio a = -^ to describe how elon- 
gated this ellipsoid is, and determine how a vortex-anti- 
vortex pair is torn apart by the rotation as the angular 
frequency is increased. As in Section [1111 metastable vor- 
tex configurations can often be found, so we will consider 
transitions between different local minima of the vortex- 
energy function. Recall that the effect of the rotation on 
the superfluid energy is expressed in terms of the angular 
momentum L^ by an extra term — fiL^ in Eq. (|40[) which 
must be evaluated for (at least) a pair of opposite-signed 
vortices to satisfy the topological constraints imposed by 
the spherical topology of the surface. We find, in anal- 
ogy to Eq. (|45)l . that the extra contribution to the vortex 
energetics is additive and reads 
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where A{a) represents the area of the ellipsoid out to a 
distance a from the north pole while ai and CT2 represent 
the positions of the two defects (cti < tT2)- If ni = 1, the 
energy in Eq. (|40p is decreased by moving the positive 
vortex closer to the north pole and the negative one closer 
to the south pole. To understand this, note that the sense 
of rotation of the superfluid around a vortex is defined by 
an observer facing the surface. Hence, a negative vortex 
on the southern half of the surface rotates in the same 
direction as a positive one on the northern half (relative 
to the positive z-axis), and both agree with the sense of 
rotation of the substrate. The rotational force on vortex 
1 derived from Eq. (PP)) is 
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As 57 increases, pairs of positive and negative vortices 
will appear in this geometry. As each vortex pair is cre- 
ated, the positive vortex will move to the top side of the 



surface, and the negative one to the bottom. There is a 
critical frequency, 17;,, at which a single pair of vortices, 
once created, can exist metastably in a configuration 
symmetric about the oiy-plane. As the angular frequency 
is increased, the vortices are gradually pulled apart until 
at the frequency f7a, they reach the poles. (The latter 
transition is analogous to the center/off-center transition 
for a single vortex described in Section IIILBI ) When 
57a > 57 > f7b, the equilibrium condition, obtained by 
balancing the forces in Equations (|139p and (I145P , reads 
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Since r and z are connected by the equation of the ellip- 
soid, the vortices are located at heights 
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The vortices first become metastable when force bal- 
ance is achieved with both vortices close to the equator. 
Upon substituting the equatorial value ri — R into Eq. 
(|146p an estimate of fib is obtained 
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When the pair first appears, there will actually be a non- 
zero defect separation, although substituting Eq. (|148p 
into Eq. (|147p suggests otherwise. Imagine slowing 
the rotation speed through 57b . The vortices will ap- 
proach each other gradually; within the large vortex- 
separation approximation of Eq. (|139p . the attraction 
between them will decrease as they become closer be- 
cause r(z) increases. However, when the vortices become 
close enough, the attraction between them starts increas- 
ing and the vortices are suddenly pulled together. The 
minimum z-coordinate for metastable vortices is derived 
along these lines in Appendix |E] (which also discusses 
what happens at 57 = 57a) and reads 
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The transition through 57^ is illustrated pictorially in Fig. 
1291 which shows how the local minimum in the energy 
function disappears as the frequency decreases. 



E. Interactions on a closed surface 

To understand interactions between vortices on an ar- 
bitrary deformed sphere one must come to terms with the 
neutrality constraint on the total circulation of a flow. 
On any compact surface. 
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This constraint on the sum of the circulation indices {nt} 
always holds: if the surface is divided into two pieces by 
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FIG. 29 The rotational and fluid energy (units of K) as a 
function of cti = atot — o"2 = o" (units of R) for H — 3.5R 
and ^ = .49, .61,.74_R"^. The middle curve, roughly at 
u — ujt, shows the last position where the vortex is stable as 
u is decreased. 



a curve, the sum of the quantum numbers on the top and 
bottom half must be equal and opposite (because they are 
both equal to the circulation around the dividing curve) . 
As we shall see, this relation implies that there are mul- 
tiple ways of splitting up the energy into single-particle 
energies and two-particle interaction energies, despite the 
fact that the total energy is well-defined. The behavior 
of the one-particle and interaction terms depends on how 
the splitting is carried out. To illustrate this ambiguity, 
multiply Eq. (|150p by 4:n'^Knif{ui) and separating out 
the i = 1 term, to obtain 

4n^nlf{ui) = - ^ Air^mnJiui). (151) 

Hence, a portion 4:TT'^Kf{ui) of the "geometrical energy" 
of vortex 1 can be reattributed to this vortex's interaction 
with all the other vortices. This can be seen explicitly by 
checking that the net energy according to Eqs. (PTjl and 
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is not changed by the following transformation: 



r'(ui, U2) = r(ui, U2) - /(Ui) - /(U2) (154) 

[/^(u) = [/G(u)+47r/(u) (155) 

This flexibility is reflected in the possibility of choos- 
ing different Green's functions r'(ui,U2) for the covari- 
ant Laplacian on a deformed compact surface. A detailed 
discussion of Green's functions is given in Appendix |D] 
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FIG. 30 A) The process of folding a hemisphere into a pointed 
sphere, bounded at the north and south poles by two 180° 
disclinations. The north and south poles move outward along 
the axis, while the latitudes stay horizontal. The Gaussian 
curvature is invariant because the decreasing curvature of the 
lines of longitude is compensated by the tighter curvature 
around the lines of lattitude. B) A top view of vortices dur- 
ing the furling-up of one hemisphere. The first stage shows 
both the (0 < a < tt) hemisphere that is furled up and the 
other hemisphere together with two defects and their images. 
The hemisphere is furled up into the pointed sphere so that 
the left and right halves of the cut (which appears as a hori- 
zontal diameter in this top view) are brought together to form 
the seam on the pointed sphere; simultaneously, the defects 
Qi,Q2 move to positions ui and U2 on the pointed sphere. 
The furling process leads to a continuous flow pattern on the 
pointed sphere. For example, the two points marked with an 
X on the cut in the original sphere are sealed together. Both 
points feel a strong flow, the one on the right because it is 
close to the vortex at Qi and the one on the left because it is 
close to this vortex's image. 



Here we higlight this ambiguity by performing explicit 
calculations using two distinct choices of Green's func- 
tions on a model surface formed from a unit sphere. First 
cut the sphere in halves along a great circle. Choose one 
of the hemispheres and bring opposite sides of the great 
circle bounding it together and glue them. The result is 
a pointed sphere resembling a conchigliette noodle (i.e., 
a shell noodle) sealed shut (see Fig. ISOK). The surface 
closes up smoothly since opposite sides of the seal have 
matching curvatures; the surface also turns out to be 
rotationally symmetric. (A similar shape forms when a 
pollen grain w ith a weak sect or (such as the pollen of a 
lily) dries out (JKatifori et alh .) 

The geometrical interaction on this surface can be 
given an appealing interpretation in terms of the method 
of images from electrostatics. When one uses the method 
of images to study charges in a half-space bounded by 
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a conducting plane, one completes the space, with the 
other half-space. Then one introduces charges into this 
fictitious region to ensure that the right boundary con- 
ditions (orthogonality of the field lines to the original 
boundary plane) are satisfied. For the pointed sphere, 
we first complete the surface by opening it up again and 
adding back the second hemisphere. We can describe 
points on the pointed sphere by {(f), a) (where a is the 
distance from the north point and cj) is the azimuthal 
angle, as for the general case described in Sec. IV. Pp . 
Points on the completed sphere are naturally labeled by 
the standard spherical coordinates (a, 9) that represent 
the azimuthal and polar angles respectively. The map- 
ping from the pointed sphere that returns each point to 
its position on the original sphere is given by 
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As (p ranges from to 2tt, a goes from to tt, across the 
hemisphere used to make the pointed sphere. Note that 
= a since the pointed sphere is formed by isometri- 
cally bending the hemisphere (the angle 9 is equal to the 
geodesic distance to the north pole of the hemisphere) . 

Now for each vortex at u = (</>, (t) on the pointed 
sphere, we introduce two vortices on the sphere, one at 



Q with (a, 



(|,cr) and one at Q* 



f -f 7r,CT). The 
latter is the image vortex of the former, obtained by ro- 
tating Q by 180° around the z-axis. 

The energy of a defect configuration on the pointed 
sphere is derived by halving the energy of the flow pattern 
produced by the doubled set of vortices on the full sphere. 
In analogy with the electrostatic problem, the purpose 
of situating the image defects in the way just described 
is to preserve the continuity of flows across the seam. 
Imagine drawing the flow pattern of all the vortices on 
the sphere. Focus on the hemisphere < a < tt. Because 
the vortices are placed symmetrically about the sphere's 
axis, the flow near a = will match the flow near a — tt 
when the surface is sealed. (See Fig. [HOB.) 

The flow pattern on the pointed sphere results from 
rolling up half of the flow pattern on the sphere. Once 
the positions of the image defects are chosen, the flow 
pattern is found by deriving it from the stream function 
x(u) introduced in Sec. |TT1 The stream function at a 
point u on the pointed sphere can be expressed in terms 
of the Green's function of the sphere, according to Eq. 
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where Q is the point on the sphere corresponding to u 
on the pointed sphere and Qi and Q* are the locations 
of the i^^ pair of vortices on the sphere. 

The energy of the vortices on the sphere takes up the 
familiar electrostatic form of a sum of the interactions 
between all pairs of defects and/or their images. The 



energy stored in the flow pattern on the pointed sphere 
(which is half as large as on the complete sphere) reads 
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In the second expression, we note that the terms in the 
first line are equal in pairs, so that a factor of ^ cancels. 
This energy is given the same form as Eq. ()153|) by sep- 
arating out the part which depends on the positions of 
two vortices, proportional to 

r,(ui, U2) = r^p''"'^^(Qi, Q2) + r^p^^'-^(Qi, q;) (i59) 

and the part which depends on one vortex, 

C/,(u) = -27rr^P"^'-^ (Q,Q*). (160) 

The function in Eq. (|159p is a Green's function for the 
pointed sphere, as the placement of the images guar- 
antees that this function is well-defined on the pointed 
sphere. It appears in the stream function, Eq. (|157p . 
as well as in the energetics, as expected for a Green's 
function. 

The potential which describes the single-particle en- 
ergy of a vortex becomes singular as the vortex ap- 
proaches the apex of the cone at the north or south pole, 
since then the vortex Q approaches its image Q* . This 
is in accord with the result, Eq. (|9ip . that the Gaussian 
curvature is the source of the single-particle energy since 
the pointed sphere has delta- function concentrations of 
curvature at its north and south poles: 



G(u) = 1 + nSNiu) + TrSsiu). 
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where <5Ar(u) and (55(u) are the appropriate delta fun- 
nctions. The geometric repulsion from the positive cur- 
vature points arises from the repulsion between vortices 
and their images! We can check step- by-step that Us is 
sourced by the Gaussian curvature, 

C/,(u) =. - /"/"r,(u,u')G(u')d'u'. (162) 

We substitute for G(u') from Eq. (|16ip and for F^ from 
Eq. (|159p which can be written in the form, 

r,(ai, (/-i; ^2, (/.2) = F^P''--(ai, ^; a2, ^) 
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In the last line, we have evaluated the Green's function 
for the sphere by writing the chordal distance between, 
e.g., Qi and Q2 in terms of the spherical coordinates 
(cti.2,01.2 ), D^ = 2[1 — cos^i cos 02 — sin 01 sin02cos(ai — 
02)] (see (JLubenskv and Prostl . ll993 )). and then combin- 
ing the two terms together. To evaluate the integral in 
Eq. (I162|) , we have to note that the area element of this 
integral is d^u = i sin adadcf). The area of a region on the 
pointed sphere is the same as the area svnOdOda of the 
corresponding region on the original sphere, and the fac- 
tor of ^ results from how the angles are related, a = ^, 
see Eq. (|156l) . Now the integral on the right-hand side 
of Eq. (|162p can be shown to be equal to the left-hand 
side using the identities 
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= 27rln— if B > A. 
2 

We have now derived one formulation of the energetics 
in terms of T g and t/, , the corresponding geometric po- 
tential. Let us contrast this isometric mapping method 
with the conformal mapping method in order to illustrate 
how different approaches can naturally lead to different 
delineations between vortex-vortex and vortex-curvature 
interactions. (The net result is of course the same from ei- 
ther point of view.) As a result of the isometric mapping 
each point is doubled, whereas the distance-distorting 
conformal mapping transforms each point on the pointed 
sphere to one point on the reference sphere. 

We first use Eq. (|142p to find that the conformal map 
is given by 
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Comparing the conformal mapping results, Eqs. (j88p and 
([55)1 to the Green's function formulation, Eqs. (pij) and 
([90)1 suggests the following identification of the interac- 
tion potential (or Green's function) and single-particle 
potential: 

r,(ui, U2) = r«p''^'^^(e(ai), <j,r, 0(^2), cj,^) 

, ^ 2 sin (J 

1 + cos^ a 



These expressions differ from Equations (|159|) and (|160p . 
Nevertheless, as promised, the net energy is the same 
whether the pairs {Ys-,Vs) or (Fc, t/c) are used in place 
of r and Uq- In fact, 



rc(ui, U2) = rs(ui, U2) - /(ui) 

C/,(u)-C/,(u)+47r/(u) 
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where /(u) = — ^ln(l -I- cosher). This transforms the 
energy from the single-particle to the interaction terms 
consistently as described at the beginning of the section. 
Appendix [P] shows that the Green's function formulation 



is generally equivalent to the conformal mapping result 
derived in Section IV. A[ even when there is no method of 
images that can be used to determine the Green's func- 
tion explicitly in general. 



VI. LIMITS ON THE STRENGTH AND RANGE OF 
GEOMETRICAL FORCES 

Geometrical forces are limited in strength due to the 
nonlinear relation between the curvature and the geo- 
metric potential. Curvature affects both the source of 
the geometrical force and the force law, as illustrated in 
the examples of Sees. IV. CI and IV. Dl As a consequence, 
even on a wildly distorted surface (with planar topology), 
there is a precise limit on the strength of the force on a 
single vortex. This result has the character of a geomet- 
rical optimization problem, like maximizing the capaci- 
tance of a solid when the surface area is given. Consider 
a vortex located at the center of a geodesic disk of radius 
R. Assume that the Gaussian curvature is zero within 
the disk, but may be different from zero elsewhere. Then 
the force F due to the curvature satisfies 
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where ni is the number of circulation quanta in the vor- 
tex. This relation between R and F is proven in Ap- 
pendix |f1 

If one warps a surface in a vain attempt to overcome 
the hmit, but the force gets diluted because the distortion 
of the region around the curvature pulls the force-lines 
apart, as we can understand from the simple example of 
vortices on cones. 

A cone of cone-angle 9 is obtained by taking a segment 
of paper with an angle 9 and gluing the opposite edges of 
the angle together. This is most familiar when 9 < 2tt. 
li 9 — 2'Km + (3, such a cone can be produced by adding 
TO extra sheets of paper, as illustrated in Fig. [311 We 
slit the m sheets of paper and put them together with 
an angle of size f3 cut out of an additional sheet. By 
gluing the edges of the slits together cyclically, a cone of 
arbitrary angle 9 is made. 

A cone has a delta function of curvature at its apex, but 
no Gaussian curvature elsewhere because the surface can 
be formed from a flat piece of paper without stretching. 
The weight 211 — 9 of the delta function is expressed, 
according to the Gauss-Bonnet theorem, as an integral 
of th e Gaussian cur vature in any region containing the 
apexdKaimenl, |2002| ) 
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where k is the geodesic curvature along the boundary of 
the region and s its arc length. Apply this formula to 
the circle of radius D centered at the apex of the cone. 
Imagine the circle as it would appear on the original 
sheets of paper, as in Fig. 1311 Its measure in radians 
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FIG. 31 How to form cones of negative curvature. One com- 
plete sheet of paper is slit and an angle is cut out of an addi- 
tional sheet of paper. The edges labeled 1 are taped together, 
and then the edges labeled 2 are taped. The circular arcs 
join together to form an extra-large circle. The cone angle is 
9 = 271 + (3, and cones with even larger cone angles can be 
formed by using additional sheets. 



is l3 + 2Trm = 9 since it consists of m complete circles 
together with an additional arc. The length is therefore 
S = DO. The geodesic curvature of the circle does not 
change when the cone is unfolded, so it is equal to -k. 
Upon substituting in Eq. (|168p . we obtain 



G{u)d?u. ^2tt- S— = 2iT-e 
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When 9 > 211 the curvature is negative. 

Now imagine a vortex (with ni = ±1, say) at a dis- 
tance D from the cone point, on the circle of circum- 
ference S just considered. The arbitrarily large negative 
curvature which is possible by making m large seems to 
defy the general upper bound on the geometric force. 
According to Newton's theorem, applied to the radius D 
circle centered at the cone's apex and passing through 

the vortex, the force on the vortex is F = - — ^ -. 

Since the circumference S — D9 is larger than it would 
be in the plane, the force is diluted; substituting the inte- 
grated curvature from Eq. (|169p . we find that it is given 
by 



F = 7r 



K2'K~9 
D 9 ■ 



(170) 



This satisfies Eq. (|167p for all negatively curved cones 
{9 > 2n) ; even when 9 ^ oo the magnitude of the force 
is less than 47r^ because the large circumference in the 
denominator of the Newton's theorem expression cancels 
the large integrated curvature in the numerator. 

In the opposite limit 9^0, the theorem described by 
Eq. (|167p is still correct of course. One has to be careful 
about applying it, however. The force on a vortex at 
radius D (given by Eq. dTTD])) is not bounded by ^ 
with R set equal to D when 9 is small enough (in fact, 
for an extremely pointed cone, <C 1, the force given 
by Eq. (|170p diverges), but this does not contradict the 
inequality because the circle of radius D centered at the 



vortex is pathological: although it does not contain any 
curvature, the circle wraps around the cone and intersects 
itself. Taking R to be the radius of the largest circle 
centered at the defect which does not intersect itself, one 
finds that the inequality is satisfied, with room to spare, 
for all values of the cone angle 9 (see Appendix |F| . One 
can describe a more awkwardly shaped surface such that 
the force on a singly-quantized vortex is arbitrarily close 
to the upper bound ^^ (see Appendix IFJI . 

One can also provide limits to the strength of the geo- 
metric force from a localized source of curvature. Rota- 
tionally symmetric surfaces such as the Gaussian bump 
have force fields that do not extend beyond the bump, 
since the net Gaussian curvature is zero, and Newton's 
theorem says that only the integrated Gaussian curva- 
ture can have a long range effect for a rotationally sym- 
metric surface. To get a longer-range force, one must 
focus on non-symmetric surfaces, like the saddle surface 
of Sec. III.BI The integration methods of Appendix |A] 
can be used to show that this surface's potential has 
a quadrupole form at long distance. Let us consider, 
more generally, a plane which is flat except for a non- 
rotationally symmetric deformation confined within ra- 
dius R of the origin. (The result will not apply directly 
to the saddle surface since its curvature extends out to 
infinity.) In this case, the total integrated Gaussian cur- 
vature is zero, implying that the long-range force law can- 
not have any monopole component. A dipole component 
is not ruled out by this simple reasoning, but Appendix 
IfI shows that the limiting form of the potential is at least 
a quadrupole (or a faster decaying field). 



Eir) 



, ^2COS(20- 72) 



(171) 



where r and cf) are the polar coordinates of the vortex 
relative to the origin, and fi2 and 72 are constants that 
depend on the shape of the deformation in the vicinity of 
the origin. As in the previous case, there is an upper limit 
on the quadrupole moment ^2 ; no matter how strong the 
curvature of the deformation is: 



^2 < ttKR^ 



(172) 



For electrostatics in the plane, the maximum quadrupole 
moment of N particles with charge 27r and N with charge 
—2tt in a region of radius R is at most of the order of 
KNR^, which has the same form as the bound in Eq. 
(|172p . except for the factor of N. Because of the non- 
linearity of the geometrical force and restrictions on how 
much positive and negative curvature can be separated 
from each other, the quadrupole moment is bounded no 
matter how drastically curved the surface is. 

These results describe key physical differences (result- 
ing from the fact that the curvature cannot be adjusted 
without changing the surface) between the geometri- 
cal forces discussed in this work and their electrostatic 
counterparts despite the close resemblance from a formal 
viewpoint. 
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VII. CONCLUSION 

In this article, we have laid out a mathematical formal- 
ism based on the method of conformal mapping that al- 
lows one to calculate the energetics of topological defects 
on arbitrary deformed substrates with a focus on appli- 
cations to superfluid helium films. The starting point of 
our approach is the observation that upon a change of co- 
ordinate the metric tensor of a complicated surface can 
be brought in the diagonal form gab = e'^'^^^^Sab- This 
corresponds to the metric of a flat plane which is locally 
stretched or compressed by the conformal factor e^'^*-"'. 
Many of the geometric interactions experienced by topo- 
logical defects on curved surfaces are simply determined 
once the function uj{u) is known. Vortices in thin helium 
layers wetting a curved surface are a natural arena to 
explore this interplay between geometry and physics but 
our approach is of broader applicability. 

The curved geometry results in a modified law for de- 
fect interaction as well as in a one body geometric po- 
tential. On a deformed plane, the latter is obtained by 
solving a covariant Poisson equation with the Gaussian 
curvature as a source. Table [J presents a summary of the 
general form that the defect interaction (first row) and 
the geometric potential (third row) take up in curved 
spaces with the topology of a deformed plane (first col- 
umn), disk (second column) and sphere (third column). 
These results can be derived starting from the differen- 
tial equations that the geometric potential satisfies or the 
appropriate Green's functions that we list in the second 
and fourth row respectively for each of the three surface 
topologies. The fifth row of Table [l] directs the reader to- 
wards the relevant sections and appendices of the paper 
where he will be able to find some concrete applications 
of the formalism and technical derivations. 

For example, the geometric potential of an Enneper 
disk (a minimal surface with negative curvature described 
in Sec. IV. Bp is given by the conformal factor a;(u) eval- 
uated at the point P = {ui,M2} where the vortex is lo- 
cated combined with an "electrostatic-like" interaction 
with an image defect located at the inverse of P with 
respect to the circular boundary. The geometric poten- 
tial satisfies the Liouville (non-linear differential) equa- 
tion that reduces to the Poisson equation derived for the 
plane in the limit of an infinitely large disk. In the case 
of deformed spheres, we showed in Appendix [P] that one 
can make a convenient choice of Green's function so that 
all the geometric effects are included in the defect-defect 
interactions without introducing a one-body geometric 
potential explicitly. An interesting application naturally 
arises on vesicles deformed into an elongated shape, like 
a zucchini. The range of the defect interaction becomes 
much longer and its functional form different from the 
logarithmic dependence expected in flat two dimensional 
spaces. 

We hope that the discussion of the geometric effects 
presented in this work may pave the way for their ob- 
servation in thin superfluid or liquid crystal layers on a 



curved substrate. A useful starting point could be the de- 
sign of experiments to detect the geometric potential by 
balancing it with forces exerted on the defects by exter- 
nal fields or rotation of the sample as discussed in Section 
mil Such experiments should focus on single vortices, or 
on situations where the separation between vortices is 
comparable to the length scale of the geometry. Signa- 
tures of the geometric interactions described here may 
also survive in defect pinning experiments ca rried out in 
some bounded three dimensional geometries ( Voll et all . 
I2006D . 
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APPENDIX A: Nearly Flat Surfaces 



The calculations in Section riI.B | are based on pertur- 
bations about near flatness (see ( Davidl . Il989[ ) and refer- 
ences therein). The perturbation theory will be in powers 
of an aspect ratio, a, which measures the ratio of surface 
height to width of the landscape features. (We imag- 
ine that the height of the surface is given in the form 
am{x,y) where to is a fixed function.) The leading cor- 
rections to the flat space energies are second-order in a. 
There are two of these; one is the geometric potential. 
When there are at least two vortices present, there is 
also a second order correction to the Green's function, 
which ought to be retained since it is comparable to the 
geometric potential. The latter could be calculated by 
expanding the metric in Eq. (J13p in powers of a. Never- 
theless, because the perturbations are singular, we prefer 
to use conformal mapping for this step just as we use 
in Sec. IV. Al to derive the geometric potential. Our cal- 
culations are limited to the case of an infinite deformed 
plane. 

We use the x and y coordinates of a plane parallel 
to the surface for our coordinate system (the "Monge 
Gauge" ) . The metric is then ds"^ — dx^ + dy"^ + dz"^ — 
(I -I- hl)dx'^ + (1 + hl)dy'^ + 2hxhydxdy. Subscripts on h 
indicate derivatives, so that hxx = d^h etc. Upon calcu- 
lating the curvature tensor we find the Gaussian cur va- 
ture in the second order approximation (JDavidl . Il989f ) 



G{x, y) = hxxhyy - hly. (AI) 

The geometric potential is found by approximating the 
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Kn,nj[An^r{ui,Uj)-2TTlnf\ 


A-r"^ KrLifijV D{vii, Uj) 


4Tr^Kn,njrp{u,,Uj) 


Green's 
Function 


1 ^^\u,-u,\ 

Zir a 


1 , (i+n',n]-2n,n,co^(^.i-^,)\ 


-2^1n^ + 4^[^('^':)+'-'K): 


Geometric 

Potential, Bi(ui) 


AV[//rf2ur(u„u)G(u)+lnf] 


TrKln^^-irKLoiu) 
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Applications 
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Sees. II.A-II.C; 
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Appendix A 
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Sec. V.B 
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Neutrality and 

Green's function flexibility: 
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TABLE I An outline of vortex interactions on curved surfaces. The net energy of a set of vortices on a surface with the 
topology of a plane, disk, or sphere is given by J^^ nfi?i(ui) + X]i< E2{ui, Uj), where simple expressions for the single-particle 
(or geometric) potential and two-particle potentials are given in the table. A conformal mapping is necessary for evaluating 
some of these expressions. For example. Hi (in the expression for the Green's function on a deformed plane) is the Cartesian 
coordinates of the conformal image of vortex i. 



exact expression Eq. (|29p . The Green's function may be 
replaced by the planar one since G is already quadratic 



Uoix^y) 



E 



dx'dy'Tfiat{x, y; x' , v'){Kxhy 



"xyJ 



(A2) 



To do a conformal mapping to an equivalent flat space 
problem, we must solve the curv ed space gener alization of 
the Cauchy-Riemann equations (JDavidl . Il989l) which de- 
fine isothermal coordinates, X ,y, namely 



v°(y) = -7b"v''(x). 



(A3) 



where 7^ — g"''^^/gf-cb- (In words, the gradients of X 
and y are at right angles to each other and have equal 
magnitudes at every point.) We insert the expression 
for the metric in terms of h into Eq. (|A3|) and expand 
to second order under the assumption that X = x + 
S^,y ^ y + t] where the deformation parameters ^ and 77 
are second order in a: 



Vx 



Vy-^x 






(A4) 
(A5) 



By taking the derivative of Eq. (|A4p with respect to x 
and Eq. (jASp with respect to y and adding the results, 
we obtain 



^jlatV ■■ 



hy'^jlath, 



(A6) 



which may be solved by means of the Green's function 
(keeping in mind the boundary condition that the con- 
formal map must approach the identity at infinity; i.e., 
?7 ^ 0), giving the result 



JJ dx'dy'TfiatihyWJi^th). (A7) 



vix,y) 



Similarly, we may solve for ^ (and then check that (jA5|) 
is satisfied). 

We can use our expression for the conformal map- 
ping in conjunction with the conformal invariance of 
the Green's function Eq. ([9^ . which implies that 

r{xi,yi;x2,y2) = ^fiat{xi+(,i,yi+r]i;x2+^2,y2+V2)- 
Upon expanding the flat space Green's function in ^ and 
77, we have 



r(x,2/;x',j/')~-^ln(Ax2 + A2/') 



A^Ax + ArjAy 
27r(Ax2 



Ay^y 
(A8) 

where AQ indicates taking the difference between the 
values of Q at the two points where the Green's function 
is evaluated. Eq. (|A7p . and its analogue for ^ now show 
that Eq. (jASP gives the interaction eneragy to the same 
order as the geometric potential. 

Two vortices of opposite signs that are very near to 
one another (a distance / such that a << I << ro for 
a Gaussian bump, say) cannot tell whether they are in 
curved or flat space. Indeed, the flow fields cancel one 
another outside a range moderately greater than I. li I 
is much less than the curvature scale, the effects of cur- 
vature are negligible. The energy of the vortices should 
therefore be 27rK In - as in fiat space. This conclusion 
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can be checked explicitly in the small aspect ratio approx- 
imation by combining the expressions for the geometric 
interaction and the correction to the Green's function 
Eq. ((X2|) and Eq. ((M|) . (One can take the Ax,Ay-^Q 



limit of 



A^Ax+ArjAy 



with the help of Eqs. 



bh and 



iT]).) As expected, the dependence on the surface pro- 
file h cancels. This consistency check was behind our 
original suspicion about the existence of a geometric in- 
teraction. The energy of two vortices at Ui,U2, without 
the geometric interaction included is — 47r^ii'r(ui, U2). 
This energy differs from the flat space energy even when 
the vortices are very close to one another by a position 
dependent contribution: 

£;„t(ui,U2) w 27rXlnsi2 - ATr^Kg{-{u2 + u^)). (A9) 

This cannot be the correct expression, since as just ar- 
gued, the energy should be the same as in flat space. 
Single particle energies give a simple resolution. If the 
total energy were 

E = ~4:7r^KT{ui,U2) + 27r^Kg{ui)+2TT^Kg{u2), (AlO) 

then all the 5's will canc el when Ui ^ U2. The Gr een's 
function calculations in ( Vitelli and Nelsonl . |2004[ ) and 
the conformal mapping calculations in Section IV.AI show 
that this is actually the correct resolution, and that 

9H = ~^. 



APPENDIX B: The Saddle Surface's Potential 

For the saddle surface with a small aspect ratio (see Eq. 
pop , we may determine the entire geometric potential 
analytically as a function of position. We will only outline 
the procedure here. We would like to evaluate 



Upir) 



1 

2^ 



\n\r-r'\p{r')dx'dy' (Bl) 



when p{x,y) — G{x,y) is the curvature of the surface 
(at the point vertically above {x,y); we are using the 
small-aspect ratio approximation of Appendix [X)l. Some 
thought shows that the curvature given by Eq. (jAip for 
the surface Eq. (|30p takes the form of a polynomial times 

Go — e '0 . We will therefore discuss how to evaluate 
the potential Up for "charge" distributions of the form 



p{x,y) = P{x,y)e ^0 



(B2) 



We start with p = Go; as discussed above, the az- 
imuthal of this distribution symmetry allows its potential 
to be determined by Gauss's Law: 



^Ug„ 



27rr2 



2nr'Go{r')dr' 



(B3) 



This integral is elementary and Ugq can be evaluated by 
one further integration, although this cannot be done in 
closed form. 



Conveniently, the potential due to a distribution of the 
form (|B2p can be determined from the special case of 
p = Gq by differentiation. (Intuitively, derivative charge 
distributions such as dxGo are superpositions of inflnites- 
imally shifted copies of Gq. We can therefore apply su- 
perposition to find potentials for such distributions. This 
is analogous to finding the electric fields of multipoles by 
differenting the monopole field.) To this end, we rewrite 
Eq. (jBip for the special case p = Go as 

f ^ln\A\Goir-A)dA,dAy. (B4) 



Ucoir) 



where (A^,, Aj,) are the components of A = r — r'. It 
follows that 

d^^dl^Uaoir) = - / ^ In |A|a^9^"Go(r - A)dA,dA,. 

(B5) 
The right hand side represents the potential correspond- 
ing to the source in Eq. (|B2p with a special degree k 
polynomial in place of P. This polynomial, obtained by 
multiple differentiations of a Gaussian, is very compli- 
cated, but we will show that polynomials of this specific 
form can be superimposed to give any desired polyno- 
mial (including the degree 8 polynomial appropriate to 
our Gaussian saddle-surface). We will then, in principle, 
be able to express Ug as a superposition of Ugo and its 
derivatives. 

The expansion of the charge distribution of Eq. (|B2p in 
terms of the derivatives of Go can be carried out with the 
help of Fourier integrals. Our goal is to find an expression 
of the form 



(B6) 



pix, y) = P{x, y)e-- -y = Q{d,,dy)e-- 



where P{x,y) is the polynomial appearing in Eq. (|B2p . 
We have to determine a polynomial operator Q{dx, dy) = 
Sn m InmdJ^d^ SO that Eq. (jB6p is true. We have set 
To = 1 for convenience. Applying the Fourier transform 
to both sides of Eq. (jB6p gives 



or (by substituting u = ipx,v = ipy), 



Q{u,v) 



= .-(- 



'^^P{~du,~d,)e^^^^ 



(B8) 



The operator Q{dx,dy) which satisfies Eq. (|B6p can be 
produced by working out the derivatives in this expres- 
sion and replacing u and v by dx and dy. Now the po- 
tential can be worked out using 



UG{r)^Q{dx,dy)UGo{r). 



(B9) 



In fact, multiplying Eq. (jB5P by g„„j and summing 
over n and m shows (with the help of Eq. (|B6p that 
Q{dx,dy)UGo{r) = -^ n\n\S\p{r - A\)dAxdAy. 

Since all derivatives of Ugo can be calculated analyt- 
ically starting from Eq. (|B3p . Eq. (jB9p will yield an 
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analytic expression for [/g, provided we can show that Q 
has no constant term. To show this, we integrate both 
sides of (IB6I) to see that 



7rQ(0,0) = / P(x,2/)e-(^'+^')da; dy ^ G{x,y)dxdy. 

(BIO) 
That is, Q's constant term is proportional to the net 
Gaussian curvature; since the net curvature is zero for 
any surface which flattens out at infinity, Q has no con- 
stant term. 

The potential of the saddle surface can thus be deter- 
mined in closed form by the following procedure: expand 
the curvature to determine the polynomial P. Calcu- 
late Q from (jB8[) . Since (jB10[) guarantees that Q has no 
constant term, we may calculate the geometrical poten- 
tial by differentiating (jB3[) repeatedly. This method is 
not much more practical for human calculations than is 
numerically integrating (jBip by hand. A computer pro- 
gram, like Mathematica (which produced 272 terms), can 
use Eqs. (|B9p and (jBSp to calculate the values rapidly 
and make the graphs shown in Figs. [TT] and [5] There 
is one comprehensible consequence of these calculations: 
the long distance potential is dominated by a quadrupole, 
attracting vortices from some directions and repelling 
them towards others. Hence there are four additional 
local minima outside of the central trap! 



J 



APPENDIX C: Van der Waals Attraction on a Curved 
Surface 



Because the van der Waals force is very short-ranged, 
falling off like -\, one can approximate the integral ex- 
pression for the disjoining pressure n(x) in Eq. (|68p by 
corrections depending only on the local curvature of the 
substrate. The integral is the total van der Waals inter- 
action energy between a point P at x which is above the 
helium film and all the atoms in the substrate: 



n(x) 



/- 



;(|x-x'|)d3x' 



(CI) 



We now choose a simpler coordinate system (see Fig. 
\^ by rotating space so that the tangent plane to the 
substrate at the point of the substrate closest to P be- 
comes horizontal. Let us take the point of tangency to 
be the origin of our new coordinates, (^1,^2,^3) In this 
coordinate system, P is the point (^1,^2,^3) — (0,0,1?) 
(where D is the thickness of the film at P). The ro- 
tated substrate can be described by its height above the 
new "horizontal plane" (an arbitrary plane parallel to the 
tangent plane) using the equation ^3 = /irot(?i,C2)- The 
disjoining pressure is 



n(p) = / / dfid^ 



'l™t(«l:?2) 



rf6^vcf+^2 + (^-e3)2) = - 



6D^ 



ot(O.O) 

dCid^2 I d^M\/ei 



il 



(C2) 



where we have first integrated over the entire region be- 
low the plane ^3 = 0, thereby getting the van der Waals 
interaction between a point and a flat substrate as the 
first term. We then subtract the surplus energy that has 
been included by integrating over the shaded region (see 
Fig. I 



Since the force is short-ranged, we use the quadratic 

approximation to Kot, hrot{^i,^2) = ^rot(0,0) - 3-K1 - 

^K2 where we have assumed the axes to be aligned with 
the principle curvatures. Finally since ki^2-D ^ 1 this 
wedge-shaped region is extremely thin close to the origin 
and the remainder term can therefore be approximated 
by ignoring the dependence of w on ^3 : 




FIG. 32 The disjoining pressure is the interaction between P 
and the substrate, which is below the surface S in the figure. 
This can be calculated by using the result for the interaction 
with an imaginary flat substrate (below the tangent plane 
Stan) and subtracting the very narrow shaded region. 
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diidU^n, + |«2)t.(y^ 



il 



(C3) 
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This integral can be evaluated in polar coordinates: 

rdrd(t) ki ^ 2 , , '«2 2 • 2 
( — r cos H r sm c 



An = a 



(r2+ 152)3 V 2 ■ — - ■ 2 



Kl + K2 



Hit 



^dr 



(r2 + D2)3 



4D2' 



(C4) 



since the mean curvature H is given by ^(ki + K2)- Upon 
combining this expression with the flat substrate result, 
we obtain Eq. ((7^ . 



APPENDIX D: Consumer's Guide to Green's Functions on 
Compact Surfaces 

The ambiguity in the one- vortex energy (Eq. (|155p ) on 
the sphere also implies that there is no particularly nat- 
ural choice of a Green's function on the sphere. With so 
many choices out there, you'll be greatful for this friendly 
guide to help you focus on the important features and 
possible pitfalls of these different functions. 

The first point you need to know is that all of them 
work pretty much just as well, provided they are used 
consistently; one should not use the single-vortex energy 
Eq. ([9T|) designed to work with a different Green's func- 
tion from the one used to calculate the pair interaction 
Eq. (PH)) . The general definition of of a Green's Func- 
tion, broadened from Eq. (J13p . is that it is a symmetric 
function of two points on the deformed sphere satisfying 
the equation 



V2r(x,y) = -^(x,y)+F(x) 



(Dl) 



The only restriction on the function F is that its integral 
over the deformed sphere must equal 1. (Integrating the 
Laplacian on the left shows that there is no solution un- 
less the right-hand side integrates to zero.) The Green's 
function on a sphere, — tt- In — — — has -^ — 5(X, y) as 

*~n I \ \ r 

its Laplacian (Lubenskv and Prostill992i i. Eq. (|D1|) is a 
more versatile vision of what a Green's function should 
be, using a function F in place of the constant. 

That Eqs. ([9l|) and (|90p give the correct net energy 
follows from the result proven in Sec. IV.AI bv conformal 
mapping to the unit sphere: 



(D2) 



i<j 



The interaction potential in this equation 

rc(x,y) = --^ln|A'-3^| 
Ztt 



(D3) 



satisfies 



V^r,(u,u') = e^^f") 



-(5t(u,u') 



^ U^ sphere\yi ,1^ ) 
1 

g2a;(u) 



+ 7Z) 



47r 



(D4) 



In the first step, the scale factor is introduced to com- 
pensate for the change from the reference to the target 
surface. In the second step, the Laplacian of the sphere's 
Green's function — T^ln^^ is substituted. In the third 
step, the (5-function is transformed back to the target 
surface. The last line shows that Tc is a Green's func- 
tion as set out by Eq. (JD1|) . which we call the "con- 
formal Green's function." The F-function that goes with 
this Green's function gets its spatial dependence from the 
conformal factor. 

The single particle potential lo in Eq. (jD2p satisfies 



VuW = Gt(u) -e 



2w(u) 



(D5) 



which follows from Eq. (|105p with the curvature of the 
unit sphere, Gr — 1, substituted. 

Now any Green's function V can be used to solve Pois- 
son's equation for any net-neutral function p on the tar- 
get surface, 



d2u'r(u,u')p(u') = -/9(u). 



(D6) 



This follows from Eq. (|Dip . It can be used to derive Eqs. 
([^I and (|9ip from their special case, the energy derived 
by conformal mapping. We first use the Poisson-like in- 
tegral to "solve" two special cases of Poisson's equation, 
Eqs. (|D4p and Eq. ([D5]) in terms of the arbitrary Green's 
function F. Regarding u' as a constant in the former 
equation, we find that 

g2aj(u") 



Fe(u,u')=F(u,u')- //dV'F(u,u" 

+ /(u') 



47r 



(D7) 



where /(u') is the constant left undetermined by the 
Poisson equation. Since both F and Tc are symmetric 



in u,u', /(u') = -//d2u"F(u',u")^ 



2^(u") 



Ci where 



Ci is a constant. Again applying Eq. (|D6p . this time to 
Eq. (|D5p . implies that 



/"/"d2u"F(u,u")GT(u")-47r/(u) + C2. (D8) 



uj{m) 



(This is not really a solution of the nonlinear Eq. (|D5 
since w still appears on both sides of the equation. 
Rewriting the previous equations implies that 



F(u,u')=F,(u,u')-./(u)-/(u') 
;7G(u)=c^(u)+47r/(u)-C2, 



Ci 



namely that Uq and F are related to oj and Tc according 
to the energy-shuffling transformation Eq. (|155p so that 
the more general expressions of Eqs. (|?1]) and ([50)1 can 
be used in place of Eq. (|D2p to determine the energy. 
The sum of the energies from Eqs. |9T|) . ([90]) is equal 
to the correct energy, Eq. (|D2p up to a constant. The 
arbitrary Green's function can also be used to find the 
flow pattern according to the formula 



N 



huj 



-^ — ' m 



(D9) 
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Calculability 


Neutralizer 


Limit 


Conformal 


+ 


- 


- 


Pair 


+ 


+ 


- 


Standard 


- 


+ 


+ 



TABLE II The advantages and disadvantages of the Green's 
functions, as far as their ease of calculation, simplicity of the 
neutralizing function F, and limiting behavior in case the 
deformed sphere is stretched into a deformed plane. 



There are some advantages and disadvantages of dif- 
ferent choices for F in Eq. (|D1[) . Let us focus on the 
most popular choices. The "standard Green's function" 
is defined with F — ^ {A is the area of the surface) and 
is simply related to the eigenfunctions of the Laplacian, 



(x,y) = ^-*A(x)**A(y). 



(DIO) 



The "pair Green's function" is defined via conformal 
mapping, 

rp(x,y) = Tsphere{X,y) + -^ (l^(x) + c^(y)) (Dll) 

47r 

and incorporates all the single-particle energy into the 
interaction energy, so that Upair = 0. This Green's func- 
tion satisfies the most elegant differential equation. 



V2rp(x,y)=-(5(x,y) 



G(x) 
47r 



(D12) 



Last, the "conformal Green's function" (which was our 
starting point) has F — ^ (see Eq. (|D2p '). 

If you are looking for style in your Green's functions, I 
would choose the pair Green's function. It is easy to cal- 
culate by conformal mapping (Eq. (jDlip ) but it can be 
defined without referring to oj, Eq. (jD12p . much prefer- 
able to the haphazard looking Eq. (jD4p defining the con- 
formal Green's function. The standard Green's function 
is stodgier and does not handle well. The methods for 
finding the standard Green's function, Eq. (jD10|) are 
more limited and, if one wants to use it, the best option 
might be to derive it by using conformal mapping any- 
way: 



rs(x,y) =rc(x,y)- 



1 



At 



re(x,u)+r,(y,u)d2u+C3 



(D13) 
(This equation is derived analogously to Eq. (jD7p .) On 
the other hand, there are always advantages to familiar- 
ity. In particular, in the limit where part of the deformed 
sphere is stretched out to infinity so that it actually be- 
comes a deformed plane, Fj converges to the ordinary 
Green's function of a noncompact surface, since 4 tends 
to zero. For a short summary of all the Green's functions 
features and failings, see Table [TTl 



APPENDIX E: Approximations for Long Surfaces of 
Revolution 

Let us start by determining the conformal map from 
the surface of revolution defined by the equation r = r{z), 
Zs l£ z < Zn, to the unit sphere. We use the coordinates 
0, a introduced in Sec. IV.DI to parameterize the surface; 
a is given by: 



dz 



(El) 



The Cartesian coordinates are x — r(a) cos (j),y 
r{a) sin(/>, z == ^(cr), and hence the metric is 



ds^ = dx'^ + dy^ + dz^ = da"^ + r{(Tf 



(E2) 



If the map {a, (f)) -^ (8, $) to the unit sphere is to be 
conformal, then according to Eq. ((86|) . 

da^ + r{afd(t? = e-'^'^{dQ^ + sin^ Qdcf^). (E3) 

By syrr imetry, $ = (^ and 9 = 9 (g) and is independent of 
(/) (see (jVitelli and Nelsonl . [2003) for the analogous use of 
symmetry on a rotationally symmetric bump on a plane). 
By matching the coefficients of del) and da one finds that 

da 



jT^, or (after integration): 



r 

sin 6 = sech( / 

J (7 



da' . 
r{a')' 



(E4) 



where ctq can be an arbitrary arc length. According to 
Eq. (JE3|, w = ln^ — 



da ' 



= ln 



1 



r(cr) 



sech 



da' 
r{a'y 



(E5) 



To determine the energies and flow patterns on a 
rotationally symmetric surface, we use the "Pair Green's 
function" Eq. (jDlip . the Green's function which incor- 
porates all of the energy into interaction-energy terms. 
This Green's function can be found using Eqs. (jE5P and 
(jE4p ; adding cj at the sites of the vortices to the Green's 
function on the sphere, 

-^ In ^/2[1 - cos 9i cos 82 - sin6i sin62 cos(<i>i - $2)], 
and rearranging, gives 




2r{ai)r{a2) 



[cosh/J^ 



3(01 - 02)] 

(E6) 



The energy of a set of vortices is simple using the 
pair Green's function (see the previous appendix), E = 
J2i<j 47r^ninji^Fpair(ui,Uj). As an example the energy 
of a vortex-antivortex pair at opposite sides of a circle of 
latitude (0i = 02 + vr and ai = a2 = a) is 



£■ = 27r In 



2r{a) 



(E7) 



showing that the energy grows logarithmically with the 
distance between the vortices in this case, as in Eq. (|136p . 
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To prove the azimuthal symmetry of the flows, note 
that accordmg to Eq. (|D9p . the flow velocity at u is 



eflipsoid can be expressed in terms of a in the form: 



n,;/i. 



m 



(E8) 



Now if the vortices are aU far from u, then the integral in 
(|E6P is very large. Since ln(yl + e) fa A + -j for large A, 
the cosine term, the only one which depends on the az- 
imuthal angles, gives exponentially small contributions. 
Therefore the flows can be calculated as in Sec. IV. Dl by 
using the circulation quantization and the approximate 
azimuthal symmetry to determine the flow speeds. Al- 
ternatively, we may calculate the velocity directly from 
Eq. (|E6p with the help of the further approximation that 
cosh a; « ^e'^', yielding 



r(u,uO«-- 



1 

An' 



da' 1_ r(cr)r(o-,) 

r At: a?- 



(E9) 



The flrst term gives the flow pattern of Eq. (|141[) , after 
a brief calculation using the neutrality constraint, while 
the second term, when summed as in Eq. (jE6|l . cancels 
out also by neutrality. 

For a surface (such as an ellipsoid) where the xy-plane 
is a plane of symmetry, our results will simplify if we 
make the choice erg = <^eq in Eq. (|E4p where a^q is the 
arclength corresponding to the equator, at z = 0. In 
this case, the conformal map takes pairs of antipodal 
points on the deformed surface to antipodal points on 
the sphere. (Since antipodal points (cri,0i), (172, ^2) are 
points at opposite ends of a diameter of the surface, (72 = 
Ideq — (Ti, (^2 = TT -I- (/!)i.) If wc cousidcr the interaction 
energy of a pair of antipodal points, we find according to 
Eqs. dMl) and 



(ElO) 



^antipodal c\ y c\ / \ 
= ZTT m ZTTLjiZi) 

K a 



Whether the two vortices are at opposite tips or at oppo- 
site ends of the equator, their image vortices are always 
at the same distance on the unit sphere, so the first term, 
the interaction energy of the images, is a constant. This 
gives another illustration of the folly of making a strict 
separation between intervortex and curvature-vortex in- 
teractions. One would like to think that the growth of the 
energy as the two vortices are separated on an elongated 
surface is due to the attraction between them. But Eq. 
(lElOp shows that it can also be interpreted as resulting 
from the single particle potential lo. 

Let us now turn to the problem of describing the equi- 
librium positions of a pair of vortices on a rotating el- 
lipsoid. Both the transitions at fia and fib can be un- 
derstood only with a more accurate version of the force 
than the band- force approximation, Eq. (|139p . The er- 
ror in the approximation is important when the vortices 
are near the poles of the ellipsoid or, as just illustrated, 
near each other. We will assume that the aspect ratio of 
the ellipsoid a = -p is very large. The equation for the 



r = R\ I 



Y^R^' 



(Ell) 



The energy of a vortex-antivortex pair according to Eq. 
(gni) is 



E =^ Er 



En. 



(E12) 



Here the energy of the flow pattern, or "resting energy," 
is the energy of the vortices on a stationary ellipsoid, 
Brest = ~^'^^KTpair{ai,( j)i;a2 ,4>2)- The "rotation en- 
ergy" Eq is given by Eq. (I144p . Both energies are func- 
tions of a single variable, the distance s between the two 
vortices along the surface, if we assume that the vortices 
are at iai,(f>i) = {aeq - f,0) and (0-2, ^2) = (ceq + f,0). 
These relationships assume that the vortices are situated 
symmetrically about the xy-plane; the vortices will have 
equal azimuthal angles in order to minimize the energy 
of the flow pattern. 

The equilibrium position of a pair of vortices is deter- 
mined by balancing the rotational force and resting force 
acting on one of them. The resting force Frest on vortex 1 
is derived from the kinetic energy of the flow pattern and 
is positive since it pulls the vortices toward each other in 
order to decrease the width of the band of moving fluid 
between them. The rotational force Frot is negative since 
it pulls the vortices toward the poles of the ellipsoid (and 
away from each other) in order to increase the total angu- 
lar momentum of the flow. At equilibrium the rotational 
and resting forces on the vortices balance, as we can see 
by differentiating Eq. (|ET2)) to obtain - ^^ 
or equivalently 



ds 



dEn 
ds 



Frest (s) 



-Fn{s) 



(E13) 



where Frest and Fq are the resting and rotational forces 
on the vortices. (A short calculation shows that the force 
on one of the vortices — 4^ is equal to 4^ , since the en- 

dai ^ ds ' 

ergy change produced by moving one vortex an inflnites- 
imal distance is the same as the energy change produced 
by moving both vortices half the distance.) The equilib- 
rium positions can be found by graphing Frest and — -Fh 
as in Fig. [33] and finding the intersection points. The 
exact expression for the resting force can be found by 
differentiating Eq. (jE6p to obtain 



^ rest — 



t:K 
r(cri) 



Dth[ 



'"■=' da' 



^^"'^ y^^(¥)^ 



(E14) 
The rotational force is given exactly by Eq. (|145p . 

Fig. [33] shows the resting force and minus the rota- 
tional force on one of the vortices for Vt = Vti,, Q, = Q.a 
and for an intermediate value of the frequency. 

Stability of the equilibria illustrated in Fig. [33] can 
be determined by considering the direction in which the 
resting force curve crosses the rotational force curve. The 
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FIG. 33 The resting force (dashed hne) and minus the rota- 
tional force (solid lines) between two vortices on an ellipsoid 
of aspect ratio 3.5. The rotational forces are illustrated for Qa 
(top curve), Sli, (the bottom curve) and an intermediate value 
of SI. For smaller values of SI, the attraction between the vor- 
tices always overcomes the rotational force, causing them to 
annihilate. For larger values of Q, the rotational force over- 
comes the attraction, causing the vortices to move to opposite 
poles. 



middle point of the three equilibrium points at the in- 
termediate frequency is a stable equilibrium because the 
resting force curve crosses the rotational force curve from 
bottom to top. This implies that if the vortices fluctu- 
ate away from each other (increasing s) , then the resting 
force becomes stronger than the rotational force and pulls 
them back together. 

Let us consider how the stable equilibrium disappears 
at Slfc. As SI is lowered the stable and unstable equi- 
librium come together and then "annihilate" when the 
rotation-force curve detaches from the resting force curve, 
as illustrated by the lowest curve in Fig. [33l which cor- 
responds to fl = fib- Since the rotation force and resting 
force curves are tangent at Qb, the frequency Qb and sep- 
aration of the vortices Sb at this transition point can be 
determined by solving Eq. (|E13p simultaneously with 



F^,,t{sb) ^ ~F[,isb). 



(E15) 



When a 3> 1, we will be able to avoid solving simultane- 
ous equations since the value of Q,b is already determined 
by Eq. (|148p . Using this result, we will be able to solve 
Eq. (|ET5l) for Sb- 

The simple band approximation to the force, Eq. 
(|139p . suggests that the vortices move continuously to- 
ward one another as f2 is decreased, annihilating at the 
equator. Substituting the expression for the critical fre- 
quency that is implied by the band model, Eq. (|148p . into 
Eq. (|147p in fact implies that Sb = 0, which is incorrect. 
The band approximation fails because it implies that the 
force between the vortices decreases monotonically as the 



vortices approach one another. Though in conflict with 
our intuition from the plane, this result is correct over 
the large middle range of the resting force curve in Fig. 
[551 As the rotational conflnement weakens, the vortices 
get closer together, and the resting force weakens too, 
preserving the equilibrium. However the resting force 
starts increasing strongly as the vortices approach one 
another, because the vortices start to feel one anothers' 
asymmetric flow fields. This force will certainly overcome 
the rotational force when the rotational conflnement de- 
creases further. (Actually, Eq. (JE15P implies that Sb does 
not correspond exactly to the maximum of Frest because 
the rotational conflnement is not a constant force fleld.) 
We can derive the corrections to the force from Eq. 
(|E14p ; if a is large, we may neglect the second term and 
assume that 



da' 
r(cr') 



R 



(E16) 



since the radial profile of the ellipsoid, Eq. (jElip . is 
slowly varying. We then obtain the approximation that 
is valid when the vortices are close (compared to ^, the 
characteristic distance for variation of the radius) . 



ttK , , s 

Frest « -r^COth -— , 

r(tTi) 2R 



(E17) 



Notice that the force diverges as ^2^ when the vortices 
are close together (as in the plane) and approaches Eq. 
(|139p exponentially fast as the vortices move apart; this 
generalizes the band model approximation to the case 
where the two vortices may be close together. As we will 
see, for a large value oi a, Sb » R at the moment when 
the vortices annihilate. We therefore simplify Eq. (|E17|) 
by making another approximation, coth x « 1 -|- 2e~^^. 
Then an approximate version of Eq. (jElSp that is derived 
from Eqs. (|E17p . pi5)) reads 



ttK 1 dr 
~R^d^^'"' 



2ttK- 



i?2 



TrQbfiPs dr 
m aai 



(E18) 



The first term describes the decrease of the resting force 
due to the variation in r(z). The second term results from 
the exponentially decaying portions of the flow fields. 
(We are replacing ^(cti) by R whenever that is accurate 
enough since the width of the ellipsoid is slowly varying. 
Of course, the slow variation of r(cri) is important in 
some terms; the resting force initially decreases as s de- 
creases because the band approximation to the force de- 
creases with increasing circumference.) Using Eq. (jl48p 
for Slfc in Eq. (|E18P gives 



dr 
dui 



= 2e^~ 



(E19) 



In order to evaluate the left-hand side, we note that ui — 
o'eq ~ f ~ '^eq — z aside from terms of order -ij since the 
sides of the ellipsoid are nearly vertical near the equator. 
Therefore Eq. (|ETT|) implies that ^ = 3^. 
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Rearranging Eq. IE19I now gives 



Sh — R\'a- 



Sb 



(E20) 



which can be solved by substituting it into itself. The 
first iteration gives 



Sb 



RlnAa'^ - Rlnln 



ia'^R 

Sb 



(E21) 



Since the second term has two logarithms in it, it is 
smaller than the first in the limit where a — > oo, so finally 



Zb ~ i?lna, 



(E22) 



(since Zf,, the distance from a vortex to the equatorial 
plane, is approximately half the distance between the 
vortices). We have justified Eq. (|149p . Two iterations 
of Eq. (|E20l) give Zb = R\na - ^\n^-^; the error for 
this approximation actually approaches for large a. 
The exact result can be found by computer, but the 
approximate result is reasonable even at a = 5, where 
§ = 1.8 wln5 = 1.6. 

The height Zb depends only logarithmically on a be- 
cause the extra short-distance vortex-vortex interaction 
decays exponentially and would not be strong enough to 
pull the vortices together if Zb were very large. (Check 
this by substituting our final result, Eq. (|E22p . into Eq. 
(jElSp . All the terms, the ones from the band model as 
well as the exponential correction, have the same basic 
dependence on a.) To see that the approximations we 
have made are valid, one has to calculate F^^^^{s) from 
the exact expression Eq. (|E14|1 . The resulting expres- 
sion can be simplified by dropping various terms, which 
mostly have a relative size of -^ and (i2-«)2. ^.j^g reason 

is that ^ — is^ so the vortices are proportionally very 
close to the equator, and again r{z) can be replaced by 
R. (This also justifies the approximation in Eq. (|E16|) 
where the integrand is replaced by a constant.) One par- 
ticularly large term, resulting from the second term of 
Eq. (jEUp . has been neglected in Eq. (|E18|1 . but the 

K-Trr" (a-i) 



2r(ai 



is still of relative order -r-^ . 

In a 



neglected term, 

Now we turn to the critical frequency fla where the 
vortices move to the poles. The band model also requires 
a correction in order for this transition to be described 
correctly. In fact Eq. (|147[) would imply that the vortices 
never exactly reach the tips of the ellipsoid even as 51 — > 
oo, and thus fia = cx3. In fact, the band force on the 
left hand side of Eq. p46p , which approaches infinity at 
the poles, cannot be balanced by the rotational force at 
a finite frequency. Of course, the exact force approaches 
zero rather than infinity at the pole (see Fig. [55]) . The 
value of Oa may be derived from the condition that the 
actual resting force curve and the rotation force curve 
have to be tangent at the origin, as for the uppermost 
(fi = fla) curve in Fig. [331 We therefore have to find 
when s = 2cre5 satisfies Eq. (jElSp . Linearizing Eq. (|E14p 




FIG. 34 Illustration of the combined rotational and resting 



energies near Qa, for H = 3.5R and 



0,1.2, 3.1, 4.9i?" 



The first of these is the uppermost curve. The third curve 
corresponds to fi = Qa- Although the second curve looks 
practically flat on this scale, it has a curvature of about ^^ 
at its off-center minimum. 



near cri = to find the derivative of the force implies that 
fla is given by 



mfln a^K? 



1 rj rR/ da 

-I p~^ Jo y'd7 

^ R^ 



(E23) 



where k is the curvature at the tip of the ellipsoid. The 

2 

critical frequency is larger than the result fi^, derived 
in Section IIII.Bi for a bump with the same curvature 
because the rotational confinement must overcome the 
mutual attraction of the vortices as well as the repulsion 
of the vortices from the curvature. For an ellipsoid with 
a large value of a, the correction term is unimportant, so 



fla 



y^n 



AmR^ 



(E24) 



The transition can be visualized using the energy curves 
illustrated in Fig. I34( where the local minimum of the en- 
ergy function moves away from the axis as 57 is decreased 
through Q,a- 

There is an aspect ratio Uc below which there are no 
off-center local minima, for any rotation speed. That is, 
when the angular velocity is decreased enough, a vortex- 
antivortex pair initially at the poles immediately moves 
to the equator and annihilates. This situation is illus- 
trated for a sphere in Fig. [351 The value of ac can be 
determined numerically, and is 1.33. One simply graphs 
the total energy at fi = f^a (as given by the exact expres- 
sion, Eq. (|E23[) ) and checks whether there is an energy 
barrier or not. At O > O^, a pair of vortices at the 
poles will be stable. If there is no barrier, as in Fig. [351 
slightly decreasing fi will cause these vortices to leave 
the poles and annihilate each other. If there is a barrier, 
as in Fig. [331 slightly decreasing 57 will create an off- 
center local minimum. This can be seen from the energy 
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FIG. 35 The combined rotational and kinetic energy for vor- 
tices on a sphere, where there are no stable off-axis posi- 
tions. The energy is graphed as a function of cti = 2aeq — (72 
in units of 7?. From the top, the curves correspond to 
fl — 0, .4, .52, 1.5i?~^. The third of these corresponds to 



curves: there is a local maximum at the origin, because 
fl < Qa, and there is also a local maximum at the top of 
the barrier. Therefore the vortex can find a local mini- 
mum somewhere in between. 



APPENDIX F: Derivations of Bounds Valid Even for Strong 
Distortions 

The results of Sec. IVII can be derived from theorems on 
"univalent" analytic functions. We will state these the- 
orems here and d erive the limi ts on the geometric force 
from them. (See ( Rudinl . [l987t) for the proofs.) An ab- 
stract example of the type of question these theorems 
address is the following. Let f{t) be an analytic function 
defined by the following series: 



/(t) = t + a2t^ + aaf^ + 



(Fl) 



Suppose this series converges out to radius 1, at least. If 
one of the coefficients, maybe ag, is much larger than the 
rest, then the function is dominated by the t^ behavior, 
and most points in the range of the function will occur 
six times as values of the function. Therefore, if one is 
looking for a univalent function (a function which is one- 
to-one inside the unit circle) then there will be upper 
limits on the sizes of the a„ . A challenging mathematical 
problem is "What are the maximum sizes for the a„'s?" 
The answer (proved by De Branges) is that |a„| < n, and 



that the function 



IT^tp 



attains the maximum value for 



every Taylor series coefficient simultaneously. To find the 
upper bound on the vortex force in a flat disk, we will 
use only the bound 



which has a simpler proof teudinl . UQSTt) . Note that the 
conditions of this theorem do not require that the func- 
tion remains one-to-one outside the unit circle. For ex- 
ample, the fmiction t + At^ satisfies the conditions of the 
theorem although it takes on the value zero at t = 0, —10. 
The analyticity of f{t) is allowed to break down as well 
beyond a radius of 1. 

Similar problems can be stated for functions g(t) de- 
fined outside of the unit circle, with expansions of the 
form 



g{t) ^ i + ^ + -i 



(F3) 



To make the predictions about the quadrupole force due 
to a bump in a plane we will use the Area Theorem 
(JRudinl . Il987l ) which states that, if g is one-to-one and 
analytic outside the unit circle, and 



then 






ai < 1. 



(F4) 



(F5) 



To prove Eq. (|167p . one just realizes that the assump- 
tion means that a part of the surface has the same ge- 
ometry as a radius R disk in the plane with a vortex at 
the center. We can introduce a coordinate system on this 
portion of the surface by introducing Cartesian coordi- 
nates u, V (with w = u + iv) on the disk in the plane, 
and then mapping these coordinates isometrically to the 
surface. This mapping is different from the conformal 
mapping C used to calculate vortex energies. To relate 
them, let Z ^ X + iy where X, y are the coordinates 
of the conformal image of the surface. Then Eq. ([5^ 
takes the form du^ + dv^ = er^'^idX^ + dy'^) and it fol- 
lows that 2{w) is a conformal map from a portion of the 
plane to itself, hence an analytic function of w on the 
circle of radius R (say Z — ciw + C2W^ + ...). Further- 
more, rewriting the expression for the scaling of lengths 



as \dw\^ 



\dZf, we see that 



In I 



aw 



(F6) 



We now define 



f{t) = 



Z{Rt) - Z(0) 
Rci ■ 



Then / is a one-to-one analytic function on the unit circle 
(which is scaled by i — > i?t into the radius R circle). Since 
/(O) = 0,/'(0) = 1, / has the form of Eq. dFT]), so Eq. 
(jF2| implies 



Cl 



(F7) 



02 < 2 



(F2) 



Now the force on the vortex is 7rXV(jj(0) which can be 
expressed in terms of the coefficients of Z's Taylor series 
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FIG. 36 Construction of the largest circle centered at a point 
on a cone with 6 < n. The cone is cut open and flattened so 
that the center is on the bisector of the angle. 



by means of Eq. ^■. F = 2t:K{^^, -3ff ). The upper 
bound, Eq. plT)) . follows from Eq. ' (|F7)) . 

Let us now see whether the bound just proven can be 
improved at all; i.e., whether the ratio of the force on 
a singly quantized vortex to -„■ can ever be as big as 
47r. For example, for vortices on cones, the ratio of the 
force to ^ is maximal in the limit where the cone angle 
6* ^ 0. To find this ratio, we must take R to be the 
radius of a disk centered at the vortex which is flat and 
non-self-intersecting; taking R = Rmax, the radius of a 
disk which is as large as possible, maximizes the ratio 
we are interested in. The radius Rmax can be found by 
imagining the disk expanding out from the vortex. If 
6 > TT, Rmax—D bccausc the first calamity that befalls 
the disk as it expands is that it starts overlapping the 
cone's apex. But if 6* < tt, then the disk overlaps itself 
before this as one can see on the unfolded version of the 
cone illustrated in Fig. 1361 Some simple trigonometry 
shows that Rmax — D sin | . Eq. (|170p shows that for 

small 9, -Ljj > tt^, which is a little less than 47r. 

There is a surface that saturates the original bound, 
though; this surface is illustrated in Fig. [37l The surface 
is obtained by folding a disk in half and sealing it shut ex- 
cept for a very small opening at one end of the diameter. 
This opening is then connected to an infinite plane. The 
top part of the substrate is a semi-circular slab with the 
superfluid layer laminating both sides so that the helium 
spreads out to the plane. The topology of the helium film 
is still that of a plane. A vortex placed at the center of 
the disk, B, saturates the bound; it is attracted by the 
negative curvature of the neck joining the plane and the 





FIG. 37 (a) The surface which contains an isometric disk of 
radius R and has the maximum geometrical force. A semi- 
circle (with films on both sides) is connected by a neck A to 
a plane. If a single vortex is placed at B, the force on the 
vortex approaches ^^ as the edges of the surface becomes 
sharper, (b) An unfolded image of the flow pattern. 



disk and is repelled by the positive curvature at the top 
of the fold. 

To show that this surface (the "calzone surface") sat- 
urates the bound, we will find the force on the vortex 
using conformal mapping. Instead of mapping the entire 
surface to a reference plane, we can just map the folded 
disk portion of the surface. The flow patterns on the two 
portions of the surface are uncorrelated when the neck be- 
comes infinitely small, aside from requirements imposed 
by the circulation's invariance. The circulation around 
any curve on the plane enclosing the neck will be — be- 
cause the vortex in the disk region is inside it, and the 
flow pattern on the planar base will not be sensitive to 
the location of this vortex because the neck is so narrow. 
It will consist of a set of concentric circles, representing 
a flow whose energy is independent of the position of the 
vortex. The force does not depend on this portion of the 
flow, so the two portions may be dismantled at the neck. 
As illustrated in Fig. [371 the neck now turns into the 
core of a second vortex, at point A of the folded disk. 
The folded disk now has the topology of a sphere, sat- 
isfying the neutrality condition because the two vortices 
are equal and opposite. The map Z = ir-Z)'^ *^^ ^^^ 
radius R disk can be used to relate the folded disk to a 
reference plane, since the points on the circle which fold 
together, w = Re^^'^, both map to the same point of the 
real axis in the plane. Since the vortex at B maps to 
infinity, the force on the vortex at A can be calculated 
from the geometric potential alone (without any inter- 
action terms), giving t^K 2„L} = ^^. Also the flow 
pattern illustrated in the figure can be found by map- 
ping the concentric circles centered around the origin in 
the Z plane to the disk using the function w{Z). 

The result Eq. ()172|) about the long range force due 
to a bump contained inside of a radius R but with an 
arbitrary height and arbitrary curvatures follows (by an 
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argument similar to the one used for the first inequahty) 
from Eq. (|F10p . The conformal mapping takes the fiat 
part of the surface (a plane with a radius R hole parame- 
terized by the complex variable w) in a one-to-one fashion 
to a reference plane with a hole of some distorted shape. 
As above, this function is analytic and cj = In |^|. By 
rescaling one can ensure that Z ^ w at infinity. Applying 
the area theorem to g{t) = -^Z{Rt) shows that 



Ziw) 



R'^ 



(F8) 



where ai < 1. Expand E — ^ ^^^^^ ^^ 
find the large distance form of the energy 

Er^ -ttKR^^^; 
w 



^ Sft In 4§ for large w to 



(F9) 



it follows that ^2 = '!TKR^\ai\ and 72 = arg(ai) + tt in 
Eq. (|17ip and the bound on the quadrupole moment 
M2 < ttKR^ follows from the bound on ai. 

Now we can also ask what type of bump maximizes the 
quadrupole moment. It turns out that the value ttKR^ 
cannot be attained by any surface which is flat outside a 
circle of radius R. There is a surface which consists of a 
bump surrounded by a surface isometric but not congru- 
ent to K , the plane with a circle of radius R removed. 
This surface is gotten from K by sealing opposite sides 
of the circle together to make a mountain ridge. 

The reason this surface has the biggest quadrupole mo- 
ment is because its conformal mapping to the plane is the 
function that maximizes ai . According to the area theo- 
rem the only one-to-one analytic function on K for which 
ai = 1 is 



Z{w) 



w 



(FIO) 



This function maps both points i?e**'^ to the same point 
in the reference plane, so any flow pattern on the target 
plane will still be continuous when the two edges of K 
are sealed. 

This quadrupole-maximizing surface does not contain 
a flat copy of K. Hence an open question is to find the 
largest value of fi2 for a bump in a plane which is actually 
flat outside a radius of i?, as well as the shape of the bump 
which has this maximum quadrupole moment. 



References 

Anibegaokar, V., B. I. Halperin, D. R. Nelson, and E. D. 

Siggia, 1978, Phys. Rev. Lett. 40(12), 783. 
Andelman, D., J. F. Joanny, and M. O. Robbins, 1988, Euro- 

phys. Lett. 7, 731. 
Batchelor, G. K., 2001, An introduction to fluid dynamics 

(Cambridge University Press, Cambridge). 
Bausch, A. R., M. J. Bowick, A. Cacciuto, A. D. Dinsmore, 

M. F. Hsu, D. R. Nelson, M. G. Nikolaides, A. Travesset, 

and D. A. Weitz, 2003, Science 299, 1716. 
Bowick, M., D. R. Nelson, and A. Travesset, 2000, Physical 

Review B 62, 8738. 



Campbell, L. J., and R. M. Ziff, 1979, Physical Review B 20, 

1886. 
Caspar, D., and A. King, 1963, in Basic Mechanisms in Ani- 
mal Virus Biology (Cold Spring Harbor Symposia on Quan- 
titative Biology XXVII). 
David, F., 1989, in Statistical Mectianics of Membranes and 

Surfaces, edited by T. P. D. R. Nelson and S. Weinberg 

(World Scientific, Singapore). 
DeVries, G. A., M. Brunnbauer, Y. Hu, A. M. Jackson, 

B. Long, B. T. Neltner, O. Uzun, B. H. Wunsch, and 

F. Stellacci, 2007, Science 315, 358. 
Dinsmore, A. D., M. F. Hsu, M. G. Nikolaides, M. Marquez, 

A. R. Bausch, and D. A. Weitz, 2002, Science 298, 1006. 
Dubrovin, B. A., A. Fomenko, and S. P. Novikov, 1992, 

Modem Geometry- Methods and Applications, volume 1 

(Springer- Ver lag. New York). 
Earnshaw, S., 1842, Proc. Trans. Cambridge Phil. Soc. 7, 97. 
Fernandez-Nieves, A., V. Vitelli, A. Utada, D. R. Link, D. R. 

Nelson, and D. A. Weitz, 2007, Phys. Rev. Lett. 99, 157801. 
Glaberson, W. I., and R. J. Donnelly, 1966, Physical Review 

141, 208. 
Guyon, E., J. P. Hulin, L. Petit, and C. D. Mitescu, 2001, 

Physical Hydrodynamics (Oxford University press, Ox- 
ford). 
Halperin, B. I. Halperin, private communication. 
Hashimoto, N., R. Goto, H. Yano, K. Obara, O. Ishikawa, 

and T. Hata, 2007, Physical Review B 76, 020504. 
Hexemer, A., V. Vitelli, E. J. Kramer, and G. H. Fredrickson, 

2007, Phys. Rev. E 76, 051604. 
Hexemer, A. Hexemer, PhD thesis. University of California 

Santa Barbara, 2006. 
Hyde, S., K. Larsson, Z. Blum, T. Landh, S. Lidin, B. W. 

Ninham, and S. Andersson, 1997, The Language of Shape 

(Elsevier, Amsterdam). 
Israelachvili, J. N., 1985, Intermolecular and Surface Forces 

(Academic Press, New York). 
Jeans, J. H., 1927, The Mathematical Theory of Electricity 

and Magnetism (Cambridge University Press, Cambridge), 

fifth edition edition. 
Kamien, R., 2002, Rev. Mod. Phys. 74, 953. 
E. Katifori, E. Cerda, D. R. Nelson and J. Dumais (to be 

published). 
Lubensky, T. C, and J. Prost, 1992, Journal de Physique H 

2, 371. 
Nelson, D. R., 2002, Defects and Geometry in Condensed Mat- 
ter Physics (Cambridge University Press, Cambridge). 
Nelson, D. R., 2002, NanoLetters 2, 1125. 
Nelson, D. R., and J. M. Kosterlitz, 1977, Physical Review 

Letters 39, 1201. 
Park, J.-M., and T. C. Lubensky, 1996, Phys. Rev. E 53, 

2648. 
Pearl, J., 1964, Applied Physics Letters 5, 65. 
Polyakov, A. M., 1981, Physics Letters B 103, 207. 
Polyakov, A. M., 1987, Gauge Fields and Strings (Harwood 

Academic Publishers, Chur, Switzerland). 
Quere, D., J. M. di Meglio, and F. Brochard-Wyart, 1989, 

Europhys. Lett. 10, 335. 
Rudin, W., 1987, Real and Complex Analysis (McGraw-Hill 

Company, New York), third edition edition. 
Saffman, P. G., 1992, Vortex Dynamics (Cambridge Univer- 
sity Press, Cambridge). 
Santangelo, C, V. Vitelli, R. D. Kamien, and D. R. Nelson, 

2007, Phys. Rev. Lett. 99, 017801. 
Schwarz, K. W., 1981, Physical Review Letters 47, 251. 



50 



Schwarz, K. W., 1990, Physical Review Letters 64, 1130. 

Scott, W. T., 1959, American Journal of Physics 27, 418. 

Tilley, D. R., and J. Tilley, 1990, Superfluidity and Supercon- 
ductivity (Institute of Physics Pubhshing, Bristol), third 
edition. 

Tinkham, M., 1996, Superconductivity (McGraw Hill, New 
York). 

Trugman, Stuart Trugman, private communication. 

Vinen, W. F., 1963, in Proceedings of the International School 
of Physics "Enrico Fermi, " Course XXI: Liquid Helium 
(Academic Press, New York). 

Vinen, W. F., 1969, in Superconductivity, edited by R. D. 



Parks (Dekker, New York), volume 2. 
Vitelh, v., J. B. Lucks, and D. R. Nelson, 2006, Proceedings 

of the National Academy of Sciences 103, 12323. 
Vitelh, v., and D. R. Nelson, 2004, Phys. Rev. E 70, 051105. 
Vitelh, v., and D. R. Nelson, 2006, Phys. Rev. E 74, 021711. 
Vitelh, v., and A. M. Turner, 2004, Physical Review Letters 

93, 215301. 
Voll, P., N. apRoberts Warren, and R. J. Zieve, 2006, cond- 

mat/0601395. 
Yarmchuk, E. J., and R. Packard, 1982, Annals of Physics 46, 

479. 



